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A    USER'S    GUIDE 
TO    THE    OA36  60   APL   WORKSPACE 

by 

F.  Russell  Richards 
Naval  Postgraduate  School 
Monterey,  CA   9  3940 

I .   INTRODUCTION 

This  report  provides  documentation  for  the  OA36  6  0  APL 
workspace  that  was  developed  during  the  summer  and  fall  quarters 
of  1977  in  conjunction  with  the  offering  of  the  OA3660  course, 
Data  Analysis,  at  the  Naval  Postgraduate  School.   The  OA3660 
workspace  contains  APL  functions,  variables,  and  data  sets  that 
are  useful  for  interactive  exploratory  data  analysis.   The  APL 
functions  include  selected  procedures  from  the  STAT101  public 
library  [6],  and  from  Donald  R.  McNeil's  textbook,  Interactive 
Data  Analysis:   A  Practical  Primer  [5]  (with  some  modifications) 
The  OA3660  workspace  also  includes  functions  written  by  myself 
and  students  in  my  class  and  various  utility  functions  extracted 
from  other  public  library  workspaces.   Documentation  for  each 
function  is  contained  in  this  report  and  in  the  OA3660  Public 
Library  Workspace. 


The  data  sets  include  data  contained  in  examples  and  exercises 
in  the  above  mentioned  book,  data  from  John  Tukey ' s  Exploratory 
Data  Analysis  [8],  data  generated  from  test  scores  in  the 
OA3660  class,   and  data  extracted  from  a  few  other  textbooks. 
The  data  are  cross  referenced  with  the  sources  in  a  variable 
named  DATAMAP. 

This  report  is  intended  for  the  user  who  has  some  experi- 
ence with  APL  and  who  is  taking  a  course  in  Data  Analysis  or  one 
who  is  already  familiar  with  the  basic  techniques  of  interactive 
exploratory  data  analysis  as  described  in  Tukey  [8],  McNeil  [5], 
and  Mosteller  and  Tukey  [7] .   This  report  will  not  attempt  to 
explain  the  analysis  techniques  themselves,  nor  will  it  discuss 
detailed  APL  concepts.   It  will  give  a  brief  description  of  the 
construction  of  data  arrays  and  transformations  of  data  since 
those  are  key  elements  in  interactive  exploratory  data  analysis. 
It  also  includes  a  brief  discussion  of  basic  APL  commands  such 
as  logging  into  the  system,  loading  workspaces,  etc.   For  further 
details  on  APL  and  the  NPS  time-sharing  system,  CP/CMS,  the 
reader  is  referred  to  the  Naval  Postgraduate  School  Technical 
Note  No.  0141-33,  APL (CMS) -An  Introduction  [1];  the  IBM  report 
#GH20-0906-l,  APL\360-OS  and  APL\360-DOS  User's  Manual  [2]; 
APL-An  Interactive  Programming  Language  [3]  by  Gilman  and  Rose; 
and  APL  Programming  and  Computer  Techniques  [4]  by  Katzan. 


Each  function  contained  in  the  OA3660  workspace  is 
discussed;  the  function  syntax  is  given;  the  function 
parameters  (if  any)  are  described;  APL  listings  of  the  functions 
are  given;  and  at  least  one  example  is  included  to  illustrate 
the  use  of  each  function .   An  attempt  has  been  made  to  assure 
that  every  function  in  the  workspace  is  completely  debugged. 
However,  should  any  problems  be  experienced,  the  user  is  urged 
to  notify  me  of  the  problems.   Users  are  also  encouraged  to 
submit  to  me  interesting  data  and  APL  functions  that  are  use- 
ful for  interactive  exploratory  data  analysis. 


II.   APL  PRIMER 

This  chapter  is  written  primarily  for  the  user  who  is 
not  familiar  with  the  APL  programming  language.   We  attempt  to 
provide  only  that  material  that  the  user  needs  to  know  about 
APL  to  use  successfully  the  OA3660  APL  Public  Library.   There- 
fore, we  discuss  in  this  chapter  login  procedures,  basic  APL 
workspace  management,  error  recovery  procedures,  creation  and 
storage  of  data,  data  transformations,  function  syntax,  and 
logout  procedures.   A  user  conversant  with  APL  should  proceed 
to  Chapter  III. 

A  good  understanding  of  APL  would  enhance  the  user's 
facility  for  working  with  data  arrays  and  performing  interactive 
exploratory  data  analysis.   Therefore,  we  encourage  th  user  to 
seek  out  more  detailed  information  on  APL.   References  [1,2,3,4] 
are  all  recommended. 

A.   Getting  Started:   Login 

In  order  to  use  the  OA3660  APL  Public  Library,  the 
user  must  have  access  to  a  computer  terminal  that  has  APL 
capability  linked  to  the  Naval  Postgraduate  School  IBM  360/6  7 
computer.   There  are  many  different  types  of  remote  terminals 
available  for  use,   and  each  has  its  unique  features  of  oper- 
ation.  Therefore,  the  user  should  check  out  the  operating 


* 

The  IBM  2  741  terminals  with  a  special  APL  typing  ball,  the 
Intertec  terminals,  and  some  of  the  CRT  terminals  have  APL 
capability.   In  addition,  there  are  several  school  owned 
portable  terminals  with  APL  capability. 


instructions  for  each  terminal  to  augment  the  general  instructions 
given  here . 

If  the  terminal  is  wired  directly  to  the  computer,  the 
user  need  only  turn  it  on  to  access  the  NPS  CP/CMS  Time  Sharing 
System.   If  not,  the  user  must  link  to  the  computer  via  an 
acoustic  hookup.   First,  turn  the  terminal  on  to  the  correct 

settings.   Then  dial  the  appropriate  telephone  number  for 

* 
connection  to  CP/CMS.    When  the  shrill  audible  tone  is  received, 

place  the  telephone  receiver  into  the  acoustic  coupler  which 

is  connected  to  the  terminal.   If  all  of  the  terminal  settings 

are  correct,  and  if  the  time  sharing  system  is  in  operation, 

the  user  should  receive  the  message  "CP-67  ONLINE".   At  this 

point,  the  user  must  log  onto  the  system  by  typing  (with  the 

non-APL  character  set) : 

LOGIN  XXXXPYY  (CPl 


where   XXXX   is  the  user's  identification  number  assigned  by 

the  computer  center,  P   indicates  a  private  user  ID  (type   G 

if  a  general  user  without  private  disk  space) ,  YY   is  the  terminal 

number,  and  (CPJ  indicates  a  carriage  return.    If  the  identifi- 


cation number  is  valid,  the  system  will  request  the  user's  password, 


* 

The  telephone  number  and  the  terminal  settings  depend  on  the 

terminal  being  used.   For  ASCII  terminals  (nearly  all  terminals 
at  NPS  other  than  the  IBM  2741  terminals)  dial  either  x2611 
(with  speed  setting  at  110)  or  x3025  (with  speed  setting  at  300) 
For  EBCDIC  terminals  (IBM  2741) ,  dial  x2701. 

On  some  terminals  the  user  must  depress  simultaneously  the  keys 
CONTROL  and  S  in  place  of  a  carriage  return  even  though  there 
may  be  a  key  so  labelled.   We  will  use  (CR)  to  indicate  a 
carriage  return  for  all  terminals. 
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The  assigned  password  should  be  typed  followed  by  a  carriage 
return : 


PASSWORD  (CR, 

If  the  password  is  incorrect,  the  system  will  reject  the  login 
and  ask  the  user  to  start  over.   If  the  password  is  correct, 
the  system  will  request  the  user's  four-digit  project  number 
to  be  followed  by  a  four-character  cost  center  code.   The 
project  number  must  be  assigned  by  the  computer  center.   The 
cost  center  code  is  the  user's  section  identifier  or  department 
code . 


PPPPCCCC  (CR 

If  the  project  number  is  acceptable  the  system  will  respond  with 
a  ready  message  indicating  that  the  user  is  in  the  CMS  subsystem 
of  the  CP/CMS  Time  Sharing  System.   Otherwise,  the  user  must 
repeat  the  entire  procedure. 

Once  the  user  has  entered  CMS  he  should  switch  on  the 
APL  character  set*  and  type 


APL  (CR 


* 

Switching  to  the  APL  character  set  may  consist  of  flipping  a 

switch,  changing  a  type  ball  or  impact  print  set,  or  issuing 
a  sequence  of  program  instructions  depending  on  the  type  of 
terminal . 


to  enter  the  APL  subsystem.  The  system  will  then  respond  with 
the  message: 

A*P*L\N*P*S 

LIBRARY  DOCUMENTATION  SYSTEM...    ) LOAD  1  LIBDOC   TYPE  DESCRIBE.* 
The  user  is  now  in  APL,  and  all  of  its  powerful  features 
are  available  to  him.   He  may  use  APL  in  the  calculator  mode 
somewhat  as  he  would  use  a  hand  held  calculator;  he  may  define 
his  own  functions  to  perform  a  sequence  of  operations;  or  he 
may  access  the  public  libraries  which  consist  of  commonly  used 
preprogrammed  functions  that  perform  a  variety  of  useful  compu- 
tations . 

In  the  next  section  we  describe  briefly  the  use  of  APL 
in  the  calculator  mode  to  generate,  manipulate  and  store  data. 
In  later  sections  we  describe  the  use  of  APL  to  transform  the 
data,  and  the  use  of  APL  Public  Library  Workspaces.   The  user 
is  directed  to  references  [1,2,3,4]  for  information  about  writing 
functions  in  APL,  and  for  more  detailed  information  about  primitive 
operations  in  APL. 


* 

If  any  symbol  other  than  the  right  parenthesis  )  appears  before 

the  word  LOAD,  the  system  is  not  properly  translating  the 
character  set  into  the  required  APL  characters.   If  this  happens, 
the  user  should  type  the  symbol  that  appears  (represented  here 
by   ^)   followed  by  the  word   OFF   as  follows: 


t^OFF   (CR 

This  will  get  the  user  off  the  system  at  which  time  he  should 
consult  with  the  computer  center  staff  as  to  the  special 
requirements  for  his  terminal.   Do  not  proceed  in  APL  until 
the  character  ")"  appears  before  the  word  LOAD  in  the  message 
typed  when  APL  is  entered. 
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B.   Generating  and  Storing  Data 

Generation  of  data  arrays  is  very  simple  in  APL  because 
one  need  not  be  concerned  with  format.   One  must,  however,  become 
familiar  with  the  APL  syntax.   First  is  the  assignment  symbol 
■*-     which  plays  a  role  somewhat  like  the   =   sign   in  FORTRAN. 
It  means  to  take  the  expression  on  the  right  and  assign  it  to 

the  variable  named  on  the  left.   Suppose,  for  example,  one  wants 

* 
to  generate  a  vector  of  data  named   X   consisting  of  the  four 

observations:   3.15   12.5  7   8   6.00  3.   One  need  only  type: 

X  <-    3.15   "12.57    8    6.003   (cr) 

with  the  decimal  typed  if  needed  and  with  one  or  more  spaces 
serving  as  delimiters  for  the  separate  observed  values.   (Real 

numbers  may  be  expressed  in  scientific  format  by  use  of  the   E 

5 
notation.  2.5E5  means  2.5  x  10   or  250000.)   Extra  spaces  before 

and  after  the  assignment  arrow,  ■*- r    are  not  needed,  but  they  will 

not  hurt  anything.   Additional  values  can  be  added  to  the  vector 

by  using  the  catenate  operator  (the  comma)   as  follows: 


Just  about  any  name  can  be  used  for  variables  as  long  as  the 
first  character  is  alphabetic  and  the  other  characters  are 
alphabetic  or  numeric.   The  character   A   is  also  acceptable. 

One  must  take  care  to  strike  the  characters  for  the  APL 
character  set.   Several  of  the  APL  characters  are  identical 
in  appearance  on  the  terminal  keyboards  to  non-APL  characters, 
but  different  results  are  obtained.   Among  the  characters 
that  appear  the  same  are   )'  (  -  +  x  -j-  /  :  ;  *. 
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X   «-   X,    5.76         4      1.47     (CR 

This    statement   says    to   catenate    the    three    indicated   values    to 
the   end  of    the   old   data   vector      X      and    to    call    the    resulting 
vector      X.      The   entire    vector   can   be    viewed  by    typing      X 
followed  by   the     (cr) 


X     (CI 

3.15       "12.57      8      6.003      5.76       "4      1.47 

In  the  above  vector,  the  negative  signs  are  typed  using  the 
negative  symbol  (upper  shift  2) ,  not  the  minus  sign  (upper 
shift   +) . 

If  one    is    entering  many   data   values    into    a    vector  he 
may    require    several    lines   of   input    to    do   so.      Additional    values 
could  be    inserted,    as    above,    using    the    catenate   operator.      How- 
ever,   a   simpler  way    is    to    type    a   comma    followed  by    the   quad 
symbol      □       (upper   shift   L)    and  a   carriage    return    at   the   end  of 
a   line    if   more    values    are    to   be    continued   on    succeeding   lines. 
The   next   line   will   automatically   begin  with    the   quad   symbol 
prompting   the    user   to   continue   entering   values.      This    is 
illustrated  below: 


DATA   «-    4       6       8      12      6         2         4001       22       36       29       1 
D:    4      2      12      9       "3      "17      9         9      2       "1      4 


,  D    (CR 
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If  a  value  is  typed  incorrectly  and  discovered  before 
the  (CR)  is  struck,  the  value  can  be  corrected  by  backspacing 


* 
to  the  incorrect  entry  and  hitting  the  line  feed  key.    This 

will  erase  the  incorrect  value  and  everything  to  the  right  of 

the  value.   Then  type  the  correct  values. 

Data  correction  or  modification  after  a  line  has  been 

terminated  by  a  carriage  return  can  be  accomplished  several  ways 

One  way  is  to  determine  the  index  of  an  incorrect  entry  and 

assign  a  new  value  to  that  specific  element  of  the  data  array. 

For  example,  suppose  one  wants  to  change  the  third  element  from 

8  to  4.81  leaving  all  the  other  elements  alone.   This  can  be 

done  by  typing: 


X[3]  <-  4.81  (CR) 

Multiple  corrections  can  be  made  simultaneously  as  follows: 

X[3,  4,  5]  <-  4.81   5.9   "3.26   (S) 

The   X  vector  would  then  contain: 

3.15   "12.57   4.81   5.9   "3.26   1.47. 

The  index  operator   \   (upper  shift  I)  is  convenient  for 
determining  the  index  of  a  given  value  (or  indices  of  a  set  of 
values)  in  an  array  of  data.   If   X   is  the  data  array  and  B 
is  a  set  of  values  in   X,  then   X  \  B  will  generate  the  set  of 


* 

On  some  terminals  the  line  feed  key  is  labelled  LF. 
On  others,  the  attention  or  break  key  must  be  struck. 
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indices  of  the  values  of   B   in   X.   If  a  value  appears  in  more 

than  one  place  in   X,  only  the  subscript  of  its  first  appearance 

is  given.   For  example,  if   X  +    3.5   "12.57   4.81   5.9   "3.26   "4   1.47 

the  index  of  the  value  5.9  is  found  as  follows: 

X  i  5.9 
4 

Let   B  «-  "12.5  7   "3.26.   The  indices  of   B   in   X   are 
given  by : 

X  i  B 
2   5 

This  operation  could  be  useful  for  altering  selected  values  in 
an  array.   To  illustrate,  assume  that  we  want  to  change  the 
values   "12.57   and   "3.26   in   X   to  80  and  83.   The  operation 
below  will  accomplish  this  change: 

X[X  i  "12.57   "3.26]  +  80   83 
X 
3.15   80   4.81   5.9   83   "4   1.47 

Additional  values  can  be  inserted  at  any  position  in  the 
vector,  and  values  can  be  deleted  from  the  vector  by  using  the 
take  "t" ,  drop  "+",  and  catenate  ","  operations.   The  take 
operation  with  syntax 

r  t  X, 
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selects  the  first   r   elements  from  the  vector   X   if   r   is 
positive  and  the  last    |r|   elements  if   r   is  negative.   The 
drop  operation  with  syntax 

r  +  X 

deletes  the  first   r  elements  of   X   if   r   is  positve  and  the 
last   |r|   elements  if   r   is  negative.   If   |r|   is  greater  than 
the  number  of  elements  in   X,  r  +  X  will  insert  zeros  to  the 
right  of  the  elements  of   X   until   |r|   elements  are  obtained 
and   r  4-  X  will  result  in  the  empty  vector.   The  examples  below 
demonstrate  the  take  and  drop  operations. 

X 
3.15   "12.57   4.81   5.9   "3.26   "4   1.47 

3  t  X 
3.15   "12.57   4.81 

"2  +  X 
"4   1.4  7 

9  t  X 
3.15   "12.57   4.81   5.9   "3.26   "4   1.47   3.15   "12.57 

2  I  X 
4.81   5.9   "3.26   "4   1.47 

"4  +  X 
3.15   "12.57   4.81 
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The  take  and  drop  operators  can  be  used  in  conjunction 
with  the  catenate  operator  to  edit  data  vectors.  The  examples 
below  illustrate  some  ways  this  can  be  done. 

Let 

X+-12   34   5   6   78   9   10. 

Create  a  new  vector   Y   consisting  of  the  first  four  elements 
of   X   followed  by  the  elements   11   12   and  then  the  last 
three  elements  of   X;  i.e.   Y  =  1   2   3   4   11   12   8   9   10. 


Y  +    (4  t  X)  ,  11   12,  ("3  t  X)   (CR 
or       Y  <-    (4  f  X)  ,  11   12,  (7  4-  X)    (c 

Let  W  «-  0   ~1   ~2   ~3   ~4   ~5.   Create  a  new  vector   Y   con- 
sisting of  the  first  eight  elements  of   X,  followed  by  the  number 
25,  followed  by  the  middle  four  elements  of   W. 


Y  <-   (8  t  X)  ,  25,  (1  I  (5  t  W)  )    (CR 

A  few  comments  about  the  above  operation  are  in  order.   APL 
always  operates  from  right  to  left  except  when  parentheses  are 
used  to  override  the  standard  order  of  operation.   Thus,  (5  t  W) 
is  the  first  operation  executed.   This  results  in  the  vector 

0  "1   ~2   ~3   "4.   The  next  operation  is   (1+0   "1  ~2      "3   "4) 
which  gives   ~1   ~2   ~3   ~4 .   This  is  catenated  to  25  yielding 

25    1   ~2   ~3   ~4    which  is  then  catenated  to   (8  t  X)   giving 

1  2   3   4   5   6   7   8   25   "1   ~2   "3   "4. 
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In  data  analysis  one  often  wants  to  compare  different 
groups  of  data  or  to  relate  one  group  to  another  group.   Many 
of  the  functions  in  OA  36  60  require  the  data  arrays  be  matrices 
whose  columns  represent  the  different  groups  and  whose  rows 
represent  the  various  observations  for  each  group.   Let  us  see 
how  a  data  matrix  can  be  created.   Suppose  for  example  that  one 
wants  to  create  the  matrix: 


A  = 


1 
3 
5 

7 


The  easiest  way  to  create   A   is  to  create  a  vector  consisting 
of  the  eight  elements  and  then  reshape  the  vector  into  the 
desired  4x2  array.   The  reshape   operator   p  can  be  used  to 
reshape  a  vector  into  any  specified  size.    The  matrix  A   is 
created  by  the  operation: 


A  -*■  (4,2) 


12345678 


The  general  syntax  for  the  reshape  operation  is   A  -*-  (N,M)pB 
where   B   is  a  vector,  or  scalar,  M   is  the  number  of  columns  of 
the  array   A,  and  N   is  the  number  of  rows  of  A.   The  first   M 
elements  of   B  will  be  used  for  row  1,  the  second  M   for  row  2, 
etc.   If  there  are  less  than   M*N   elements  of   B,  the  elements 
will  be  repeated  from  the  beginning  as  many  as  are  needed.   If 
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B   contains  more  than   M-N   elements,  only  the  first   M»N   will 
be  used.   Try  a  few  examples  to  see  what  reshape  does. 

Remember  that  the  reshape  operator  fills  up  matrices 
row  by  row.   Suppose  that  the  user  entered  the  data  into  the 
vector  by  columns;  e.g.  B  -^  1   3   5   7   2   4   6   8.   How  can  the 
4x2  matrix  A  be  generated?    (4,2) pB  will  not  work  since 
that  would  give 


However,  if  one  requested   (2,4)pB   one  would  obtain 

13   5   7 
2   4   6   8 

What  we  need  now  is  to  transpose  this  to  obtain   A.   The  transpose 

* 
operator  is   §   (upper  shift  0  overstruck  with  upper  shift  /  ) . 

This  sequence  of  operations  can  be  combined  as  follows : 


A  <-  $(2,4)  pB      (CR 

The  operator   p  is  also  useful  for  determining  the 
dimension  or  shape  of  an  arbitrary  variable.   When  used  for 
this  purpose  it  is  called  the  shape  operator.   The  syntax  is 

PZ   (cr" 


* 
Type  shift  0  (oh) ,  backspace,  and  type  upper  shift   / 
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Table  2.1  illustrates  the  results  of  the  shape  operator  when 
applied  to  selected  variables.   The  first  entry  is  the  scalar  3. 
APL  considers  all  scalars  to  have  no  dimension  (not  dimension  0) 
Thus,   pZ   is  the  empty  vector.   In  the  next  to  last  case,  the 
variable   Z   is  defined  to  be    i0  (iota  0) .   (This  represents 
a  second  use  of  the  index  operator   i   where  the  operator  has 
only  a  right  argument.   If   K   is  any  nonnegative  integer,   iK 
is  the  vector  consisting  of  the  first   K   positive  integers. 
For  example,   i3  would  be  the  vector  1   2   3.)   The  use  of 
the  index  operator  in  Table  2.1  is  a  rather  special  one.   APL 
interprets   i0   as  a  vector  with   0   elements,  hence  the  result 
0  for   pZ.   This  probably  is  puzzling  to  the  reader,  but  it  is 
really  quite  handy.   The  reader  should  simply  commit  to  memory 
the  facts  that  scalars  have  no  dimension  and  the  empty  vector 
has  dimension  0.   In  the  last  entry  of  Table  2.1,  Z   is  defined 
as  the  1x3  array  consisting  of   1   2   3.   This  array  appears 
to  be   identical  to  the  vector   Y  «-  1   2   3,  but  APL  makes  a 
distinction  between  the  two. 


VARIABLE 

pZ 

Z    «-    3 

- 

Z    *-    2       3       5 

3 

Z    +    (2,2)p       1       2 

3 

4 

2      2 

Z    «-    iO 

0 

Z    +    (l,3)p               2 

3 

5 

1       3 

Table  2.1:   THE  SHAPE  OPERATOR 
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The  last  APL  operators  discussed  in  this  section  are 

the  grade  up  and  grade  down  functions.   These  two  functions, 

* 
denoted  h      and   ?,   give  the  indices  of  the  elements  in  a 

vector  in  order  of  magnitude:   ascending  order  for  grade  up 

and  descending  order  for  grade  down.   This  is  very  useful  for 

data  analysis  since  much  of  our  work  requires  us  to  sort  the 

data.   These  operations  accomplish  this  easily.   Here  are  some 

examples : 

A   <-    9       12      6       4       ~2       3       7       "5       11 
£A 

856437192 

291734658 

In  the  grade  up  of   A   the  first  element,  8,  says  that  the 
eighth  value  of   A   should  be  taken  first;  the  second  element, 
5,  says  that  the  fifth  element  of   A   should  be  taken  next, 
etc.,  to  sort  A   in  increasing  order.   Similarly  for  grade  down 
to  sort  in  decreasing  order.   The  sorted  arrays  can  be  written 
as  follows: 

AUA] 

~5       "2      3      4      6       7      9      11      12 
A[tA] 

12      11      9       7      6      4       3      ~2      ~5 


* 
These  symbols  are  upper  shift  H   and   G,  respectively,  over- 
struck  with  upper  shift   M. 
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We  have  only  scratched  the  surface  of  APL  features 
useful  for  generating,  modifying  and  storing  data.   The  collection 
of  operations  discussed  here  will,  nevertheless,  enable  a  user 
to  do  most  of  the  things  that  he  requires  to  generate,  correct, 
modify  and  store  data.   In  the  next  section  we  discuss  arithmetic 
operations  in  APL  so  that  the  user  can  perform  calculations  and 
transform  data  arrays. 


C.   Arithmetic  Operations 

In  this  section  we  discuss  some  of  the  basic  arithmetic 
operations.   We  restrict  attention  to  those  operations  most  fre- 
quently used  in  data  analysis. 

The  symbols  used  for  the  basic  operations  of  adding,  sub- 
tracting, multiplying,  and  dividing  are  the  standard  ones  used 

* 
in  mathematics:   (+,  -,  x,  v).    However,  in  APL  these  operations 

are  more  powerful  than  their  equivalents  in  most  programming 

languages  since  the  operations  can  be  applied  to  entire  vectors 

or  matrices.   When  used  with  vectors  or  arrays,  the  operations 

are  applied  componentwise.   This  requires  that  the  vectors  or  arrays 

have  the  same  size.   An  exception  is  made  for  the  case  where  one 

of  the  arguments  is  a  scalar.   When  this  happens  APL  adds, 


* 
The  subtraction  sign  is  upper  shift   +.   Do  not  confuse  this 
with  the  negative  sign  which  is  upper  shift  2. 
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subtracts,  multiplies  or  divides  the  scalar  to  each  element  of 
the  vector  or  array.   Of  course,  the  operation  can  be  applied 
to  two  or  more  scalars  as  with  any  other  programming  language . 
The  examples  below  illustrate  the  basic  operations: 


A  <-  3   5   7   8     (CRJ 

B  «-  ~1   "2   1   2   (CJ 

A  +  B  (Q) 
2   3   8   10 

A  -  B  (S) 
4   7   6   6 

A  x  B  (CR) 
~3   "10   7   16 

A  v  B  (S) 
~3   "2.5   7   4 

2  x  a 

6  10   14   16 

3  +  B 

2  14   5 

C  +    (2,2)p 

3  5 

7  8 
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D    ^    (2,2)pB        (C 
D       (CR 


C    +   B 


1  2 

1  2 

2  3 
8      10 

C 
1.5      2.5 
3.5      4.0 

27 

2 

The    last   two  examples   are    included   to   re-emphasize    the 
right-to-left  order   of   operations    in   APL .      Thus,    5    is    added   to 
4    to    give    9   which   is    then   multiplied  by    3   to   yield   27.      Similarly, 
in    the    last   case,    2    is   multiplied  by    3    to   give    6   which    is    then 
divided   into    12    to   yield   2.      Unlike    some   other   languages    there 
is   no  hierarchy   of   operations   other    than    right    to    left.      Of 
course,    parentheses   may  be    used  as    in   algebra    to   change    the 
priority   of  operation. 

(3x4)    +5       (cr) 
17 

(12    v    3)     x    2    (cr) 
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In  addition  to  the  four  basic  operations,  data  analysis 
frequently  requires  power,  log,  reciprocal,  and  exponential 
transformations.   As  with  the  four  basic  operations  these  can 
be  applied  to  scalars,  vectors,  or  matrices. 

Table  2.2  presents  the  syntax  for  these  operations. 


Trans  formation 

Syntax 

xa 

Jin  X 

1/X 

X 
e 

X  *  a 
©X  l 

tX 

*X 

Table  2.2.   SYNTAX  FOR  DATA  T RAN SFORMAT IONS 


Here  we  give  examples  of  these  transformations 


X+-12345   (CR) 

X  *  0.5   (cr) 
1   1.41421   1.73205   2   2.23607 

X  *  2    (6r) 
1   4   9   16   25 

X    *    ~3        (CJ 
1      0.125      0.037037 


The  symbol  for  the  logarithm  is  the  upper  shift  0  (oh)  over- 
struck  with  the  asterisk  (upper  shift  P) . 
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OX        (CR) 

0  0.693147      1.098612      1.386294      1.609438 

fX        @ 

1  0.5      0.333333      0.25      0.2 

X    *    t    "3        (£5) 

1  0.793701      0.693361      0.629961      0.584804 

(Recall  that  APL  operates  from  right  to  left  so  that  the  operator 
above  raises  each  element  in   X   to  the   -1/3   power.) 

*X   @ 
2.71828   7.389056   20.085537   54.598150   148.413159 
The  logarithm  of  a  number  N   to  an  arbitrary  base   B   can  be 
determined  by  typing   B  ®  N ;  e.g., 

10  ®  1000 
3 

Three  other  useful  arithmetic  operations  for  scalars, 
vectors,  and  matrices  are  the  ceiling  [,    floor   L,  and 
absolute  value   | .   All  of  these  can  be  used  with  a  single  right 
hand  argument  or  with  both  left  and  right  arguments. 

The  examples  below  illustrate  the  use  of  these  operators: 

Monadic  (one  argument) 
T2.75   6   0.08   "3.6 
3   6   1   ~3  (Gives  smallest  integer  >^  argument.) 

L2.75   6   0.08   "3.6 

2  6      0      ~4  (Gives    largest  integer    <_  argument.) 

|2.75      6      0.08      "3.6 
2.75      6      0.0  8      3.6         (Gives    absolute   value    of   each   element. 
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Dyadic  (two  arguments) 

3   5   T   2  8 

3   8  (Gives  the  maximum  of  each  component.) 

3   5   L   2  8 

2   5  (Gives  the  minimum  of  each  component.) 

3   |   7   6   8   "2 

* 
10   2   1    (Gives  the  remainder  after  dividing  each  element  by  3.) 

Clearly,  the  floor  and  ceiling  operators  are  useful  for 
rounding  values.   One  can  select  any  number  of  significant  digits. 
If  one  wants  to  round  values  to  the  nearest  integer,  one  should 
type   L0.5  +  N.   If,  say,  four  significant  digits  are  wanted, 
one  could  type  (L.5  +  N  x  1E4)  f  1E4 .   The  examples  below 
illustrate  the  use  of   L   for  rounding: 

N  «■  0.0835126   12.51877623    1.33333333   5.25 

L  .5  +  N 
0   13   1   5 

(L.5  +  N  x  1E4)  t  1E4 
0.0835   12.5188   1.3333   5.2500 

/Another  operator,  the  reduction  operator   /,  allows 
the  operations  discussed  above  to  be  applied  to  all  the  elements 
of  a  vector  or  to  the  rows  or  columns  of  a  matrix.   The  syntax  is 


*         .  i  .... 

All  remainders  of   K|N  are  expressed  as  positive  integers  in 

the  set   0,  iK-1.   If  a  remainder  is  negative,  K   is  added  to 

it. 


25 


f/A     where      f      can   be    any   of   the    arithmetic   operations    discussed 
above   and     A      is    a   vector   or   a  matrix.      The   result  of   this   oper- 
ation when   applied   to   a   vector      A     of   size      n      is    the    scalar 
A[l]    fA[2]f     •••       fA[n].      When    applied   to    a   matrix      M      having      r 
rows    and      c      columns,    the   result   is    the   vector      B      of   dimension    r 
where      B[i]    =   A[i;l]    fA[i;2]    f    •••    fA[i;c]       for      i    =    l,2,...,r 
(column    reduction) .      The   matrix      M      can   be    reduced  over   its 
rows   by   typing      f/[l]M.      The  examples   below   illustrate    the    reductio 
operation : 

+   /l      3      7      4 
15  (the   sum  of   all    the   elements) 

x/1      3      7      4 
84  (the   product  of   all    the   elements) 

r/1      3      7      4 
7  (the    largest   element) 

M  «.    (3,2)p3      12      6      5      4 

M 
3      1 
2      6 
5      4 

-/M 
2    4   1      (differences  of  elements  in  rows  1,  2,  and  3) 

L/M 
12   4       (minimum  values  in  rows  1,  2,  and  3) 
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+/[1]M 
10   1     (sums  of  elements  in  columns  1  and  2) 

T/[1]M 
5   6      (maximum  values  in  columns  1  and  2) 

Finally,  we  describe  the  inner  product  operation  and 
the  matrix  inverse.   The  inner  product,  like  the  reduction 
operation  can  be  applied  with  any  general  APL  operators .   The 
syntax  is   Af.gB,  where   A   and   B   are  vectors  or  matrices 
(which  must  satisfy  certain  size  restrictions)  and   f   and   g 
are  any  general  APL  operators.   For  vector  arguments,  A   and   B 
must  be  the  same  size,  say   n,  and  the  result  is  given  by 
f/A[l]gB[l]   A[2]gB[2]    ...    A[n]gB[n].   For  matrix  arguments 
the  number  of  columns  of   A  must  be  the  same  as  the  number  of 
rows  of   B.   The  result  is  a  matrix  of  size  n  x  m  where 
pA  =  (n,k)   and  pB  =  (k,m).   The  (i,j)th  element  of  the  result 
is   f/A[i;l]gB[l; j]   A [i ; 2 ] gB [2 ; j ]   ...   A [i ;k] gB [k ; j ] .   The 
reader  should  recognize  the  operation  above  as  matrix  multi- 
plication when   f   is   +   and   g   is   x.   The  examples  below 
illustrate  the  inner  product. 

X  <-  1   4   3   2 

Y  <-  2   3   1   2 

X+  .  x  y 
21   ((1x2)  +  (4x3)  +  (3x1)  +  (2x2)) 

X  x  .  +  Y 
336  ((1  +  2)  x  (4  +  3)  x  (3+1)  x  (2  +  2)) 
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A  +■  3 
A 


1   2 

3   4 


B<-2   3p35   2   311 

B 
3   5   2 
3   11 

A  +  .  x  B 
7    7    4     (the  matrix  product  of  A   and  B) 
21   19   10 
33   31   16 

The  domino  operator   !   (type   D,  backspace,  *)  is  used  to 
solve  for  the  matrix  inverse  of  a  nonsingular  square  matrix, 
The  inverse  of  a  nonsingular  square  matrix  C   is  found  by 
typing   !C.   For  example, 

110   2 


C   «-    2      2      p 

C 

1 

1 

0 

2 

@C 

1 

-0.5 

0 

0.5 

C    +    .     x    @c 

1 

0 

0 

1 
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D.   Workspace  Management 

In  the  previous  sections  we  have  seen  how  data  are  input 
and  variables  are  transformed.   In  this  section  we  describe  how 
the  user  manages  his  APL  workspace  so  that  he  can  save  data  for 
use  from  one  session  to  the  next;  he  can  load  public  library 
workspaces;  and  he  can  use  functions  available  in  the  public 
library  workspaces.   For  information  about  writing  functions  the 
user  should  see  references  [2,3,4]. 

When  a  user  types  APL  he  is  put  into  a  clear  APL  work- 
space.  In  this  clear  workspace  data  can  be  created  and  any  of 
the  operations  described  in  the  earlier  chapters  can  be  performed. 
If  the  user  wishes  to  maintain  the  data  or  results  for  future 
use  he  must  save  the  workspace.   This  is  done  by  typing 

) SAVE  WSNAME 

where  WSNAME  is  an  arbitrary  name  (first  letter  alphabetic, 
eleven  characters  or  less)  that  the  user  selects  for  the  workspace 
This  private  workspace  can  subsequently  be  loaded  into  the  user's 
active  workspace  by  typing 

)LOAD  WSNAME 

This  causes  the  active  workspace  to  be  cleared  and  a  copy  of 
the  named  workspace  to  be  written  into  the  active  workspace. 
Alternatively,  if  the  user  does  not  want  to  clear  out  the  contents 
of  the  active  workspace  to  bring  in  the  named  workspace,  he  can 
type 
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) COPY  WSNAME 

This  will  simply  augment  the  existing  contents  of  the  active 
workspace  with  the  contents  of  workspace  WSNAME.   However,  if 
an  existing  variable  or  function  in  the  active  workspace  has  the 
same  name  as  a  function  or  variable  in  workspace  WSNAME  the 
latter  will  replace  the  former. 

Any  modifications,  additions,  deletions,  corrections, 
etc.  that  the  user  makes  to  a  copy  of  a  workspace  will  affect 
the  active  workspace,  but  will  not  affect  the  permanent  copy  of 
the  workspace  maintained  on  his  private  files  unless  the  user 
saves  the  so  modified  active  workspace.   This  can  be  accomplished 
by  typing    ) SAVE  WSNAME   (or  simply    ) SAVE   if  workspace  WSNAME 
were  LOADed) . 

The  user  can  determine  the  name  of  the  active  workspace 
at  any  time  by  typing    )WSID   (for  Workspace  IDentification) . 
Similarly,  he  can  change  the  name  of  the  active  workspace  by 
typing   )WSID  NEWNAME .   A  list  of  all  the  user's  private  APL 
workspaces  is  obtained  by  typing   )LIB.   Entire  workspaces  can 
be  permanently  destroyed  by  typing   ) DROP  WSNAME.   (Be  careful 
with  this  one  I) 

Within  an  active  workspace,  a  user  can  obtain  a  list  of 
all  of  the  functions  contained  in  the  workspace  by  typing  )FNS. 
The  functions  will  be  listed  alphabetically  so  that  all  functions 
from  those  beginning  with  a  certain  letter,  LETTER,  onward  can 
be  obtained  by  typing   ) FNS  LETTER.   Similarly,  a  list  of 
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variables  in  the  workspace  can  be  obtained  by  typing   ) VARS   or 
) VARS  LETTER.   Typeout  can  be  terminated  at  any  point  by 
hitting  the  BREAK  or  the  ATTN  key.   Variables  and/or  functions 
can  be  deleted  from  the  active  workspace  by  typing   ) ERASE  LIST 
where  LIST  is  a  single  function  or  variable  or  a  list  of  functions 
and  variables  to  be  deleted.   Names  in  the  list  should  be 
separated  by  one  or  more  blanks.   Two  other  system  commands 
allow  the  user  some  control  over  his  APL  environment.   These 
commands  control  the  number  of  digits  printed  out  and  the  width 
of  a  typed  line.   They  are   ) DIGITS  N   and   ) WIDTH  N,  where   N 
is  the  number  of  digits  to  be  printed  or  the  desired  line  width, 
respectively. 

In  addition  to  his  own  private  APL  workspaces,  the  user 
also  has  access  to  all  of  the  APL  public  libraries  available 
at  NPS.   These  libraries  have  numbers  between  1  and  999  and  are 
intended  to  hold  workspaces  of  general  interest.   See  [1]  for 
a  list  of  the  public  library  numbers.   The  contents  of  a  public 
library  can  be  displayed  using  the   )LIB   command  followed  by 
the  library  number.   For  example, 
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)LIB  1 


A-DISK 

R/0 

IOFNS 

17.25 

7/09 

PLOTFORM 

17.25 

7/09 

NEWS 

15.54 

8/2  7 

WSFNS 

17.25 

7/09 

TEXTEDIT 

17.26 

7/09 

FORMAT 

17.26 

7/09 

CATALOG 

15.17 

7/14 

MAILBOX 

12.50 

7/14 

MULTIPLO 

14.06 

7/18 

FILEFNS 

17.18 

8/14 

The  list  contains  the  workspace  name  and  the  time  and  date  that 
it  was  last  modified. 

A  copy  of  a  public  library  workspace  can  be  put  into  the 
user's  active  workspace  by  typing   ) LOAD  n  WSNAME   or 
) COPY  n  WSNAME   where   n   is  the  library  number  and  WSNAME  the 
workspace  name.   (Recall  the  differences  in  LOAD  and  COPY;  LOAD 
will  first  clear  the  contents  of  the  active  workspace.)   Selected 
functions  and  variables  from  a  workspace  (either  public  or 
private)  can  be  copied  by  typing 


) COPY  n  WSNAME  OBJECT 
or 

) COPY  n  WSNAME  GROUP 
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where  OBJECT  is  a  single  variable  or  function  that  the  user  wants 
to  copy  into  his  active  workspace  and  GROUP  is  a  group  of  variables 
and/or  functions  that  has  previously  been  defined  for  WSNAME.* 

When  a  public  library  workspace  is  brought  into  the 
user's  active  workspace,  documentation  can  usually  be  obtained 
by  typing  the  word  DESCRIBE.   It  is  somewhat  standard  procedure 
to  document  public  library  workspaces  with  a  DESCRIBE  variable 
giving  general  information  and  with  each  function  in  the  work- 
space being  documented  by  a  "HOW"  variable.   Type  the  name  of 
the  function  followed  by  HOW  with  no  intervening  spaces.   For 
example,  documentation  for  the  function  STEMLEAF  in  public 
library  2  OA3660  is  obtained  by  typing  STEMLEAFHOW.   "HOW" 
variables  generally  describe  the  function  syntax,  parameters, 
input/output  requirements,  etc.   Of  course,  the  user  can  determine 
if  a  workspace  contains  this  sort  of  internal  documentation 
by  typing   ) VARS . 

When  the  user  has  completed  his  work  session  he  can  logout 
by  typing   )OFF.   If  he  wants  to  get  out  of  APL  but  not  logout, 
he  should  type  ) OFF  CMS.   Then,  he  should  switch  back  to  the 
standard  keyboard. 


* 
The  group  structure  for  the  OA3660  workspace  is  discussed  in 
Chapter  III. 
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E.   Error  Recovery 

This  section  describes  recovery  procedures  that  the  user 
can  employ  if  he  receives  error  messages  during  execution 
of  functions  contained  in  a  workspace.   Usually,  error  messages 
received  when  functions  in  public  library  workspaces  are  executed 
are  the  result  of  improper  syntax  or  problems  with  the  shapes 
of  the  arguments.   As  soon  as  an  error  is  detected,  execution 
of  the  function  is  suspended,  the  number  of  the  line  containing 
the  error  is  typed,  a  caret  is  inserted  at  the  position  in  the 
line  that  the  error  occurred,  and  an  explanatory  error  message 
is  printed  out.   When  the  function  is  suspended  the  values  of 
all  variables  determined  up  to  the  point  of  suspension  can  be 
obtained  by  simply  typing  the  variable  name.   All  of  this  makes 
error  discovery  and  correction  quite  simple.   From  the  point  of 
suspension  the  user  can  branch  to  any  line  of  the  function  by 
typing  ■*■   n,  where   n   is  the  desired  line  number.   If  n   is 
omitted  a  branch  is  made  outside  of  the  function.   Many  types  of 
calculations  can  be  performed  while  a  function  is  suspended, 
including  the  execution  of  other  functions.   However,  the  function 
will  remain  suspended  until  a  branch  is  made  as  described  above. 
Suspended  functions  tend  to  clutter  up  the  user's  workspace  so 
that  suspended  functions  should  not  be  left  pending.   A  list 
of  all  suspended  functions  can  be  displayed  at  any  time  by 
typing   )SI   (state  indicator).   If  several  functions  are  listed 
as  being  suspended,  the  user  should  type  as  many  branch  arrows 
(+) ,  one  per  line,  as  there  are  asterisks  displayed  in  the  list. 
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Since  the  functions  contained  in  the  public  library 
workspaces  have  been  tested  extensively,  most  of  the  errors 
that  are  encountered  result  from  improper  function  syntax  or  the 
use  of  arguments  that  are  not  conformable  or  of  improper  size. 
The  documentation  contained  in  the  next  chapter  shows  the  proper 
syntax  for  each  function,  and  gives  the  requirements  placed  on 
the  arguments. 
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III.   THE  OA3660  WORKSPACE 

This  chapter  provides  documentation  for  the  public 
library  workspace,  2  OA3660.   General  workspace  documentation 
is  contained  in  the  DESCRIBE  variable  and  in  the  lists  of 
functions,  variables,  and  groups.   Short  writeups  are  given 
for  each  function  in  the  workspace  via  "HOW"  variables,  and 
examples  of  the  use  of  each  function  are  provided.   Data  are 
contained  in  the  workspace  to  provide  easy  illustration  of  the 
functions.   The  data  are  described  in  the  DATAMAP  variable. 
Finally,  the  actual  APL  programs  are  displayed. 

A.   General  Documentation 

The  OA3660  APL  workspace  is  contained  in  public  library  2. 
Therefore,  the  user  must  type   ) LOAD  2  OA3660   or   ) COPY  2  OA3660 
to  create  a  copy  of  the  OA36  60  workspace  in  his  active  workspace 
area.   The  functions  and  variables  contained  in  OA3660  are 
displayed  by  typing   ) FNS   and   ) VARS   as  shown  below. 


)FNS 

AND  ALU  OVA         A3ti  A3RSR         BOXPLOT    CHI  SQUARE  CODERVS 

COMPARE    CONDENSE  CONTINGENCY  CORRELATION  FILL  P.J 

INPUT         US  LINE  LIT  LSLINE      MEDPOLISti  rtSTATS      NUi 

ONEti  0NE3  PARTIAL    REGRESS    SCAT  SHOtfRES    SPLIT  STATISTICS 

STEMLEAF  SUMSQ         TWICE         UTCOND 
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)VARS 

AGES  ANDHO'W       ANOVAHOW  AZRHO'W       A3RSRH0W  BOXPLO;   .     . 

CtilCHHTS  CHISQUAREHO'W  COAL  CODERESilOW  COMPAl      . 

C  0  N  U  E  N  St;  HOW  CO  N  S  U  MP  TION  CON  T I NG  ENCY  H  0  W      C  0  R  h  fi  LAI  J  OH  ii  0  W 

CRIMES       DATAmAP    DEATHS       DEP  DEPTH  DESCRIBE  DISCO. 

DRAPE*       EPSILON    GhADESl    GRADES2  H  Y DROP LANTS  ID/SECTS    KSHOW 

LASD AREAS  LENGTH       LIMEHOW  LSLT.NEtiOW  ''-iEDPOLI 'SiitiOW  -iISS 

MORTALITY  MSTATSHOW  NDIVX         NOISY  VGAP  "To 

NUM  NUMSU14H0H  ONEUUOW  0NE3H0W    OPTION      PARTI ALHOW 

PRECIPITATION         PRES RATING  REGRESS  HOW  REST  OS       RIVERS 

RSELECT    SCALE         SCATHOW    SHOWRESHOW  SPLITHO'W  STATISTIC: 

STEmLSAFHOW  TESTX         TESTY  TWICEHOW  USPOP         VOLCA  ■ 

WARPBREAKS  WID  WIDTH  MNTERCEPT  CH 


Notice  that  there  is  a  variable  called:  DESCRIBE.   This  variable 
gives  general  documentation  about  the  workspace. 


fiTRE  DESCRIBE  VARIABLE  PROVIDES  GENERAL  DOCUMENTATION 
f\RGR  THE  ENTIRE  workspace.  SIMPLY  TYPE  'DESCRIBE'  TO 
^OBTAIN    THE    INFORMATION . 


U ESCRIBE 


0A36SQ-*DATA    ANALYSIS 


DATE:       JANUARY     1978 

PROGRAMMER'.       P.     hUSSELL    RTCtiARDS 

THIS    WORKSPACE    CONTAINS    FUNCTIONS    AlvD    VARIABLES       USEFUL       t 
EXPLORATORY    INTERACTIVE    DATA    ANALYSIS.       FOR      INFORMATION       ON       .   . 
USE    OF    THE    FUNCTIONS    TYPE    THE    FUNCTION    NAME    FOLLOWED    BY        'HOW'. 
FOR    EXAMPLE,     TYPE    ST  EM  LEAF  HOW    FOR    DOCUMENTATION    ON       THE       FUNCTION 
STEMLEAF.       SOME    FUNCTIONS    ARE    USED    ONLY    AS    SUBPROGRAMS    FOR      OTh 
FUNCTIONS    AND    WILL    in  Of    HAVE    A     'HOW'     VARIABLE.       FOR     iOhE       DETAT, 
DOCUMENTATION    SEE    THE    NFS    TECHNICAL    REPORT,     A    USER'S    GUIDE    TO 
6M3660  APL      WORKSPACE,        BY       E.       RUSSELL       RICHARDS       Oh       DONALD 
MCNEIL'S    BOOK,    INTERACTIVE    DATA    ANALYSIS. 

THE    WORKSPACE    ALSO    CONTAINS       VARIOUS      SMALL       DATA       SETS       FUh 
ILLUSTRATION       OF       THE       FUNCTIONS.  TYPE       UATAnAP         TO         OBTA 

IDENTIFICATION    OF    THE    DATA    AhRAYS. 

'THE    USER    CAN    SAVE    SPACE    IN    HIS    ACTIVE    WORKSPACE       BY       CO,    tlNG 
ONLY    THOSE      FUNCTIONS       AND      VARIABLES       iEEDED      TO         PERFORM 
iEC&SSARX       TASKS.  FUNCTIONS       AND       VARIABLES       ARE       CONVENIENTLY 

COLLECTED       INTO       GROUPS       WHICH       ARE       IDENTIFIED       IN  T H  ' 

VARIABLES.       FOR    EXAMPLE,     IF    THE    USER    WANTS    TO       GENERATE      STt     . 
OR    BOXPLOT    DISPLAYS    HE    CAN    COPY    THE    DISPLAY    GROUP    AS    FOLLOa     : 

)COPY    2  0/13660  DISPLAY 

THE    GROUP    CONTAINS    NOT       ONLY       THE       FUNCTIONS       BUT       ALSO 
REQUIRED    PARAMETERS    AND    DOCUMENTATION . 
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There  is  also  a  variable  named  DATAMAP  which  describes  the 
data  contained  in  OA3660  and  provides  a  reference  to  the  data  source 
This  is  shown  below. 


fiTRE    VARIABLE    DATAMAP    PROVIDES    INFORMATION    ABOUT    TdK 
fiUATA    SETS    CONTAINED   IN    THE    WORKSPACE,       IT    GIVES    THE 
ftSOVRCE    Of    EACH    DATA    SET    AND    TdE    SIZE    OP    Taa    DATA    ARRAY. 
fiSInPLX    TYPE     'DATAMAP'    FOR    DATA    DOCUMENTATION. 


DATAMAP 

THIS    WORKSPACE    CONTAINS    SEVERAL    DATA    SETS    THAT      C An       bE      USED       TO 
ILLUSTRATE    THE    FUNCTIONS    OR    TO    SERVE    AS    HOMEWORK.    PROBLEMS,  THE 

ENTIRE    COLLECTION    OP    DATA    ARRAYS    IS    CONTAINED    IN    THE    GnOUP      rjAriED 
DATA.       THE    DATA    C AhE    FROM    THE    FOLLOWING    SOURCES: 

(l)DRAPER,     N.h.     AND    H.    SMITH,     APPLIED       REGRESSION       ANALYSIS, 

WILEY    AND    SONS. 

(2)MCNEIL,    D.R.,    INTERACTIVE    DATA    ANALYSIS,    ADDISON^WESLEY . 

{3)MILLER,    I.    AND    I.E.    FREVND,       PROBABILITY      Ain)      STATISTICS 

FOR    ENGINEERS,    PRENTIC E-rt ALL . 

(h)RICHARDS,    F.R.,    CLASSROOM    GRADES. 

(S)TVKEY,    JOHN*    EXPLORATORY    DATA    ANALYSIS,     ADD1 'SON  - JES LEY . 
THE    TABLE    BELOW    LISTS    EACH    DATA    ARRAY,    THE    DIMENSION,    THE      SOURCE 
(USING    AUTHOm    ABBREVIATIONS),     AND    PAGE    OR    CHAPTER    NUMBERS. 
DATA  qDATA  SOURCE 

AGES  4  2  MC(P.H) 

CHICKWTS  14  6  MCiP, 30) 

COAL  49  TV 

CONSUMPTION  5  5  MCiP.100) 

CRIMES  50  4  MC (P. 132) 

DEATHS  5  4  MC {P. 94) 

DISCOVERIES  100  MC (P. 121) 

DRAPERX  13  4  DS { P. l^S) 

DRAPERY  13  DS( P. 17  8) 

GRADE SI  15  5  RI 

GRADES2  15  5  RI 

HYDROPLANTS  3  4  Tu(CH.3) 

INSECTS  12  6  MC(P. 12, 36) 

LAN  DARE  AS  4  8  MC(P.S) 

MORTALITY  5  4  TV (CH. 16) 

PRECIPITATION  6  9  MC(P,3) 

PRESRATTNG  114  MC (P. 126) 

RIVERS  141  MC(P.IH) 

TESTX  10  2  MF (P. 253) 

TESTY  10  MF(P.2S3) 

US  POP  19  2  TU(CH.Q) 

VOLCANO  219  TV 

WARPBREAKS  9  6  MC ( P . 2 8 ) 


Finally,  notice  that  there  are  "HOW"  variables  for  most 
of  the  functions  contained  in  the  workspace.   (The  only  functions  with- 
out  "HOW"  variables  are  utility  functions  which  are  used  by  other 
functions,  but  which  are  transparent  to  the  user.)   These  variables 
provide  documentation  on  the  use  of  the  functions.   For  example, 
documentation  on  the  function  SCAT  is  obtained  by  typing 
SCATHOW. 


SCATHOW 
SYNTAX:       SCAT    ARRAY 
PARAMETERS: 

(1)  WID~    CONTROLS    THE    HORIZONTAL    SIZE       OF       THE       DISPLAY 
(DEFAULT=3Q    CHARACTERS) . 

(2)  DEP*    CONTROLS    THE      VERTICAL      SIZE       OF      TtiE       DISPLAY 
(DEFAULT=15    LINES). 

(3)  NDIVX,    NDIVY-    NUMBER    OF    UNITS      OH      X-      AND      Y-AXES, 
RESPECTIVELY    ( DEFAU LT=^ ,4) . 

GROUP:       RELATIONS ,    SMOOTH,    COMPARISONS 

DESCRIPTION:  SCAT  PRODUCES  A  SCATTER  PLOT  OF  THE  DATA  CONTAINED 
IN  ARRAY.  THE  ARGUMENT  ARRAY  CAN  HE  A  VECTOR  OR  A  MATRIX  JITH  AS 
MANY  AS  9  COLUMNS.  IF  A  VECTOR  OF  SIZE  U ,  THOSE  VALUES  ARr, 
PLOTTED  VS.  THE  INTEGERS  1  TO  N  ;  IF  A  MATRIX,  THE  SECOND,  THIRD, 
ETC.  COLUMNS  ARE  PLOTTED  VS.  COLUMN  1  ON  THE  SAME  AXES. 
DIVISIONS  ON  THE  AXES  OF  THE  PLOT  ARE  NOT  EXPLICITLY  PRINTED, 
EXCEPT  AT  THE  EXTREMES  OF  THE  PLOT.  THE  USER  CAN  CONTROL  THE 
RESOLUTION  OF  THE  PLOT  tiY  MODIFYING  THE  PARAMETERS  WID,  DEP, 
NDIVX    AND    NDIVY.  PRINTING       TIME      INCREASES       DRAMATICALLY       ;U?H 

RESOLUTION;    THEREFORE,    OTHER    PLOT    PROGRAMS    SHOULD    BE    USED    IF    HIGH 
RESOLUTION    IS    DESIRED.    FOR    A    SINGLE    GROUP   OF    DATA,    THE      NU!4BEh 
(2<N<3)    MILL    BE    PRINTED    IF    N    POINTS    LIE      CLOSE      TOGETHER      ON       THE 
DISPLAY.       FOR    MULTIPLE    PLOTS    ON    THE    SAME    DISPLAY,       THE      LETTER      A 
REPRESENTS    GROUP    1,    b    REPRESENTS    GROUP    2,       ETC.  FOR      2  POINTS 

CLOSE    TOGETHER    THE    LETTER    WILL    BE    PRINTED   tflTR       AN      UNDERSCORE. 
THE    DISPLAY    CANNOT    HANDLE    3  OR      MORE      POINTS      CLOSE      TOGETHER      IF 
THERE    ARE    MULTIPLE    PLOTS. 


The  format  of  SCATHOW  is  followed  for  every  function.   The 
function  syntax,  the  user  controlled  parameters,  the  group  (or 
groups)  containing  the  functions,  the  subroutines  used  by  the 
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function,  and  a  brief  description  of  the  use  of  the  function 
are  displayed. 

The  functions  and  variables  are  grouped  into  APL  groups 
for  ease  of  handling  and  to  conserve  space.   The  group  structure 
is  described  below. 


)GRPS 
COMPARISONS  DATA  DISPLAY    DOC  ESSENTIALS  GFIT 

RELATIONS  SMOOTH.       STATS         TwOWAY 


)GRP    COMPARISONS 
CONDENSE  CONDENSEHOW 

SCATHOW    LINE  LINEHOW    NDIVX 

COMPARE    COMPAREHOW  FILL 

ANDHOW       LSLINE      LSLINEHOW 


NUM 

MISS            UTCOHD 

FMT              3C 

NDIVY 

WID              DEP 

DEPTH            iCi 

NUMSUM 

NUMSUMHOW 

RS ELECT    AN 

)GRP    DATA 

AGES  CRICKETS  COAL  CONSUMPTION  CRIMES      DEATHS 

DISCOVERIES  DRAPERX  DRAPERY    GRADES1    GRADES2    HYDROPLA'NTS 

INSECTS    LAN  DARE  AS  MORTALITY  PRECIPITATION         PRESRATIl 

RISERS       TESTX         TESTY  USPOP         WARPbREAKS  VOLCANO 


)GRP  DISPLAY 
STEMLt'AF  STEMLEAFHOW  WIDTH         SCALE         LIT  hrxpr< 

bCXPLCTHCW  LENGTH       FILL  NUMSUM       Hl'hSUkRM  USS 

CR 


.7 


)GRP    DOC 

AN  DHOW       ANOVAHOW  A3RH0W  A3RSRH0W  BOXPLOTHOW 

CHISUUAREHOW  CODERESHOW  COMPAREHOW  CONDENSEHOW 

CONTINGENCY  HOW       CORRELAT  ION  HOW  DATAMAP    DESCRIBE  K.SHOW 

LINEtWW    LSLINEHOW  MEDPOLISHHOW  MSTATHOW  NUMSUMHOW 

ONEHHOW    0NE2H0W    PARTIALHOW  REGRESSHOW  SCAT  HOW    SHOWRSSHOa 

SPLITHOW  STATISTICSHOW  STEMLEAFHOW  TWICEHOW 
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)GRP    ESSENTIALS 

AND  ANOVA         A2R              A3RSR         bOXPLOT  CdlSQUARE 

COMPARE  CONDENSE                     CONTINGENCY  CORRfLATIi 

DEPTH  EPSILON    FILL            FMT              InPJT  KS                 Lb,                 '  ' 

IS  LI  NO.  tiEO  POLISH.                   hISS           NDIVX.  VOIVY 

toUilSUM  U:ihld            Otii3            OPTION       REGRESS  RESIOS       bCAt           ' SC 

S RON  RES  SPLIT         ST  EM  LE  At                      SUMSQ,  TJICb          JTC     .         ,v; 

Co.  MS  TATS      PARTIAL    STATISTICS  RS ELECT    R~             '*.PJ 


J       HI 


)GRP   GFTT 
KS  HSdON         CUISQ'JARE  CHISHUAREHOW 


)GRP    RELATIONS 

LliiE  LINEHCW    LSLINE  LS  LI'VEHC'N  REGRESS    REGRES. 

ANUdON      RS ELECT    Sr AT  SCATdCW     -lib  DEP  NDIVX 

^INTERCEPT  FMT  CR 


)GRP    SMOOTH 
SPLIT         0NE3  ONER  A3R  A3RSR         TWICE         SCAT  SPLIT RON 

A3RH0N       A3RSRR0N  0NE31WN    ONEHHOW    T'NICEHO'N  SCAT HO-/     vDI'/X 

NL)I\/Y         NIO  DEP 


)GRP    STATE 
CORRELATION  rh  FMT  MET  ATS      PAnriAL    STATIST 

CORRELATICNRCN       in  ST  ATS HON  PARTIALdCN  STAT 


)GRP    TNONAY 

MEDPOLISli  MEOPOLISHHON  LINE  LINE  RON    LSLINE       LSLINE  . 

RESIOS      EPSILON    NORMEFFECTS  SdONRES    CODERES    RSELECT    SdOWRESa     ■ 

COOERESdON  AND  ANDdON  CONTINGENCY  CONTJ..-     iCYdON       ANOVA 

AfiOVAdOW  OPTION       INPUT  SUmSQ 
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The  user  can  conserve  space  in  his  active  work  area  by 
selectively  loading  or  copying  only  the  required  functions  or 
groups.   We  have  attempted  to  anticipate  the  type  of  analyses 
that  the  user  will  attempt,  and  the  group  structure  has  been 
selected  with  the  objective  to  combine  functions  and  parameters 
which  naturally  go  hand-in-hand.   A  specific  group  can  be 
selected  by  typing 

)COPY  2  OA3660  GROUPNAME 

where  GROUPNAME  is  any  one  of  the  groups  listed  above.   Two  of 
the  groups,  ESSENTIALS  and  DOC,  are  especially  useful.   The 
ESSENTIALS  group  contains  all  of  the  functions  and  default 
values  of  all  of  the  required  parameters .   Documentation  and 
data  are  not  included.   The  DOC  group  contains  all  of  the 

documentation  (DESCRIBE,  DATAMAP,  and  "HOW"  variables). 
The  user  may  want  to  load  the  entire  OA366  0  workspace,  check 
on  some  of  the  documentation,  and  then  erase  the  group  DOC 
to  make  room  for  user  generated  data  and  other  functions  or 
variables.   The  DOC  group  can  be  erased  by  typing   ) ERASE  DOC. 
Similarly,  any  group  can  be  erased  by  typing  ) ERASE  GROUPNAME. 
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The  following  sections  give  brief  descriptions  (the 
writeups  contained  in  the  "HOW"  variables)  and  examples  of  the 
functions  contained  in  the  OA3660  library.   Each  function  descrip- 
tion gives  the  syntax,  a  list  of  the  parameters,  a  list  of  the 
groups  which  contain  the  function,  a  list  of  subroutines  (other 
functions)  used  by  the  function,  and  a  narrative  about  the  use 
of  the  function.   The  narrative  tells  what  the  function  does, 
describes  any  restrictions  on  the  function  arguments,  and 
indicates  how  the  function  parameters  affect  the  output.   The 
order  of  presentation  of  the  functions  follows  what  appears  to 
me  to  be  a  natural  sequence  of  exposure  to  tools  of  data 
analysis . 

The  HOW  variables  are  listed  alphabetically  in  Appendix  A, 
and  the  APL  program  listings  are  given  in  Appendix  B. 
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B.   Function  Descriptions  and  Examples 

)GRP  DISPLAY 

STEMS  AF  STEMLEAFHOW              7 I  DT  h         SCALE         LIT              bOXP 

BOXPLOTHOW  LENGTH      FILL           Nl'MSUM      NUldSUMHON                   '<iIS£ 
CR 


STEkLEAFHOW 
SYNTAXt       STEMLEAF    VECTOR 
GROUP',       DISPLAY 
PARAMETERS i 

(1)  WIDTH-    CONTROLS    Tut:    WIDTH    {CHARACTERS    PER    LINE)    Ob 
THE    DISPLAY    {DEFAULT =10) . 

(2)  SCALE-    VARIES    THE    DEPTH       {STEM      INTERVAL)       OF      Trii 
DISPLAY    IN    UNITS    OF    1,  2,  OR    0.5  TIMES    A    POWER    CF    10. 
SELECT    Ail    INTEGER    FROM    1    TO    3  ( DEFAVLT-1 )  . 

SUBPROGRAM t       LIT 

DESCRIPTION t  STEMLEAF  GENERATES  A  STEM  AND  LEAF  DISPLAY  OF  A 
VECTOR  OF  OBSERVATIONS.  THE  FUNCTION  AUTOMATICALLY  SCALES  T'i* 
DATA  USING  A  SCALING  ROUTINE  BASED  ON  THE  RANGE  AND  SIZE  OF  Thr 
DATA  VECTOR.  BY  VARYING  THE  SCALE  PARAMETER  THE  USER  CAN  CHANGE 
THE  SIZE  OF  THE  STEM,  BUT  IT  WILL  NOT  NECESSARILY  BE  SCALED  BY 
Tub:  AMOUNT  SPECIFIED.  TRY  VARIOUS  VALUES  LIKE  1,2,  AND  3  TO  EE*. 
WHICH  CHOICE  YIELDS  THE  BEST  DISPLAY.  IF  A  LEAF  CONTAINS  TOO 
MANY    CHARACTERS    TO    BE      PRINTED      On       A       LINE,       THE       LINE       '-JILL       bE 

truncated  and  the  number  of  truncated  characters  will  be  written 
at  the  end  of  the  line.  if  the  stem  interval  is  two  times  k 
power  of  ten,  the  leaves  may  contain  the  characters  a,8,c,  ...  ,* 
in  addition  to  the  digits  0,1,..,,  9,  in  this  case,  the 
characters  a  to  >■  represent  the  numbers   10  to   19. 


SCALE 
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STtinLclAV    VCLCAiiC 

00 

25666799 

01 

0  0  01 3 6 6 7  90 

02 

00112224445  56667  7  839  0  3 

Ob 

D 1 1 2  2  4  U  5  b  b  b  6  6  6  7  8  3  3 

0  4 
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0  6 

301144  5  5  6566 7 7  7  8  3 9 

0  7 

0  00011123  3 4555678889 

0  8 

12  2  2  2  3 3 3 b 6 7 9 

0  9 

0  301233'+  4  b  b  6  7  7  9 

10 

01122  3 3  4  4  5 639 

1  1 

0112  3 3  4  6  6  9 

1  2 

11244456 

1  j 

03478 

1  4 

1  5 

3  67 

1  6 

2  5 

17 

2  9 

13 

b 

19 

0  3  3  7  9 

a   I 
u  o  [ 

io| 

12| 

l6| 

Id! 


2  5  6  6  6  7 

0  0112  2 
01112 j 

3  0114  4+ 
12  2  2  2  3 

i;llZ^O 

1  1  z4U(4 
0  OGGh 

2  5r- 

•3  /i  l>  Z>  i'i  ' 


LA' a/'  VOLCAuC 

9  9  -m  A  /i  bUGGH'-   < 

244455666  7  7  8893  -jnBBCrtFk  FPc'GGGdl- 

3333hhH&7  8,69W}$AAtiti"rCDE£Fb,GGGGGGiih,i 

5  5  6  6  6  C  7  7  7  8  8  3 «  ri  AABBBC  Db f.  F F FG d  1  i  I  ~ 

3  3  5  6  7  3  A  A  ABC uD F &  t  FG  n  il ' 

344  5 6  BSAbbCbDFGG* 

b'&AUt1.  li  I 


R  C tiAuGIiiG  t'CALt.  :?' 
r  IdCtxFASIuG  SCALE 
fl    fii1    FAlit    'JUi-iBEh, 


'ul'Ci 


:1  ". 


<U.:LFh 


,/ILL    ALWAYS    KKFUL. 
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bO 
SYnTAA  \ 
RAnETE 


XPL0TH0W 

BOXFLOT    VECTOR 
h  : 

LENGTH*     CCNTRi 


•S    TtlE    HORIZONS  AL 


'IZE 


7/ ft" 


Dm, 


hi 


GROUP i 

ii  0  G  A 
hi  FT 
ATA. 
k  QU 
RISK 


SUBFi 
USSCi 
OF    Dt 
UP  Ft 
ASTEi 


(DEFAULT=bQ      CHARACTERS)  . 
DISPLAY 
Aki       FILL 
10 N\       BOXFLOT    GENERATES    A    bOX    PLOT    DISPLAY    F 

A    RECTANGULAR    BOX   WITH    ENDS    CORRESPONDING 
ARTILES    IS    PRESENTED      •■'/ITU      Tar.       kEDIAN      MAh 


ItiTE 
CROS 

1.5 


RQUA 

SES 
INTE 


WHISKERS    ARE    DRAWN       On      EACH      SIDE 
MARKING    THE    LOWEST    AND      HIGHEST      DATA 
UTILE    DISTANCE    OF    THE    QUARTILES.       LATA 
(OUTLIERS)    ARE    MARKED   W 'IT a    CIRCLES    AND 
RQUARTILE    DISTANCES    GET    HEAVY    CIRCLES, 


OF      THE 

VALUES 

VALUES 


A       V 

LOW  Eh 

u      W I T 

BOX. 

W I T  H 1 

OUTBID 

mvRE 


ft  r  L 

h 
WI 

i 

ft  ',': 
T  tl 


..  t\ 
lib 

hi. 
i  ,'i 

at 


BOXFLOT    VOLCANO 


xxxxxxxx I 


9  9 


[xxxxxxxxxxxxxx  GGCO  C  &$©& 
I  11 


a  WE    CAN    SQUEEZE    THE    BOXFLOT    IN    OR    STRETCH    IT    CUT    Li 
n  CHANGING    LENGTH. 


LENGTH+3 0 


BOXFLOT    VOLCANO 


99 


xxxxx     COOOSS® 
i  2  j  I 


LENGTH-*-!  Q 

BOXFLOT    VOLCAn 


199 


xxxxxxxxxxxl 


xxxxxxxxxxxxxx 


O  OO    0  0  »  **« 

3  c 
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hums  una  on 

SYNTAX  i       NUMSUM    VECTOR  OR 

PARAMETER i 

MISS*    NUMBER    US  to    TO    INDICATE    MISSING    VALUE    I..     ?ar.     UATA 

ARRAY    ( DEFAULT-' 99999) . 
GROUP'.      COMPARISONS 
SUBPROGRAM i      FMT 

DESCRIPTION  i       NUMSUM    OPERATES    ON    EITHER    A    VECTOR    OR    AN    ARRAY    BUT, 
IN    ALL    CASES,    THE    DATA    ARE    TREATED    AS    A      SINGLE      BATCH,  NUMSUM 

PRODUCES    A    NUMERICAL    SUMMARY    WHICH    GIVES    THE    SAMPLE      SIZt        {APT En 
DELETION    OF    MISSING       VALUES),       THE      EIGHTHS,       THE      EXTREMES,       Tttr. 
SPREADS,    AND    THE      MIDPOINTS    IN    TABULAR    FORM.         THIS      SUMMARY      IS 
USEFUL    FOR    TESTING    THE    SYMMETRY    OF    A    DATA    SET,       AND      TO      EV ALU ATI 
THE    EFFECTIVENESS    OF    DATA    TRANSFORMATIONS    IN    PRODUCING    SYMMETRY. 


NUMSUM    ARRAY 

MISSING    VALUE 


NUmSUi-j    VOLCANO 


r-  «*  p?  r?  n  t- 


NUMERICAL   SUMMARY 

SAMPLE    SIZE    =  219 

ul DPT S  L0Q/8/MIN                  MEDIAN           UPQ/B/MAX 

6  5.00!  I                  65.00|                                  I 

66.5  0|  3  7.00     I                                    I                   96.001                   59.00 

7  2.50J  2  u . 0  0     j                                   [                121.00|                   97.00 
10  0.50J  2.00|                                    I                199.00J                197.00 


SPREADS 


a  THE    STEttLEAF,    THE    BOXPLOT ,  AND    THE    NUMERICAL    SUMMARY    AL. 

p  SUGGEST    THAT    THE    VOLCANO    DISTRIBUTION    IS    POSITIVELY 

fl  SKEWED.       TRY    A    LOG    TRANSFORM   TO    SEE    IF    THAT    JILL    MAKr. 

a  THE    DISTRIBUTION    MORE    NEARLY    SYMMETRIC. 

NUMSUM    ®  VOLCANO 


NUMERICAL    SUMMARY 


'AmPLE    SIZt.     =  219 


MIDPTS 

4.17  I 

4.09  | 

3.99  j 

2.99  [ 


LOQ/8/MIN 

J,D1 

3.18 
0.69 


MEDIAN 

4.17 


UPQ/8/MAX 

! 
4.56  i 
4.80  [ 
5.29  1 


SPREADS 

0.95 
3  .62 

4.60 
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fl  THERE   IS    NOW    A    DECREASING    TREND    IN    THE   MIDPTS    COLUMN. 

fl  THAT    SUGGESTS   THAT    THE    LOG    IS   TOO    EXTREME    A    TRANSFORMATION . 

n  TRY    SQUARE    ROOT. 


NU14SUM    VOLCANO*.  5 


NUMERICAL    SUMMARY 


SAMPLE    SIZE 

= 

219 

irn!t»wnnii»inirniinw^«nn 

"nrrnniBPffirnnrt) 

i>»iTr,<-Tnr-"TTni«n 

MIDPTS 

LOQ/8/MIN 

MEDIAN 

UPQ/Q/MAX 

SPREADS 

8.06     | 

I                      8.06     | 

1 

! 

7.94     | 

6.08 

1                                     I 

9.80     | 

3.72     | 

7.95     | 

4.90 

1                                     1 

11.00     | 

6.10     i 

7.76    | 

1.41 

1                                      1 

14.11     | 

12.69     [ 

m!Uii«iiiii»niinii»mi»iin<iiiiriniiri?ririinn''iiiiitiiiirfi 


n»nr,",nrnn'IH'»rinrnT>"nrnri 


fl  THAT    LOOKS    A    LOT    BETTER.       LET    US    OBTAIN    THE    BOXPLOT 
fl  AND   STEMLEAF    PLOT    OF   THE    SQUARE    ROOT    TRANSFORMED    VAT  A. 


BOXPLOT    VOLCANO*. 5 


1.4 


14.1 
O     O    xxxxxxxxxxxxxxxxxx 


I xxxxxxxxxxxxxxxxxx  OOO 

I  2 


LENGTH+50 


BOXPLOT    VOLCANO*. 5 


1.4 


14.1 


O    O  xxxxxxxxxxxxxj 


[xxxxxxxxxxxxx  OO 


fl  THE    BOXPLOT    LOOKS    A    LOT   MORE   SYMMETRIC.       NOW    FOR    THE 
fl  STEMLEAF   PLOT t 
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SZEviLEAJ?    VOL 


i    M. 


0  *  .  5 


01 

02 
0  3 
04 
0  5 
06 
0  7 
08 

0  9 
10 
i  I 
i  2 

1  3 
14 


02 


00 

0  0 

0 1 

1 1 

!J  u 
-j  "| 

5  5 

46 
1 


4  6 
22  3 
4  45 

112 
122 


111 
ill 

1  1  -L 

in 

78 

899 


5 5 6 7  7  7  9  9  9 

2  3  3  4  i+  4  5  6  6  7  7  8  8  9  9  9  9 

2  3  4445666666668999 

1 1 1 2  2  2  3  3  3  u  4  5  5  5  5  5  5  5  5  5  7  7  7  7  3  t 

1122222344444445  5567  77  7888  39 

12  2  3  3  4  5  5  5  5  6  6  6  7  77  7  8  8  0  9 

1222344  5  5  566 G  7  8  8  9 

2  24  5  67788 
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NUM 

MISS            UTCOND 

FkT 

SC  at 

NDIVY 

WID              DEP 

DEPTH 

'1  G  A  P 

NUMSUM 

NUMSUMHOW 

R  SELF. CT 

and 

)GRP    COMPARISONS 
CONDENSE  CONDENSEHOW 

STAT  HOW    LIME  LINEUOW    NDIVX 

COMPARE    COMPARE HOW  FILL 

AH  DHOW      LSLINE      LSLINEHOW 


ANDHOW 

SYNTAX i       X       AND      Y 

GROUPS i       RELATIONS,    TWOWAY,    COMPARISONS 

DESCRIPTION :  THE    FUNCTION     'AND1    IS    USED    TO      CREATE       A       u EW       DATA 

ARRAY    CONSISTING    OF    X    AUGMENTED    BY    Y    AS    ADDITIONAL      COLUMNS.  X 

AND    Y    CAN    BE    SrALARS,     VECTORS,    OR    MATRICES      BUT      CANNOT       BOTh       BE 

SCALARS.       IF    AN    ARGUMENT    IS    A    SCALAR,     A    COLUMN    IS    GENERATED      EACH 

ELEMENT    OF    WHICH    IS    THE   SCALAR.       IF    THE    ARGUMENTS    ARE    VECTORS      OR 

MATRICES,    THEY    MUST    BE    CONFORMABLE. 


fi  THE    FUNCTION     'AND1    IS    USEFUL    FOhDATA    ENTRY    AND    MAN  I PU LA- 
n  TION. 


X+3    k    5 

Y+9    10  12 
X    AND    Y 
3    9 
U   10 
5   12 

Z<-3  2p\6 
X    AND    Z 

3   12 

u   3   u 

5   5   6 

5  0  AND    Z 

50    12 

5  0    3    4 

5  0    5    6 

5  0  AND    X 

50    3 

5  0    h 

50    5 


r  EXCEPT    WHEN    ONE    ARGUMENT    IS    A    SCALAR,    THE    ARGUMENTS 
a  AND    MUST    BE    CONFORMABLE    IN    SIZE. 

V«-4  6  8  10 

X    AND    W 

ARGUMENTS    OF      AND      ARE    NOT    CONFORM  At LE. 
W    AND    Z 

ARGUMENTS    OF      AND       ARE    NOT    CONFORMABLE. 
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COMPAREHOW 

SYNTAX i       COMPARE   MATRIX 
PARAMETERSx 

(1)  DEPTH*       VERTICAL      HEIGHT      OF      DISPLAY       (DEFAULT =  20 
LINES) . 

(2)  MISS"      NUMBER      USED      TO      INDICATE      MISSING       VALUES 
(DEFAULT    ="99999) 

(3)  NGAP»      NUMBER      OF      HORIZONTAL      SPACES      BETWEEN      THE 
BOXPLOT    DISPLAYS    (DEFAULT=3) . 

GROUP i       COMPARISONS 
SUBPROGRAMx    FILL 

DESCRIPTION t  COMPARE  OPERATES  ON  AN  N  BY  K  MATRIX  TO  PRODUCE  K 
VERTICAL  BOX  PLOTS  PLACED  NEXT  TO  EACH  OTHER  TO  ALLOW  VISUAL 
COMPARISON  OF  THE  CENTERS,  SPREADS,  AND  OUTLIERS  OF  THE  BATCHES 
(COLUMNS)  OF  THE  MATRIX.  THE  USER  MUST  FILL  UP  THE  MATRIX  SO 
THAT    THERE    ARE    N    OBSERVATIONS    FOR    EACH       BATCH.  THE      PARAMETER, 

MISS,  SHOULD  BE  USED  TO  FILL  IN  MISSING  VALUES.  THE  USER  SHOULD 
ASSURE  THAT  MISS  IS  DIFFERENT  FROM  ALL  VALID  ENTRIES  IN  THE  DATA 
ARRAY . 


DEPTH 
20 

MISS 
"99999 

NGAP 
3 

a  WE    WILL    DEMONSTRATE    COMPARE    WITH    THE   CHICKWT    DATA 
fl  SINCE    WE    USED    0    TO    INDICATE   MISSING    VALUES    IN    TEA'. 
a  DATA,    WE    MUST    FIRST    ALTER    THE   PARAMETER   MISS. 


MISS+O 
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COMPARE    CHICKWTS 
423 
108 


O  x 

_1 
x 

_-L_      x 
1*1     _1_ 

~t~    I  | 

I  ! 

-1-  t   I       T 

.1.     I  I  1*1 

I   I      1*1  ~T~ 

02    |*|    |  I     ®      * 
I  I    I  ( 
1  i    ~T~ 

.-!■_      T       x 
*  |      x  O 

t" 

x 


A  THE    DIAGRAM    REVEALS   RATHER    LARGE    DIFFERENCES    IN    THE 

a  MIDSPREADS    OF    THE    CHICKWT    DATA.       WE   SHOULD    TRANSFORM 

a  THE    DATA    TO    TRY    TO    ACHIEVE    HOMOGENIETY    IN    THE    SPREADS. 

a  FIRST    LETS    GET    A    NUMERICAL    SUMMARY    OF    THE    DIFFERENT    COLUMNS 


CONDENSEHOW 
SYNTAX  i       CONDENSE   MATRIX    OR   R+-CONDENSE   MATRIX 
PARAMETERS i 

(1)  NUMm    CONTROLS    WHAT    STATISTICS    ARE    INCLUDED    IN      THE 
SUMMARY.       NUMel 1,2, 5,71  {DEFAULT -2) . 

NUM=1    GIVES   MEDIANS   OF   EACH    COLUMN    OF    THE    MATRIX 
NUM-2    GIVES   MEDIANS    AND   INTERQUARTILE    RANGES 
NUM-b    GIVES    MIN,QUARTILES,MAX,    AND    SIZE 
NUM=1    GIVES   MIN, EIGHTHS ,MAX    AND    SIZE 

(2)  MISS*      NUMBER      USED      TO      CODE        MISSING         VALUES 
(DEFAULT=~99999) . 

GROUP t       COMPARISONS 
SUBPROGRAMS i       UTCOND,    FMT 

DESCRIPTIONx  CONDENSE  GENERATES  SUMMARY  STATISTICS  FOR  EACH 
GROUP  OF  DATA  REPRESENTED  BY  THE  COLUMNS  OF  THE  ARGUMENT  MATRIX. 
SINCE  THE  ARGUMENT  IS  A  MATRIX,  THE  USER  MUST  FILL  UP  THE  MISSING 
VALUES  IN  THE  MATRIX  WHENEVER  THE  GROUPS  (COLUMNS)  HAVE  DIFFERENT 
NUMBERS  OF  OBSERVATIONS.  THE  PARAMETER,  MISS,  SHOULD  BE  USED  TO 
FILL  THE  MATRIX  SINCE  THE  FUNCTION  WILL  RECOGNIZE  THOSE  VALUES  AS 
MISSING    DATA    AND   WILL   IGNORE    THEM   IN    ALL   CALCULATIONS.  THE 

RESULT    FROM    CONDENSE      CAN      BE      USED      AS      THE      ARGUMENT      OF      OTHER 
FUNCTIONS    SUCH    AS   SCAT    OR    LINE    IF    AND    ONLY    IF   NUM=2. 
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NUM 
MISS 


CONDENSE    CHICKWTS 

151.5  43 

221  83 

248  72 

328  33 

263  83 

342  102 


p  THE    ABOVE    DATA    ARE    THE   MEDIANS    AND    MIDSPREADS    EOR    EACH 

A  OF   THE    6  BATCHES    (COLUMNS)    OF    CHICKWT    DATA,       ZERO    WAS 

p  ENTERED    INTO    THE    PARAMETER   MISS    EARLIER,       THE    MIDSPREADS 

R  SHOW    QUITE    A    DISPERSION ,  FROM    AS    SMALL    AS    151.5  TO    AS    LARGE 

a  AS    342. 


A  WE    CAN    GET    MORE    INFORMATION    ABOUT    THE    INDIVIDUAL    COLUMNS 
a  OF    CHICKWTS    BY    CHANGING    NUM    AS    FOLLOWS s 


NUM+5 

CON  DEN 

SE    CHICKWTS 

MIN 

LOQ 

MEDIAN 

UPQ 

MAX 

SIZE 

108.00 

13  6.00 

151.50 

179.00 

227.00 

10.00 

141 .00 

175.00 

221.00 

2  5  8.50 

309.00 

12.00 

158.00 

199.00 

248.00 

271.00 

329.00 

14,00 

226.00 

307.50 

328.00 

340.50 

423.00 

12.00 

153.00 

242.00 

263.00 

325.00 

3  8  0.00 

11.00 

216.00 

2  71.50 

342.00 

3  7  3.50 

404. 00 

12.00 

a  IN    SEEKING    OUT    AN    APPROPRIATE    TRANSFORMATION    TO    EVEN    OUT 

a  THE   MIDSPREADS,    IT    IS    OF'l  EN    USEFUL    TO    PLOT    THE    LOG    OF    Tut. 

a  MIDSPREADS    VS.    THE    LOG    OF    THE    MEDIANS, 

A  IN    MANY    CASES   WE    NEED    TO    LOOK    AT    A    SCATTER    PLOT    OF    OUh 

a  DATA,    OR    AT    SCATTER    PLOTS    OF    FUNCTIONS    OF    OUR    DATA, 

a  LET    US   INTRODUCE    A    FUNCTION    SCAT    THAT    WILL    PROVIDE    THE 

A  NEEDED    SCATTER    PLOT, 


53 


SCATHOW 
SYNTAX i       SCAT    ARRAY 
PARAMETERS: 

(1)  WID«    CONTROLS    THE    HORIZONTAL    SIZE      OF      THE      DISPLAY 
{'DEFAULT- ZQ   CHARACTERS) . 

(2)  DEPv,    CONTROLS    THE      VERTICAL      SIZE      OF      THE      DISPLAY 
(DEFAULT=15    LINES), 

(3)  NDIVX,    NDIVY ~    NUMBER    OF    UNITS      ON      X*       AND      Y-AXES, 
RESPECTIVELY    ( DEFAULT** ,4) . 

GROUP i       RELATIONS,    SMOOTH,    COMPARISONS 

DESCRIPTION:  SCAT  PRODUCES  A  SCATTER  PLOT  OF  THE  DATA  CONTAINED 
IN  ARRAY.  THE  ARGUMENT  ARRAY  CAN  BE  A  VECTOR  OR  A  MATRIX  WITH  AS 
MANY  AS  9  COLUMNS.  IF  A  VECTOR  OF  SIZE  N,  THOSE  VALUES  ARE 
PLOTTED  VS.  THE  INTEGERS  1  TO  N;  IF  A  MATRIX,  THE  SECOND,  THlhD, 
ETC.  COLUMNS  ARE  PLOTTED  VS.  COLUMN  1  ON  THE  SANE  AXES. 
DIVISIONS  ON  THE  AXES  OF  THE  PLOT  ARE  NOT  EXPLICITLY  PRINTED, 
EXCEPT  AT  THE  EXTREMES  OF  THE  PLOT.  THE  USER  CAN  CONTROL  THE 
RESOLUTION  OF  THE  PLOT  BY  MODIFYING  THE  PARAMETERS  WID,  DEP, 
NDIVX    AND    NDIVY.  PRINTING      TIME      INCREASES      DRAMATICALLY      WITH 

RESOLUTION;  THEREFORE,  OTHER  PLOT  PROGRAMS  SHOULD  BE  USED  IF  HIGH 
RESOLUTION  IS  DESIRED.  FOR  A  SINGLE  GROUP  OF  DATA,  THE  NUMBER  N 
(2<N<,9)  WILL  BE  PRINTED  IF  N  POINTS  LIE  CLOSE  TOGETHER  ON  THE 
DISPLAY.  FOR  MULTIPLE  PLOTS  ON  THE  SAME  DISPLAY,  THE  LETTER  A 
REPRESENTS    GROUP    1,  B    REPRESENTS    GROUP    2,       ETC.  FOR      2  POINTS 

CLOSE   TOGETHER   THE    LETTER   WILL    BE   PRINTED    WITH      AN      UNDERSCORE. 
THE    DISPLAY    CANNOT    HANDLE    3  OR      MORE      POINTS      CLOSE      TOGETHER      IF 
THERE   ARE   MULTIPLE   PLOTS. 


WID 
30 


15 
4 

4 


DEP 

NDIVX 

NDIVY 


A>il0 

Yl«-2  3  5  6  8  10  11  13  15  18 
Y2+-20  20  17  24  21  26  31  30  37  33 
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SCAT    X    AND    Yl 
RANGE    OF    Xi     0   10 
RANGE    OF    It  0   20 


o 
o 
o 
o 
o 
o 


SCAT    X    AND    Yl    AND    12 
RANGE    OF    Xt     0       10 
RANGE    OF    Y\     0       4  0 

B    b 
B  B 

B    B  B  A 

B  A       A 

AAA 
AAA 


q  IF  2*#/T  ARGUMENT    OF    SCAT    IS    A    VECTOR,    SCAT    'JILL    CREATE 

n  4B5CI55/1  VALVES   TO    BE   THE    POSITIVE   INTEGERS    1  7"  pARGUtiENT 

SCAT    Y2 
RANGE    OF    Xi     0   10 
fiMGF.  5F  Is  15   4  0 


o   o  o 

o 
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a  BACK   TO   THE   PROBLEM   OF    NONHOMCGENEOUS   MIDSPREADS    IN    THE 
A  CHICKWT    DATA.       LET    US    LOOK    AT    A    SCATTER    PLOT    OF    THE 
a  LOG,  MIDSPREADS    V.    LOG    MEDIANS. 

a  FIRST    WE    RESET   PARAMETER   HUM    TO    2,  SO    IT    WILL    GIVE 
a  MEDIANS    AND   MIDSPREADS . 

NUM+2 

SCAT    OCONDENSE   CHICKWTS 
RANGE    OF    Xi    5   6 
RANGE    OF    Yi     3.4   4.8 


o     o 
o 


o 


A  THE    SCATTER   PLOT    SUGGESTS    A    RELATIONSHIP    BETWEEN    MIDSPREADS 
A  AND   MEDIANS.       PERHAPS   SOME    ROOT    OF    THE    DATA    WOULD    MAKE 
a  THE   MIDSPREADS   MORE    HOMOGENEOUS.       TRY    SQUARE    ROOT. 

CONDENSE   CHICKWTS*. 5 

12.30368569  1.717184371 

14.86363273  2.8510  5  5  496 

15.74801575  2.355341653 

18.11001266  0.9195158937 

16.21727474  2.471407191 

18.4912651  2.852064545 

a  THAT    IS    BETTER    (RATIO    OF   MAX   MIDSPREAD    TO    MIN    MIDSPREAD 
a  IS    ABOUT    1.5.  MORE   WORK    NEEDS   TO    BE    DONE   THOUGH. 
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)GRP   RELATIONS 

LINE  LINE HOW    LSLINE  LSLINEHOW  REGRESS    REGRESS  HOW  AND 

AN  DHOW      RSELECT    SCAT  SCATHOW    WID  DEP  NDIVX         NDIVY 

^INTERCEPT  BUT  OR 


LINEHOW 
SYNTAX:       LINE    ARRAY  OR  Z+LINE    ARRAY 

PARAMETER: 

REELECT*    SELECTS      THE      OUTPUT      MATRIX      OF      RESIDUALS. 

RSELECT =1    GIVES    ABSCISSA    VS.       RESIDUALS       {DEFAULT-!) . 

RSELECT*!    GIVES    FITTED    VALUES    VS.    RESIDUALS 
GROUP i      RELATIONS 

DESCRIPTION:       LINE    FITS    A    STRAIGHT    LINE    TO    A    SET    OF    (X,Y)       POI.\ 
BY    DIVIDING    THE   POINTS    INTO    3  REGIONS    AND   USING      THE      MEDIANS       OF 
THE   X    AND    Y    VALUES    IN    THE    OUTER    REGIONS    TO    DETERMINE    THE    SLOPE. 
THE    INTERCEPT    IS    THE   MEDIAN    OF    THE    DIFFERENCES    Y*SLOPE*X.  LINl 

ALSO  COMPUTES  THE  RESIDUALS  AND  GIVES  AN  N  BY  2  MATRIX  CONTAINING 
EITHER  THE  XvVALUES  VS.  RESIDUALS  OR  THE  FITTEu  VALUES  VS. 
RESIDUALS  DEPENDING  ON  THE  PARAMETER,  RSELECT.  IF  THE  USER  DOES 
NOT    WANT    RESIDUALS    TYPED    OUT,    HE    SHOULD    USE    THE    SYNTAX: 

INLINE  ARRAY 
THE  ARGUMENT  ARRAY  CAN  BE  EITHER  A  VECTOR  OR  A  2«CQLC'MN  MATRIX. 
IF  THE  ARGUMENT  IS  A  VECTOR  OF  SIZE  N,  THE  Xr,VARIABLE  IS 
CONSTRUCTED  TO  BE  THE  FIRST  N  POSITIVE  INTEGERS.  IF  THE  ARGUMENl 
IS  A  IvCCLUMN  MATRIX,  THE  FIRST  COLUMN  IS  TAKEN  TO  BE  THE  SET  OF 
XnVALUES  AND  THE  SECOND  COLUMN  THE  SET  OF  Y~VALUEE.  LINE  CAlt  BE 
USED  IN  CONJUNCTION  WITH  THE  FUNCTION  SCAT  TO  PRODUCE  THE  SLOPE, 
THE    INTERCEPT,    AND    A    PLOT    OF    THE    RESIDUALS.       TO    DO    THIS,    ENTER: 

SCAT    LINE    ARRAY 
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X+2    3435676759810111211+ 
pX 


16 


SLOPE 

2 

3 

4 

3 

5 

6 

7 

6 

7 

5 

9 

8 

10 

11 

12 

14 


SLOPE 

2 

4 

5 

4 

6 

8 

9 

8 

9 

6 

12 

11 

13 

15 

16 

19 


SLOPE 
RANGE 
RANGE 


ARRAY+X    AND    Y 
LINE    ARRAY 
1.375  Y'INTERCEPT 


RSELECT+Q 
LINE   ARRAY 

1.3-7  5  Y  ^INTERCEPT 

8125  0.1 

1875  "0.1 

5625  "1.5 

1875  0.8 

9375  1.0 

3125  "2.3 

6875  "5.6 

3125  "0.3 

6875  "0.6 

9375  3.0 

4375  "0.4 

0625  3.9 

8125  2.1 

1875  "1.1 

5625  1.4 

3125  0.6 
RSELECT+l 
SCAT    LINE    ARRAY 
:  1.375  Y'INTERCEPT 
OF   X:    2   14 
OF   Y:    "6   4 


:  0. 

875 

875 

625 

125 

625 

125 

875 

125 

875 

625 

375 

375 

875 

875 

375 

875 


:  0 
875 
875 
625 
125 
625 
125 
875 
125 
875 
625 
375 
375 
875 
875 
375 
875 


0625 


0625 


:  0.0625 
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3   4 

16 

SLOPE 
1 
2 
3 
4 
5 
6 


9 
10 
11 
12 
13 
14 
15 
16 


SLOPE 
1 
2 
3 
4 
5 
6 
7 

8 
9 
10 
11 
12 
13 
14 
15 
16 


Y 

4   5   8   6 

PY 

LINE     Y 
1. 090909091 


10 


12   15   16   14   13   20 


Y~INTE 
1.  227 


136 
045 
045 
863 
22^ 

3  18 

4  09 
5 

590 
681 
22" 
136 
954 
954 
863 


X+\  16 

ARRAY ^X    AND    . 
LI i)l E    ARRAY 
1. 090909091 


Y-INTE 
1.22'7 


136 

045 

045 

863 

22^ 

3  18 

409 

5 

590 

681 

227 

136 

9  54 

954 

863 


RCEPT:     0.6818181818 
2  n  2  i  9  n 

363636 

45454545 

45454545 

636364 

27272'7 

181818 

090909 

909091 

8131818 

27  27  2  7 

363636 

545455 

5454545 

636364 


RCEPT:     0.6818181313 

2  7  2  7  2  7 

363636 

45454545 

45454545 

636364 

27  272? 

181818 

0  9  0  9  0  9 

909091 
8181818 

o  t  o  n  o  7 

363636 
545455 
5454545 
636364 


SLOPE'.     1.090909091  Y  -I  NT  a  RC  E  PT  :  0.6818181818 
RAHG&    OF    X:    0       2  0 
RAUGE    OF    Y :     "6   2 


t»  •*   t»  ■» 


1  T  5  7  »  ?  7 


r»  t*  ■?»  ■»  ■»  •»  t* 
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RSELECT+1 
Z-LINE    Yi 
SLOPE :  1.857142857  YvINTERCEPTi     17.4285  7143 


SYNTAX i       LSLINE    ARRAY  OR  Z-LSLINE    ARRAY 

GROUP i       RELATIONS 

DESCRIPTION i  LSLINE  DETERMINES  THE  LEAST  SQUARES  SOLUTION, 
Y=YI+SLOFE*X.  THE  ARGUMENT  ARRAY  CAN  BE  A  VECTOR  Oh  A  2-CCLUMh 
MATRIX.  IF  THE  ARGUMENT  IS  A  VECTOR,  THE  XvVALUES  ARE  T AKFJi  TO 
BE    THE    INTEGERS    1  TO    N    WHERE    N=p ARRAY.  IE      THE       ARGUMENT      IS       A 

2*C0LUMN  MATRIX,  COLUMN  1  CONSISTS  OF  THE  Xr-VALUES  AND  COLUMN  2 
THE  Y ^VALUES.  RESIDUALS  ARE  DETERMINED  AND  STORED  IN  A  lATRIX 
WHOSE  FIRST  COLUMN  IS  THE  SET  OF  XvVALUES  AND  WHOSE  SECOND  COLUMh 
IS  THE  SET  OF  RESIDUALS.  IF  THE  USER  WANTS  TO  SUPPRESS  PRINTOUT 
OF  THE  RESIDUALS,  HE  MUST  ASSIGN  THE  RESULTS  OF  LSLINE  TO  A 
VARIABLE;    I.E.,    TYPE    Z+LSLINE    ARRAY. 


SCAT    LSLINE    X  AND    Y2 

SLOPE i     1.933333333  Y rj  NT ERC EPT :  16.8  6  666667 
ANGE    OF    X\     0   10 
ANGE    OF   It    "4   u 


o       o 


o 


o 


o 


nIF  THE  USER  WANTS  TO  REGRESS  A  RESPONSE  VARIABLE  ON 
n  /  SET  OF  CARRIERS,  HE  CAN  USE  THE  PROGRAM  REGRESS. 
n  THIS  LEAST  SQUARES  MULTIPLE  REGRESSION 
R  VERY    POWERFUL    AND    FLEXIBLE. 
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REGRESSHOW 
SYNTAX i       Z+Y    REGRESS   X 
PARAMETER: 

^INTERCEPT*    DETERMINES   WHETHER    OR    NOT    AN    INTERCEPT    TERM 
IS    TO    BE    INCLUDED.         LINTERCEPT-1      GIVES      AN      INTERCEPT 
TERM,    AND    MNTERCEPT=0    GIVES    NO    INTERCEPT.     (DEFAULT      IS 
1.) 
GRCUPx       RELATIONS 
SUBPROGRAMS \       FMT    AND    SCAT 

DESCRIPTION t  REGRESS  DOES  A  MULTIPLE  REGRESSION  ANALYSIS 
RELATING  THE  DEPENDENT  VARIABLE  Y  TO  A  SET  OF  CARRIERS  X.  THE 
LEFT  ARGUMENT  Y  IS  A  VECTOR  OF  SIZE  N.  THE  RIGHT  ARGUMENT  X  IS 
AN  N  BY  K  MATRIX  CONSISTING  OF  N  OBSERVATIONS  ON  EACH  OF  K 
VARIABLES  OR  A  VECTOR  OF  SIZE  N  IF  K-l ♦  OUTPUT  CONSISTS  OF  Ah 
ANOVA  TABLE,  RnSQUARE,  STD.  ERROR,  REGRESSION  COEFFICIENTS  (THE 
FIRST  COEFFICIENT  IS  THE  CONSTANT  TERM  IF  MNTERC EPT=1 .) , 
TmSTATISTICS,  V ARIANCEnCOV ARI ANCE  MATRIX,  DURBIN^W ATSCN 

STATISTIC,  AND  A  VECTOR  OF  PREDICTED  Y  VALUES  AND  RESIDUALS. 
THERE  IS  AN  OPTION  THAT  ALLOWS  THE  USER  TO  INPUT  A  VECTOR  OF  X 
VALUES  AND  USE  THE  REGRESSION  EQUATION  TO  FORECAST  Y  VALUES.  THE 
USER    CAN    ALSO    OBTAIN    A    SCATTER      PLOT      OF      THE      RESIDUALS.  tfHEti 

EXECUTION  TERMINATES,  THE  PREDICTED  Y  VALUES  AND  THE  RESIDUALS 
RESIDE   IN    THE    N    BY    2    MATRIX    Z. 

^INTERCEPT 


pDRAPERX 
13   4 

*LETS   REGRESS    A    RESPONSE    VARIABLE    DRAPERY    ON    THE    FOUR 
n  CARRIERS    IN    THE   MATRIX    DRAPERX. 
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Z+DRAPERY    REGRESS    DRAPERX 


ANCVA 


SOURCE 

DF 

SUM   SQUARES 

MEAN    SQUARE 

FnRATIO 

REGRESSION 

4 

2.6679F3 

6.6697£2 

1.1148E 2 

RESIDUAL 

8 

4.7864£1 

5.9830£0 

TOTAL 

12 

2.7158£3 

R    SQUAREx       0.9823756203 
STD    ERRORx     2.446007963 

COEFFICIENTS         T    STATISTICS 


62.4054  0, 

1.5511  2, 

0.5102  0, 

0.1019  0. 

~0.1441  "0. 

YOU    WANT    A    PRINTOUT  OF 


DO 

Y 

V ARI ANCEwCOV ARI ANCE   MATRIX i 


8906 

0827 

7049 

135 

2032 

THE   V ARI ANCE «CQV ARI ANCE 


MATRIX? 


4.9099£+3  5 
~5.0507ff+l  5 
"5.0603E+1  5 
'5.i660£+l  5 
~4.9597£-*-l    5 

DURBINvWATSONi     2 

DO   YOU    WANT    TO 

20  30  15  40 

DC    YOU    WAN1 

Y 

RANGE    OF    Xt    7  0 

RANGE    OF    Yt  ~4 


0507£+l 
5468£~1 
1266£~1 
5425£~1 
0529£~1 
052596933 


0603£+l 
1266£~1 
2387£~1 
2570£~1 
1213£~1 


1660£+1 
5425£~1 
2570f"l 
6  9  59A'"l 
1688E""l 


FORECAST    A    VALUE    FOR    Y? 


TO    SCAT    RESIDUALS    VS.    PREDICTED    Y? 


120 
4 


959 
052 
121 
168 
027 


7&>1 
9/T  3 
3E~1 
8E~1 
5E~1 


f»    -»nOwntTf,r-n^»-Owr,r,^rrr.T.n^Or.f. 
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)GRP    SMOOTH 
SPLIT         0NE3  ON  EH  A3R  A3RSR         TWICE         SCAT  SPLIT  HOW 

A3RH0W       A3RSRH0W  0NE3H0W    QNEHHOW    TWICEHOW  SCATHOW    tiDIVX 

NDIVY         WID  DEP 


0NE3H0W 
SYNTAX i       0NE3    VECTOR 
GROUP i  SMOOTH 

DESCRIPTION  i  CNE3  SMOOTHS  A  VECTOR  OF  DATA  USING  ONE  PASS  OF 
RUNNING  MEDIANS  OF  3  WITH  TUKEYyS  ENDvPOINT  RULE.  WHEN  SMOOTHING 
(X,Y)  PAIRS ,  THE  USER  SHOULD  USE  AS  THE  ARGUMENT  VECTOR  THE  Y 
VALUES    ORDERED    ACCORDING    TO    THE   MAGNITUDE    OF    THE    X    VALUES t 

VECTOR+YllX] 
IT   IS    ASSUMED      THAT      THE      X      VALUES      ARE      EQUISPACED.  REPEA 

SMOOTHING    BY    MEDIANS    OF    3    CAN    BE    ACCOMPLISHED   WITH    A3R.  OUTPUT 

CONSISTS    OF    THE    SMOOTHED      SEQUENCE.  A      PLOT      OF      THE 

SEQUENCE   IS    GIVEN    BY    TYPING \ 

SCAT    0NE3    VECTOR 


X 
3237432   15   4536 

0NE3    X 
233443345453 

SPLITHOW 
SYNTAX i  SPLIT    A3R    VECTOR 
GROUP i       SMOOTH 

DESCRIPTION i  SPLIT  DOES  ONE  PASS  AT  DIVIDING  MESAS  (PAIRS  OF 
ADJACENT  POINTS  WITH  A  COMMON  VALUE  WHICH  IS  A  LOCAL  MAX  OR  MIN) 
USING  TUKEY'S  END^POINT  RULE.  IT  IS  USED  IN  CONJUNCTION  WITH 
A3R. 


SPLIT    ONE 3    X 
3   3   3   3   4   3 


A3RHCW 
SYNTAX i       A3R    VECTOR 
GROUP i       SMOOTH 
SUBPROGRAM i       0NE3 

DESCRIPTIONt  A3R  DOES  REPEATED  SM0CTH1NGS  OF  RUNNING  MEDIANS  OF  3 
UNTIL  THERE  ARE  NO  CHANGES  IN  THE  SMOOTHED  SEQUENCE  FROM  ONE 
ITERATION    TO    THE    NEXT.       SEE    CNE3H0W . 
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A3R   X 
3   3   4   4 


A3RSRH0W 
SYNTAX  i   4  3/?Stf  VECTOR 
GROUP i      SMOOTH 

SUBPROGRAMSt       A3R ,    0NE3,    SPLIT 

DESCRIPTION t       A3RSR    DOES    REPEATED    SMOOTHINGS    BY      RUNNING      MEDIANS 
OF    3    FOLLOWED    BY    SPLITTING    MESAS    AND    REPEATING    UNTIL    CONVERGENCE. 
SEE    A3R    AND   SPLIT. 


A3RSR   X 
3   3   3   3   4 


ONEHHOW 
SYNTAX i       ON EH    VECTOR 
GROUP t'    SMOOTH 

DESCRIPTION i       ONEH    DOES    A    SINGLE  'MANNING    OF    A    SEQUENCE    OF    DATA. 
IT    IS   USED   IN    CONJUNCTION   WITH    THE    OTHER   SMOOTHING    FUNCTIONS    AS    A 
FINAL    TOUCHwUP   TO    A    SEQUENCE    OF    SMOOTHED    DATA.       THE   I«TH    RESPONSE 
IS   REPLACED    BY    0 . 25*7 [Inl ]  +  0.5*YLI'\    +  0. 2 5xY[ 1+1 ] .  SEE    A3R    AND 
A3RSR. 


ONEH    X 
2.5   3.75   5.25   4.5   3   5.5   9   7   4.25   4.25 


TWICEHOW 
SYNTAX i       TWICE    VECTOR 
GROUP i       SMOOTH 

SUBPROGRAMSt       ONE3,    SPLIT,    A3RSR 

DESCRIPTICNx  TWICE  SMOOTHS  A  SEQUENCE  OF  DATA,  THEN  SMOOTHS  THE 
RESIDUALS  AND  ADDS  THE  SMOOTHED  RESIDUALS  BACK  TO  THE  SMOOTHED 
DATA  TO  OBTAIN  THE  FINAL  SMOOTH.  THE  SMOOTHING  IS  DONE  USING 
A3RSR. 
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TWICE   X 

3   3   3 


RlW//    LET    US    LOOK    AT    SOME    PLOTS    WITH    a0,-.i    !  ORE    InTER 


n  DATA    THAT    SHOWS    THE    EFFECi 


OF    SMOOTHING. 


f\WE    JILL    USE    THE    PRE  SEATING    DATA. 


W ID-US 
DEP-2D 

SCAT    PRE SEATING 
kA'aGF.    OF    X  i  0   120 
KAdiGE   OF    Yi    0   100 


o 

o  o 


o  o 

2 


2  °  2  °  2  °    ° 

0  0         o 

°  2  °    3  °  2  °  3 

o 

2 


o 
o  o  o 

°  2 


2 
o  o 


o  I   o  °  ° 

0 

2 


o 


3  2 


nIT    IS    HARD    TO    TELL    tiUrd    FROk    THIS, 
r  HELPS    AMY. 


LETS    SEE    IF 


SO  AT    A3RSh    PRESRATIhG 
RAUGE    OF    X\     0   120 
RAHGE    OF    It     20   90 


0  z 

o 


2  2  3  3  3° 

o 

o  o 

2 


o  2  ° 
23 
233°  o     22  °°23° 

oo   oo    2°2     °2 

O  2    °  2  0        0  0 

o  o  o         2 

2°  °322 

o  2  3  3  °      °  ° 
2 


§2^ 
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nTHAT   IS    A    LOT   SMOOTHER.      THE   PEAKS    Ai^U    VALLEYS    ARE    A    LOT 
r  CLEARER    SOW  *       AS    A    FINAL    TOUCHUP,    LETS    USE   TWICE    TO 
n  SMOOTH    THE    RESIDUALS    AND    ADD    BACK    TC    THE    SbiCCTHED    SEQUh 
n  AND   THEN    USE    OSEH    AS    A    FINAL    POLISH. 


SCAT    ON  Eli    TWICE    PRESRATISG 
RANGE    OF   Xt    0   12  0 
RANGE    OF    Yt  20   90 


°2  2 

23  2     2 

233°     °  2°     °  °  2  3  ° 

o  003°  2°2  °2 

2      o  o  o 

2223°  o  o  2  o  o     2 

oo  2  °  °     °  3  2  2 

ooo  233  o° 

2       2°  ° 

o  o 

o       o  o 

°2°  22 
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)GRP    TWOWAY 


MEDPOLISH 
RESIDS      EPSILON 
CCDERESHOW 
ANOVAHOW 


MEDPCLISHHCW 
"v  c 


And 
OPTION 


ANDHOW 
INPUT 


LINE  LINEHOW 

SHOWRES    CODERES 
CONTINGENCY 

SUMSQ 


LSLINE      LSLINEHOW 
RSELf.'C?    SHOWRESHCW 
CONTINGENCY  HO  j       AUC\f  A 


Z+MEDPCLISH    MATRIX 


MEDPCLISHHCW 
SYNTAX i       MEDPOLISH    MATRIX  OR 

PARAMETERSx 

(1)  RESIDS*  CONTROLS  OUTPUT  OF  RE 
RESIDUALS  ARE  PRODUCED  IN  A  TWO-WAY 
COMPARISON  VALUES,  (RE° .*CE)iTV  AN 
PRODUCED    IN    A    TWO    COLUMN   MATRIX    (DEFA 

(2)  EPSILON m       PROPORTION      BY       WHICH 
ABSOLUTE    VALUES    OF   THE    RESIDUALS      MI'S 
EACH    ITERATION    TO    CONTINUE    POLISHING 

(3)  NORMEFFECTSn      DETERMINES   WHETHER 
EFFECTS    ARE    OUTPUT,       THE    DEFAULT    OF    0 
OF   NORMALIZED    VALUES.       N0RMEFFECTS=1 
VALUES   TO    BE    PRINTED. 

GROUP i-     TWOWAY 

DESCRIPTION j  MEDPOLISH  ITER ATIV ELY  SWEEPS  OUT 
ROWS  AND  COLUMNS  OF  A  TWO-WAY  TABLE  TO  YIELD  TH 
OBS  -  MEDIAN  +  ROW  EFFECT  ^  COLUMN  EFFECT 
IT  YIELDS  RESIDUALS  FOR  TESTING  THE  ADEQUACY  OF 
SHOWRES  AND  CODE RES).  THE  FUNCTION  CONTINUES  PC 
STOPPING  RULE  CONTROLLED  BY  EPSILON  IS  ACTIVATE 
MEDIAN,  THE  ROW*  AND  COLUMN -EFFECTS ,  AND  THE 
VALUES   OF    THE      RESIDUALS      ARE      GIVEN.  THE      U 

RESIDUAL   PRINTOUT    BY    ASSIGNING    THE    FUNCTION    TO 
BY    WRITING i 

Z+   MEDPOLISH    MATRIX 


SI DUALS.  IF 

TABLE ;     IF 
D      RESIDUALS 
ULT=1) . 

THE      SUM      OF 
BE      REDUCED 
(DEFAULT =0.01) 
OR    NOT       NORMAL 

SUPPRESSES    OU 
CAUSES      NORMAL 


-1  » 

2, 

A  h  E 

THE 

AT 

* 

I  ZED 

TPUT 

'ZED 


MEDIANS      FROm      THE 
E    MODEL; 

-r  RESIDUAL 
THE      MODEL 
LI SUING    UNTIL 
D.  THE      CVE 

SUM  OF  ABSC 
SEE  CAN  SUPP 
A    VARIABLE,       I 


(SEE 

Tut. 
HALL 
LUTE 
Rr'SS 
.E.   , 


1 
0.01 

0 


RESIDS 

EPSILON 

NORMEFFECTS 
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DEATHS 

11.7 

8.7 

18.1 

11.7 

26.9 

20.3 

41 

30.9 

66 

54.3 

15 
24 
37 
5U 
71 


8.4 
13.6 
19.3 
35.1 
50 


Z+MEDPOLISH    DEATHS 
315.4 
40.3 
38.25 
37.575 

TYPICAL    VALUEt     25.6 
ROW    EFFECTS I  "15.3625   "10.1 
COLUMN   EFFECTS t    2.7   ~3.9   12 


1.4 
625 


12.6375 
3.1375 


32.875 


fiTHE   RESIDUALS    ARE    STORED   IN    Z.       THE    FIRST    VALUES 
n  PRINTED   OUT    ARE   THE   SUMS   OF   THE   ABSOLUTE    VALUES 
n  OF   THE   RESIDUALS.    ITERATION   CEASES   WHEN   TWO 
ft  CONSECUTIVE    VALUES    ARE    WITHIN    EPSILON. 


Z 
1.2375 

0.1 
0 

0.0625 
4.825 


2.3625 

0.1 

0 

3.4375 

0.275 


7.4625 
'3.825 
0.175 
3.7375 

0 


3 

2375 

7625 

3375 


37.2375 


SHOWRESHOW 
SYNTAX t       SHOW  RES   MATRIX 
SUBPROGRAMt       AND 
GROUPt       TWOWAY 
DESCRIPTIONx       SHOW  RES    AIDS 
MEDIAN   POLISHING    A    TWOmWAY 
IN    MEDPOLISH   IF    SHOWRES   IS 
A    MATRIX    CONSISTING    OF    THE 
OF   THE   RESIDUALS. 
OF    ALL    RESIDUALS, 


OR 


SHOWRES   MEDPOLISH   MATRIX 


IN  THE  ANALYSIS  OF  RESIDUALS  FROM 
TABLE.  THE  USER  SHOULD  SET  RESIDS=1 
TO  BE  USED.  THE  OUTPUT  OF  SHOWRES  IS 
SYMBOLS  COoxX  DEPENDING  ON  THE  SIZES 
LET  R  REPRESENT  THE  RESIDUAL,  M  THE  MIDSPREAD 
L   THE   LOWER    QUARTILE,    AND   U    THE    UPPER   QUARTILE. 


THE   CODES    ARE   PRINTED    AS   FOLLCWSx 
Oi    R<Lnl.5M 

Lwl.$M<,R<>LnM 

LmM<R<U-rM 

U+MZRZU+1.5M 
R>U+1.SM 


Oi 

«  I 

Xi 
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SECURES 

o 
o 
o 
o 


CODERES HOW 
SYNTAX i       CODERES    MATRIX  OR 

SUBPROGRAM :  AilD 
GROUP i      TWO J  AY 

DESCRIPTION i       CODERES    PRODUCES    A    DIAGNO 
THE    RESIDUALS    FRCk    MEDIAN    POLISHING   OF 
SHOULD   SET    RESIDS-1    In    MED POLISH      WHEN 
OUTPUT    MATRIX    CONSISTS    OF    THE    SYMBOLS    - 
OF    THE    RESIDUALS. 

-:  VALUE    BELOW    LOWER    QUAhTIL 
01  VALUE    BETWEEN    QUART I LEE 
+  t   VALUf.    LARGER    THAN    UPPER    q 


CODE  SEES    ME  DP  0  LIS  ii     -iA 


STIC 
A    TWO 
CODE 
o-t-  DE 

E 

UARTI 


ARRAY 
-WAY 
Rt<;E 
PEiJDI 


LES 


ten 
7  Ah  Lb. 

IS      U 


/  i,  Hi 

SR  u . 


Y'L'x 


H  c 


f 


CODE  RES    Z 


fi    LETS    SEE    WHAT    HAPPENS    WHEN    ht:FIuS- 


RESIuS+2 


Z+-MEDPCLISH    DEATHS 

315.  4 

40.  3 

38.25 

3  7.575 

TYPICAL    VALUE i     2  5.6 

h OW    EFFECTS i  "15.3  625   " 1 0 

COLUMN    EFFECTS i     2.7   "3.9 


j 
12 


1.4 
625 


12 
3.; 


6375 
375 


32.875 
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SCAT    Z 

RANGE    OF    Xi 

RANGE    OF    Yi 
1 

"10 
"8 

20 
6 

1   '• 

1 
1 

1   o 

o 

|  o  o 

*fn<frr°Wv322ffR'*0«tr! 

•1 

o 

r'-ti^iOfTWf 

•  1 

1 

o     1 

A  THERE    SEEMS   TO    BE    A    RELATIONSHIP    BETWEEN    THE   RESIDUALS 
fi    AND   THE   COMPARISON    VALUES,     (RE* .*CE) *TV . 


A+LINE    Z 
SLOPEi     0.6216761995  Y~INTERCEPTi     "0.1348  011389 


ANOVAHOW 
SYNTAX i  ANOVA 
PARAMETER i 

OPTION m    SELECTS   THE   TYPE    OF    ANOVA    PERFORMED.         OPTION =1 
GIVES    A    ONE-WAY    ANOVA      AND      0PTI0N=2      GIVES      A      TWO-WAY 
ANOVA    (TREATMENTS   WITH    BLOCKING).       (DEFAULTED 
GROUPi       TWOWAY 

SUBPROGRAMSx       INPUT,    SUMSQ,    FMT 

DESCRIPTION  i  ANOVA  DOES  A  ONEWAY  OR  A  TWQvWAY  ANALYSIS  OF 
VARIANCE  DEPENDING  ON  THE  VALUE  OF  THE  OPTION  PARAMETER.  ANOVA 
WILL  INTERACT  WITH  THE  USER  TO  OBTAIN  THE  REQUIRED  INFORMATION 
AND  DATA.  WHEN  ENTERING  DATA,  SEPARATE  THE  DATA  POINTS  WITH  AT 
LEAST  ONE  BLANK  AND  USE  NO  OTHER  DELIMITERS.  OUTPUT  CONSISTS  OF 
AN  ANOVA  TABLE,  ESTIMATES  OF  THE  OVERALL  MEAN,  AND  ESTIMATES  OF 
THE   TREATMENT    (AND, IF   OPTION-!,    BLOCK)    EFFECTS. 
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OPTION 


ANCVA 
ENTER    NUMBER    OF    TREATMENTS. 


enter  ieotor  of  number  of  cbs.  for  treatments   i  to   u 

J: 

3  3  3  3 

ENTER    3  OBSERVATIONS    FOR    TREATMENT    1 
J: 

4  5  4  6  51 

ENTER    3  OBSERVATIONS    FOR    TREATMENT    2 
0  i 

4 2  4  4  5  0 
ENTER    3  OBSERVATIONS    FOR    TREATMENT    3 

G « 

36  41  U3 
ENTER    3  OBSERVATIONS    FOR    TREATMENT    u 

U: 


4  9  4  7 

54 

SOURCE 

TREATMENT 

ERROR 
TOTAL 

3 

8 

11 

R "SQUARE    -    3.419 
OVERALL    MEAN    = 
TREATMENT    EFFECTS 


An  OVA    TABLE 

Q  o 

110,92 
154,00 
264. 92 


3  6.97 
19.25 


4  6.08 

1.2  5   0.7  5  "4.U2   3.92 


1.32 
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ft    LET    US   TAKE    THE    SAME    DATA    AND    DO    A    TWO    WAY    ANCVA 


0PTI0N+-2 


ANOVA 
ENTER    NUMBER   OF    TREATMENTS, 
Dt 

4 
ENTER   IE.CZO&    OF    NUMBER    OF   OBS.    FOR   TREATMENTS    1  TO    4 
Di 

3  3  3  3 
ENTER    3  OBSERVATIONS   FOR   TREATMENT    1 
Di 

45  46  51 
ENTER    3  OBSERVATIONS    FOR   TREATMENT    2 
Di 

42  44  50 
ENTER    3  OBSERVATIONS    FOR   TREATMENT    3 
Di 

36  41  48 
ENTER    3  OBSERVATIONS    FOR   TREATMENT    4 
Di 

40  47  54 

ANOVA    TABLE 
SOURCE  DF  SS  MS  F 

TREATMENT  3         60.67       20.22     5.56 

BLOCKS  2  204.17      102.08    28.05 

ERROR  6  21.83        3.64 

TOTAL  11        286.67 

RmSQUARE    =    0.9  24 

OVERALL    MEAN    =  4  5.33 

TREATMENT    EFFECTS         2.00   0.00  "3.67   1.67 

BLOCK    EFFECTS      "4.5  8  "0.8  3   5.42 
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CONTINGENCYHOW 
SYNTAX  i       CONTINGENCY    MATRIX 
GROUP i  TWOWAY 

DESCRIPTION i       CONTINGENCY    TAKES    A    TWO-WAX    TABLE       AND       PERFORM? 
TEST    OF      INDEPENDENCE      OF      ROWS      AND      COLUkNS.  CUI-l 

STATISTIC    AND    ITS    DEGREES    OF    FREEDOM    ARE    OUTPUT. 


CONTINGENCY    DEATHS 
CHInSQUARE    VALUE-    2.920832627  UF-    12 


n  COMPARE    THE    ABOVE    VALUE    WITH    THOSE    IN    CHI-SQUARE 

n  TABLES    FOR    12  OF    TO    DETERMINE    IF    THE    HYPOTHESIS 

a  STATING    THAT    THE    ROWS    AND    COLUMNS    ARE    If\iDEPENDl   .. 

a  CAN    BE    REJECTED    AT    THE    SELECTED    SIGNIFICANCE    LEVEL 
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)GRP   STATS 
CORRELATION  CR  FMT  MSTATS      PARTIAL    STATI, 

CORRELATIONHOrf      MSTATSHOtf  PARTI  ALHOW  ST  ATI  FT ICS HO  i 


FTATISTICFROW 
SYNTAX t       STATISTICS    VECTOR 
GROUP i       STATS 
SUBPROGRAM:       FMT 

DESCRIPTION t       STATISTICS    DETERMINES    THE    MEAN, 
DEVIATION,    COEFFICIENT    OF    VARIATION ,  LOWER    AND 
MEDIAN,       TRIHEAN,       MlDMEAN,       MIDRANGE,         RANGE 
DEVIATION,    INTERQUARTILE    RANGE,    SKEWNESS,    AND 
DATA    IN    VECTOR. 


FTATIFTICF    VOLCANO 
MEAN i     70.2  46575  3  4 
VARIANCE:     1850,56277b 
FTD.     DEV  .:     43. 01816796 
COEFF.    OF    VARIATION:     0.6123881165 


VARIANCE, 

STA 

NDA 

F\  0 

UPPER      Q 

UART 

ILh 

C 

'    5 

,         Pi  E  A  N 

AbS 

OLU 

TE 

KURTOFIF 

FOR 

m 

HE 

LOWER    QUARTILEi     3  7 

UPPER    QUARTILEi    9  6 

MEDIAN:     6  5 

T RIME  AN:     6  5.75 

MIDMEANi     64.47747748 

RANGE:     197 

Ml D RANGE:     10  0.5 

MEAN    ABSOLUTE    DEVIAT 

ION  : 

33.68493 

interquartile  range: 

59 

COEFF.    OF    SHEWN  ESS t 

0.83257670  5  9 

COEFF.    OF    KU PTOSIS: 

0  .  4  2  8 

1831164 
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CORRELATIONHOW 
SYNTAX:       R+C CRRELATION    W 
GROUP:       STATS 
SUBPROGRAM:       FMT 

DESCRIPTION:  CORRELATION  DETERMINES  THE  SIMPLE  PEARSON  PRODUCT- 
MOMENT  CORRELATIONS  BETWEEN  EACH  PAIR  OF  VARIABLES  REPRESENTED  BY 
THE    C    COLUMNS    OF    W.       THE    OUTPUT    IS    A    C    BY    C    CORRELATION    MATRIX. 
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15 


CORRELATION    GRADESl 


1.00 

2.0  0 

3.0  0 

4.00 

5.00 

I. 00 

1.0  0 

0.34 

0.5  0 

0.2  5 

0.6  4 

2.0  0 

0.3  4 

1  .00 

"0.07 

0.5  4 

0.68 

3.00 

0.  SO 

"0.07 

1.0  0 

0.52 

0.64 

4.00 

0.25 

0.5  4 

0.5  2 

1.0  0 

0.8  7 

5.0  0 

o.cu 

0.6  8 

0.64 

0.87 

1.0  0 

n  IN    THE    CORRELATION    OUTPUT,    THE    TOE    ROM    A  fib    Tiih    FIRST 

a  COLUMN    ARE    LABELS    FOR    THE    CORRELATIONS.       FOR    EXAhP;      , 

fl  THE    CORRELATION    BETWEEN    THE    SECOND    AND    FOURTH    CQl      . 

fl  OF    GRADES1    IS    0.54. 


pGRADESl 


MSTATSdOW 

SYNTAX:       MS  TATS    W 
GROUP i       STATS 
SUB PROGRAM t       FMT 
DESCRIPTION i 


DETERMINES      CORRELATION, 


hEANS ,       :':.:A:i :.  Ahi 

DEVIATIONS,  LONER  QUARTILES,  MEDIANS,  AND  UPPER  QUARTILES  FOR  THE 
COLUMNS  OF  TnE  ARGUMENT  MATRIX  W.  THE  CORRELATIONS  ARE  UISPLAYEL 
IN    AN    li    BY    N    MATRIX,    WHERE    N    IS    THE    NUMBER    OF    COLUMNS    Ok    V. 


THE    OUTPUT    FROM    MSTATS    IS    MUCH    LIKE    THAT    OF    CCRRELATIi 


tiSTAT 

S    GRADE 

51 

1.00 

2.00 

3.00 

4.00 

5.00 

1.00 

1.00 

0.  34 

0.50 

0.2  5 

0.6  4 

2.0  0 

0.34 

1.0  0 

"0.07 

0.  54 

0.68 

3.0  0 

0.5  0 

"  0  .  0  7 

1.00 

0.52 

0.64 

4.0  0 

0.2  5 

0.5  4 

0.5  2 

1  .00 

0.8  7 

5.  00 

0.  64 

0.68 

0.6  4 

0.87 

1  .00 

MEANS 

0.86 

0.81 

0.88 

0.81 

0.8  4 

ST  D.DEI/ 

0.09 

0.2  1 

0.10 

0.10 

0.0  7 

L.    QRTLS 

0.82 

0.69 

0.82 

0.74 

0.78 

MEDIANS 

0.85 

0.83 

0.88 

0.79 

0.83 

U.    QRTLS 

0.94 

0.92 

0.98 

0.91 

0.90 
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PARTIALHCW 
SYNTAX i       PARTIAL    W 
GROUP:       STATS 

SUBPROGRAMS'.       CORRELATION,     FMT 

DESCRIPTION i  PARTIAL  DETERMINES  PARTIAL  CORRELATIONS  OF  Ai.Y 
SPECIFIED  ORDER  FOR  THE  COLUMNS  OF  THE  MATRIX  W.  TtiE  USER  JILL  BE 
ASKED  FOR  THE  COLUMNS  TO  BE  PART I AILED  OUT.  THE  USER'S  RESPONSES 
MUST  BE  INTEGERS  BETWEEN  1  AND  C,  WHERE  C  IS  THE  NUMBER  OF 
COLUMNS    OF    W. 


PART 

IAL    GRADES1 

ENTER    CCLU 

MNS    OF 

MATRIX    TO 

BE    PART 

IALLED 

Gl 

1 

1.00 

2.0  0 

3.00 

4.00 

1.0  0 

1.00 

"0.29 

0.50 

0.  64 

2.00 

"0.29 

1.00 

0.47 

0.49 

3.00 

0.5  0 

0.47 

1.0  0 

0.95 

4.00 

0.64 

0.49 

0.9  5 

1.00 

ow. 


fl  BE    CAREFUL    WITH    INTERPRETING    THE    LABELS   IN    THE    OUTPUT 

A  OF    PARTIAL.       ABOVE,    LABEL    1  DOES    NOT    MEAN    COLUMN    ONE, 

fl  BUT    THE    FIRST    COLUMN    REMAINING    AFTER   THE    SELECTED    VARIABLE, 

a  HAVE    BEEN    PARTI ALED    OUT.       SIMILARLY,    2    REFERS    TO    TtiE 

fl  SECOND   REMAINING    COLUMN. 


PARTIALLED    OW 


PART 

'IAL    GRADE 

SI 

ENTER   COLUMNS   OF 

MA 

TRIX    TO 

BE    PA 

Qi 

2  3 

1.00 

2.00 

3.00 

1.00 

1.00 

"0.45 

0.23 

2.00 

0.45 

1.00 

0.77 

3.00 

0.23 

0.77 

1.00 

A  THE    ABOVE    EXAMPLE    SHOWS    THAT    MORE    THAN    ONE    VARIABLE 

R  CAN    BE   PARTI ALED    OUT.       LABEL    2    ABOVE    REFERS    TO    THE 

fl  SECOND    COLUMN    REMAINING    AFTER    THE    SELECTED    VARIABLES 

R  HAVE    BEEN    PARTIALED    CUT.    IN    THIS    EXAMPLE    2 

fl  COLUMN    4. 
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)GRP    GFIT 
KS  KSHOW         CHISQUARE  C  HI  SQL' AkEHOW 


CRISQUAREHCW 
SYNTAX*      T    CHISQUARE   X 
GROUP i       GFIT 

DESCRIPTION i  CHISQUARE         COMPARES         A         THEORETICAL         DISCRhTf 

PROBABILITY    MASS    FUNCTION    T    HAVING      N      POSSIBLE      VALUES      WITH       An 
EMPIRICAL    FREQUENCY    FUNCTION    X    OVER    THE    SAME    N    VALUES    TO    TEST 
WELL    THE      PROBABILITIES      T      FIT      THE       DATA       X.  THE      CHI^SQUARP. 

GOODHESSmOFvFIT    STATISTIC    IS      PRINTED      OUT      WITH       THE      NUMBER 
DEGREES    OF    FREEDOM. 


n  THE    THEORETICAL    PROBABILITY    MASS    FUNCTION    I    A»D    THE 

fi  EMPIRICAL    FREQUENCY    FUNCTION    X    MUST    BE    OF    THE    SAME 

fi  SIZE;    I.E.,     pT=pX.       THE    ACTUAL    OBSERVED    VALUES    ARE 

fi  IMPLICIT    IN    THE    POSITION    OF    THE    VECTORS. 


T+.l     .2     .3     .2     .1 
X+    5  20  38  32  12 


fi  THE    VALUES    IN    T    REPRESENT    PROBABILITIES.       THE    VALUES 

p  X    REPRESENT    OBSERVED    FREQUENCIES.       THE    ACTUAL    OBSERVAT1     ■ 

n  ARE    NOT    STATED    EXPLICITLY,    BUT    IT    IS    UNDERSTOOD    ?H* 

fl  THE    POSITIONS    OF    X    AND    T    REFER    TO    THE    SAME    VALUES.       ECR 

fl  EXAMPLE,    IN    THE    CASE    ABOVE,    THE    FIRST    POSITION    MAY 

n  TO    THE    VALUE       1,     THE    SECOND    0,  ETC. 


T    CHISQUARE    X 
CHI-SQUARE     VALUE    =     9.620872274  DF    =    4 


n 

n  THE    DATA    WERE    DRAWN    FROM    A    POPULATION    WITH    PRCbt 

p  MASS    FUNCTION    T,    COMPARE    THE    VALUE    ABOVE    WITH    THOSE 

p  IN    CHI  ^-SQUARE    TABLES    WITH    4  DF. 
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KSUOV 

SYNTAX i       F    KS    G 

GROUP i       GFIT 

DESCRIPTION i       KS    DORS    A    KOLMOGOROV vSMIRNOV    GOODNESS    OF      FIT      TEST 

OF    THE    EQUALITY         OF         TWO         DISCRETE  EMPIRICAL  PROBABILITY 

DISTRIBUTIONS,       THE      ARGUMENTS      F      AND      G      ARE      EACH       VECTORS      OF 

OBSERVATIONS.       THEY    NEED    MOT    BE    SORTED,    NOR    MUST    THEY    BE    THE    SAME 

SIZE.       OUTPUT    OF    KS    ARE,    THE    VALUES  \ 

MAX    F{X)**G(.X)    AND    MIN    FU)-G(A') 


F+3    2    5    9    10    16    3    5    7    19    20    12    6    3    8    8    9    7    7    7    7    2    10    7 
G+U    21020133674211086    3    317 


F    KS    G 
MAX    F(X)»G(X)    =0.12745 09804 

MIN    F{X)r-,G{X)     =  0.3137254902 


A  TO    DETERMINE    IF    THE   TWO    SAMPLES    OF    DATA    WERE    TAKEN    FROM 

a  THE    SAME    POPULATION     (OR    POPULATIONS    WITH    THE    SAME    CDF'S) 

A  COMPARE    THE    MAX    AND    MIN    VALUES    WITH    THE    DELTA    STATISTICS 

a  IN    TABLES    FOR    THE    TWO    SAMPLE    KS    TEST. 
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APPENDIX    A 
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AN  DHOW 
SYNTAX:   X   AND   Y 

GROUPS:   RELATIONS,  TWOWAY,  COMPARISONS 

DESCRIPTION:  THE  FUNCTION  'AND'  IS  USED  TO  CREATE  A  NEW  DATA 
ARRAY  CONSISTING  OF  X  AUGMENTED  BY  Y  AS  ADDITIONAL  COLUMNS.  X 
AND  Y  CAN  BE  SCALARS,  VECTORS,  OR  MATRICES  BUT  CANNOT  BOTH  BE 
SCALARS.  IF  AN  ARGUMENT  IS  A  SCALAR,  A  COLUMN  IS  GENERATED  EACH 
ELEMENT  OF  WHICH  IS  THE  SCALAR.  IF  THE  ARGUMENTS  ARE  VECTORS  OR 
MATRICES,  THEY  MUST  BE  CONFORMABLE. 


ANOVAHOW 
SYNTAX:  ANOVA 
PARAMETER: 

OPTION'  SELECTS  THE  TY 
GIVES  A  ONE-WAY  ANOVA 
ANOVA  (TREATMENTS  WITH 
GROUP:   TWOWAY 

SUBPROGRAMS:  INPUT,  SUMSQ,  FMT 
DESCRIPTION:  ANOVA  DOES  A  ONE- 
VARIANCE  DEPENDING  ON  THE  VALUE 
WILL  INTERACT  WITH  THE  USER  TO  0 
AND  DATA.  WHEN  ENTERING  DATA,  S 
LEAST  ONE  BLANK  AND  USE  NO  OTHER 
AN  ANOVA  TABLE,  ESTIMATES  OF  THE 
THE  TREATMENT  (AND, IF  OPTION-2, 


PE  OF  ANOVA  PERFORMED.  OPTION =1 
AND  0PTI0N=2  GIVES  A  TWO-WAY 
BLOCKING).   (DEFAULTED 


WAY  OR  A  TWO-WAY  ANALYSIS  OF 
OF  THE  OPTION  PARAMETER.  ANOVA 
BTAIN  THE  REQUIRED  INFORMATION 
EPARATE  THE  DATA  POINTS  WITH  AT 
DELIMITERS.  OUTPUT  CONSISTS  OF 
OVERALL  MEAN,  AND  ESTIMATES  OF 
BLOCK)  EFFECTS. 


A3RH0W 
SYNTAX:   A3R  VECTOR 
GROUP:   SMOOTH 
SUBPROGRAM:   0NE3 

DESCRIPTION:  A3R  DOES  REPEATED  SMOOTHINGS  OF 
UNTIL  THERE  ARE  NO  CHANGES  IN  THE  SMOOTHED 
ITERATION  TO  THE  NEXT.   SEE  0NE3H0W . 


RUNNING  MEDIANS  OF  3 
SEQUENCE   FROM   ONE 


A3RSRH0W 
SYNTAX:   A3RSR  VECTOR 
GROUP:   SMOOTH 

SUBPROGRAMS:   A3R,  0NE3,  SPLIT 

DESCRIPTION:   A3RSR  DOES  REPEATED  SMOOTHINGS  BY   RUNNING   MEDIANS 
OF  3  FOLLOWED  BY  SPLITTING  MESAS  AND  REPEATING  UNTIL  CONVERGENCE. 
SEE  A3R  AND  SPLIT. 
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BOXPLOTHOW 
SYNTAX:       BOXPLOT    VECTOR 
PARAMETER: 

LENGTH'    CONTROLS    THE    HORIZONTAL      SIZE       OF      THE      DISPLAY 
(DEFAULT=bO      CHARACTERS) . 
GROUP:       DISPLAY 
SUBPROGRAM:       FILL 

DESCRIPTION :  BOXPLOT  GENERATES  A  BOX  PLOT  DISPLAY  FOR  A  VECTOR 
OF  DATA.  A  RECTANGULAR  BOX  WITH  ENDS  CORRESPONDING  TO  LOWER  AND 
UPPER  QUARTILES  IS  PRESENTED  WITH  THE  MEDIAN  MARKED  WITH  AN 
ASTERISK.  WHISKERS  ARE  DRAWN  ON  EACH  SIDE  OF  THE  BOX  WITH 
CROSSES  MARKING  THE  LOWEST  AND  HIGHEST  DATA  VALUES  WITHIN  AN 
INTERQUARTILE  DISTANCE  OF  THE  QUARTILES.  DATA  VALUES  OUTSIDE  THE 
CROSSES  (OUTLIERS)  ARE  MARKED  WITH  CIRCLES  AND  THOSE  MORE  THAN 
1.5  INTERQUARTILE    DISTANCES    GET    HEAVY    CIRCLES. 


CHISQUAREHOW 
SYNTAX:       T    CHI  SQUARE    X 
GROUP:       GFIT 

DESCRIPTION:  CHISQUARE         COMPARES         A         THEORETICAL         DISCRETE 

PROBABILITY  MASS  FUNCTION  T  HAVING  N  POSSIBLE  VALUES  WITH  A>v 
EMPIRICAL  FREQUENCY  FUNCTION  X  OVER  THE  SAME  N  VALUES  TO  TEST  HOW 
WELL    THE      PROBABILITIES      T      FIT      THE      DATA      X.  THE      CHI-SQUARE 

GOODNESS~OF-FIT  STATISTIC  IS  PRINTED  OUT  WITH  THE  NUMBER  OF 
DEGREES    OF    FREEDOM. 


CODERESHOW 
SYNTAX:       CODE RES    MATRIX  OR  CODES  RES    MEDPOLISH    MATRIX 

SUBPROGRAM:       AND 
GROUP:       TWOWAY 

DESCRIPTION:       CODERES    PRODUCES    A    DIAGNOSTIC    ARRAY       FOR       AHALYZ1 
THE    RESIDUALS    FROM    MEDIAN    POLISHING    OF    A    TWO -WAY    TABLE.       THE    USER 
SHOULD    SET    RESIDS-1    IN    MEDPOLISH      WHEN      CODERES      IS      USED.  THlL 

OUTPUT    MATRIX    CONSISTS    OF    THE    SYMBOLS    -°  +  DEPENDING    ON      THE      SIZE 
OF    THE    RESIDUALS. 

VALUE    BELOW    LOWER    QUARTILE 

VALUE    BETWEEN    QUARTILES 

VALUE    LARGER    THAN    UPPER    QUARTILES 
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COMPAREHOW 

SYNTAX:       COMPARE    MATRIX 
PARAMETERS: 

(1)  DEPTH-       VERTICAL       HEIGHT      OF      DISPLAY       (DEFAULT=20 
LINES) . 

(2)  MISS'      NUMBER      USED      TO      INDICATE      MISSING       VALUES 
(DEFAULT    ="99999) 

(3)  NGAP-      NUMBER      OF       HORIZONTAL       SPACES       BETWEEN       THE 
BOXPLOT    DISPLAYS    (DEFAULT=3) . 

GROUP:       COMPARISONS 
SUBPROGRAM:    FILL 

DESCRIPTION:  COMPARE  OPERATES  ON  AN  N  BY  K  MATRIX  TO  PRODUCE  K 
VERTICAL  BOX  PLOTS  PLACED  NEXT  TO  EACH  OTHER  TO  ALLOW  VISUAL 
COMPARISON  OF  THE  CENTERS,  SPREADS.  AND  OUTLIERS  OF  THE  BATCHES 
(COLUMNS)  OF  THE  MATRIX.  THE  USER  MUST  FILL  UP  THE  MATRIX  SO 
THAT    THERE    ARE    N    OBSERVATIONS    FOR    EACH       BATCH.  THE      PARAMETER* 

MISS,  SHOULD  BE  USED  TO  FILL  IN  MISSING  VALUES.  THE  USER  SHOULD 
ASSURE  THAT  MISS  IS  DIFFERENT  FROM  ALL  VALID  ENTRIES  IN  THE  DATA 
ARRAY . 


CONDENSEHOW 
SYNTAX:       CONDENSE   MATRIX    OR    R+CONDENSE    MATRIX 
PARAMETERS: 

(1)  NUM-    CONTROLS    WHAT    STATISTICS    ARE    INCLUDED    IN      THE 
SUMMARY.       NUMeil,2,5,7]    (DEFAULT=2) . 

NUM=1    GIVES    MEDIANS    OF    EACH    COLUMN    OF    THE    MATRIX 
NUM=2    GIVES    MEDIANS    AND   INTERQUARTILE    RANGES 
NUM-b    GIVES    MIN,QUARTILES,MAX,    AND    SIZE 
NUM=7    GIVES    MIN, EIGHTHS, MAX    AND    SIZE 

(2)  MISS'      NUMBER      USED      TO      CODE         MISSING         VALUES 
(DEF AULT-" 99999) . 

GROUP:      COMPARISONS 
SUBPROGRAMS:       UTCOND,    FMT 

DESCRIPTION:  CONDENSE  GENERATES  SUMMARY  STATISTICS  FOR  EACH 
GROUP  OF  DATA  REPRESENTED  BY  THE  COLUMNS  OF  THE  ARGUMENT  MATRIX. 
SINCE  THE  ARGUMENT  IS  A  MATRIX,  THE  USER  MUST  FILL  UP  THE  MISSING 
VALUES  IN  THE  MATRIX  WHENEVER  THE  GROUPS  (COLUMNS)  HAVE  DIFFERENT 
NUMBERS  OF  OBSERVATIONS.  THE  PARAMETER,  MISS,  SHOULD  BE  USED  TO 
FILL  THE  MATRIX  SINCE  THE  FUNCTION  WILL  RECOGNIZE  THOSE  VALUES  AS 
MISSING    DATA    AND   WILL    IGNORE    THEM   IN    ALL    CALCULATIONS.  THE 

RESULT    FROM    CONDENSE      CAN      BE      USED      AS      THE      ARGUMENT      OF      OTHER 
FUNCTIONS    SUCH    AS    SCAT    OR    LINE    IF    AND    ONLY    IF    NUM=2. 
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CONTINGENCYHOW 
SYNTAX:       CONTINGENCY    MATRIX 
GROUP:       TWOWAY 

DESCRIPTION:       CONTINGENCY    TAKES    A    TWO -WAY    TABLE 
TEST    OF      INDEPENDENCE       OF       ROWS      AND      COLUMNS. 
STATISTIC    AND    ITS    DEGREES    OF    FREEDOM    ARE    OUTPUT. 


AND       PERFORMS       A 
THE      CHI -SQUARE 


KSHOW 
SYNTAX:       F    KS    G 
GROUP:      GFIT 

DESCRIPTION :  KS  DOES  A  KO LMOGOROV -SMI RNOV  GOODNESS  OF  FIT  TEST 
OF    THE    EQUALITY         OF         TWO         DISCRETE  EMPIRICAL  PROBABILITY 

DISTRIBUTIONS.  THE  ARGUMENTS  F  AND  G  ARE  EACH  VECTORS  OF 
OBSERVATIONS.  THEY  NEED  NOT  BE  SORTED,  NOR  MUST  THEY  BE  THE  SAME 
SIZE.       OUTPUT    OF    KS    ARE    THE    VALUES: 

MAX    F(X)~G(X)    AND    MI N    FU)-G(X) 


INLINE    ARRAY 


LINEHOW 
SYNTAX:       LINE    ARRAY  OR 

PARAMETER: 

RSELECT-    SELECTS      THE       OUTPUT      MA 
RSELECT=1    GIVES    ABSCISSA    VS.       RES 
RSELECT*!    GIVES    FITTED    VALUES    VS. 
GROUP:       RELATIONS 

DESCRIPTION:       LINE    FITS    A    STRAIGHT    LINE    TO 
BY    DIVIDING    THE    POINTS    INTO    3  REGIONS    AND    U 
THE    X    AND    Y    VALUES    IN    THE    OUTER    REGIONS    TO 
THE    INTERCEPT    IS    THE    MEDIAN    OF    THE    DIFFEREN 
ALSO    COMPUTES    THE    RESIDUALS    AND    GIVES    AN    N 
EITHER    THE    X-VALUES      VS.       RESIDUALS      OR      TH 
RESIDUALS    DEPENDING    ON    THE    PARAMETER,     RSELE 
NOT    WANT    RESIDUALS    TYPED    OUT,    HE    SHOULD    USE 

INLINE  ARRAY 
THE  ARGUMENT  ARRAY  CAN  BE  EITHER  A  VECTOR  0 
IF  THE  ARGUMENT  IS  A  VECTOR  OF  SIZE  N 
CONSTRUCTED  TO  BE  THE  FIRST  N  POSITIVE  INTE 
IS  A  2-COLUMN  MATRIX,  THE  FIRST  COLUMN  IS  T 
X-VALUES  AND  THE  SECOND  COLUMN  THE  SET  OF  Y 
USED  IN  CONJUNCTION  WITH  THE  FUNCTION  SCAT 
THE    INTERCEPT,     AND    A    PLOT    OF    THE    RESIDUALS. 

SCAT    LINE    ARRAY 


TRIX       OF       RESIDUALS. 
IDUALS       (DEFAULTED  . 
RESIDUALS 

A    SET    OF    (X,Y)       POINTS 
SING      THE      MEDIANS      OF 
DETERMINE    THE    SLOPE. 
CES    Y-SLOPE*X.  LI  ■ 

BY    2    MATRIX    CONTAINING 
E      FITTED       VALUES       VS. 
CT.       IF    THE    USER       DOES 
THE    SYNTAX: 

R    A    2 -COLUMN    MATRIX. 

THE  X'VARIABLE  IS 
GERS.  IF  THE  ARGUMt  ■ 
AKEN  TO  BE  THE  SET  OF 
-VALUES.  LINE  CAN  BE 
TO  PRODUCE  THE  SLOPE, 
TO    DO    THIS,    ENTER: 
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LSLINEHOW 
SYNTAX:       LSLINE    ARRAY  OR  Z+LSLINE    ARRAY 

GROUP:       RELATIONS 

DESCRIPTION:  LSLINE      DETERMINES      THE       LEAST      SQUARES      SOLUTION, 

Y=YI+SLOPE*X.  THE  ARGUMENT  ARRAY  CAN  BE  A  VECTOR  OR  A  2-COLUMN 
MATRIX.  IF  THE  ARGUMENT  IS  A  VECTOR,  THE  X-VALUES  ARE  TAKEN  TO 
BE    THE    INTEGERS    1  TO    N   W 'HERE    N-Q ARRAY .  IF      THE      ARGUMENT      IS       A 

2-COLUMN  MATRIX,  COLUMN  1  CONSISTS  OF  THE  X-VALUES  AND  COLUMN  2 
THE  Y 'VALUES.  RESIDUALS  ARE  DETERMINED  AND  STORED  IN  A  MATRIX 
WHOSE  FIRST  COLUMN  IS  THE  SET  OF  X-VALUES  AND  WHOSE  SECOND  COLUMN 
IS  THE  SET  OF  RESIDUALS.  IF  THE  USER  WANTS  TO  SUPPRESS  PRINTOUT 
OF  THE  RESIDUALS,  HE  MUST  ASSIGN  THE  RESULTS  OF  LSLINE  TO  A 
VARIABLE;    I.E.,    TYPE    Z+LSLINE    ARRAY. 


MEDPOLISHHOW 
SYNTAX:       MEDPOLISH   MATRIX  OR  Z+MEDPOLISH    MATRIX 

PARAMETERS: 

(1)  RESIDS-  CONTROLS  OUTPUT  OF  RESIDUALS.  IF  1, 
RESIDUALS  ARE  PRODUCED  IN  A  TWO-WAY  TABLE;  IF  2, 
COMPARISON  VALUES,  ( RE<>  .  *CE)  iTV  AND  RESIDUALS  ARE 
PRODUCED   IN    A    TWO    COLUMN    MATRIX    (DEFAULTED . 

(2)  EPSILON-  PROPORTION  BY  WHICH  THE  SUM  OF  THE 
ABSOLUTE  VALUES  OF  THE  RESIDUALS  MUST  BE  REDUCED  AT 
EACH    ITERATION    TO    CONTINUE    POLISHING    (DEFAU LT=0 .01) . 

(3)  NORMEFFECTS-  DETERMINES  WHETHER  OR  NOT  NORMALIZED 
EFFECTS  ARE  OUTPUT.  THE  DEFAULT  OF  0  SUPPRESSES  OUTPUT 
OF  NORMALIZED  VALUES.  NORMEFFECTS-1  CAUSES  NORMALIZED 
VALUES    TO    BE    PRINTED. 

GROUP:      TWOWAY 

DESCRIPTION:       MEDPOLISH    ITERATIVELY    SWEEPS    OUT    MEDIANS       FROM      THE 

ROWS    AND   COLUMNS    OF    A    TWO-WAY    TABLE    TO    YIELD    THE    MODEL: 

OBS    =    MEDIAN    +  ROW    EFFECT    +  COLUMN    EFFECT    +  RESIDUAL 
IT    YIELDS    RESIDUALS    FOR    TESTING    THE    ADEQUACY    OF      THE      MODEL       (SEE 
SHOWRES    AND   CODERES) .    THE    FUNCTION    CONTINUES    POLISHING   UNTIL      THE 
STOPPING    RULE   CONTROLLED    BY    EPSILON    IS    ACTIVATED.  THE      OVERALL 

MEDIAN,    THE    ROW-    AND    COLUMN -EFFECTS ,       AND      THE      SUM      OF      ABSOLUTE 
VALUES    OF    THE      RESIDUALS      ARE      GIVEN.  THE      USER      CAN      SUPPRESS 

RESIDUAL    PRINTOUT    BY    ASSIGNING    THE    FUNCTION    TO    A    VARIABLE,       I.E., 
BY    WRITING: 

Z<-  MEDPOLISH    MATRIX 
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CORRELATIQNHOtf 
SYNTAX  I  R+CORRELATION    W 
GROUP i       STATS 
SUBPROGRAM \       FMT 

DESCRIPTION t  CORRELATION  DETERMINES  THE  SIMPLE  PEARSON  PRODUCT 
MOMENT  CORRELATIONS  BETWEEN  EACH  PAIR  OF  VARIABLES  REPRESENTED  BY 
THE    C    COLUMNS    OF    W.       THE    OUTPUT    IS    A    C    BY    C    CORRELATION    MATRIX. 


Id  ST  AT  SHOW 
SYNTAX i  MSTATS    W 
GROUP i       STATS 
SUbPRCGRAMx       FMT 

DESCRIPTION i       MSTATS      DETERMINES      CORRELATIONS,       MEANS, 
DEVIATIONS,    LOWER    QUARTILES,    MEDIANS,    AND    UPPER    QUARTILES    Ft 
COLUMNS    OF    THE    ARGUMENT    MATRIX    W.       THE    CORRELATIONS    ARE    DISPLAYED 
IN    AN    N    BY    N    MATRIX,    WHERE    N    IS    THE    NUMBER    OF   COLUMNS    OF   V. 


PARTIALHCW 
SYNTAX i       PARTIAL    W 
GROUP:      STATS 

SUBPROGRAMS:       CORRELATION ,    FMT 

DESCRIPTION:  PARTIAL  DETERMINES  PARTIAL  CORRELATIONS  OF  ANY 
SPECIFIED  ORDER  FOR  THE  COLUMNS  OF  THE  MATRIX  W.  THE  USER  WILL  BE 
ASKED  FOR  THE  COLUMNS  TO  BE  PARTI ALLED  OUT.  THE  USER'S  RESPONSES 
MUST  BE  INTEGERS  BETWEEN  1  AND  C,  WHERE  C  IS  THE  NUMBER  OF 
COLUMNS    OF    W. 
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NUMSUMHOW 

SYNTAX:       NUMSUM    VECTOR  OR  NUMSUM    ARRAY 

PARAMETER: 

MISS-    NUMBER    USED    TO    INDICATE    MISSING    VALUE    IN    THE    DATA 

ARRAY    (DEFAULT="99999) . 
GROUP:       COMPARISONS 
SUBPROGRAM:       FMT 

DESCRIPTION:  NUMSUM  OPERATES  ON  EITHER  A  VECTOR  OR  AN  ARRAY  BUT, 
IN    ALL    CASES,    THE    DATA    ARE    TREATED    AS    A       SINGLE       BATCH.  NUMSUM 

PRODUCES  A  NUMERICAL  SUMMARY  WHICH  GIVES  THE  SAMPLE  SIZE  {AFT Eh 
DELETION  OF  MISSING  VALUES),  THE  EIGHTHS,  THE  EXTREMES,  THE 
SPREADS,     AND    THE      MIDPOINTS    IN    TABULAR    FORM.  THIS      SUMMARY       IS 

USEFUL  FOR  TESTING  THE  SYMMETRY  OF  A  DATA  SET,  AND  TO  EVALUATE 
THE    EFFECTIVENESS    OF    DATA    TRANSFORMATIONS    IN    PRODUCING    SYMMETRY. 


ONEHHOW 
SYNTAX:       ONEH    VECTOR 
GROUP:       SMOOTH 

DESCRIPTION:       ONEH    DOES    A    SINGLE    HANNING    OF    A    SEQUENCE    OF    DATA. 
IT   IS   USED    IN    CONJUNCTION   WITH    THE    OTHER    SMOOTHING    FUNCTIONS    AS    A 
FINAL    TOUCH-UP    TO    A    SEQUENCE    OF    SMOOTHED    DATA.       THE    I-TH    RESPONSE 
IS    REPLACED    BY    0.25*y[r-l]  +  0.5*Y[I]  +  0 . 2 5*7 [1+1] .  SEE    A3R    AND 


A3RSR. 


0NE3H0W 
SYNTAX:       0NE3    VECTOR 
GROUP:       SMOOTH 

DESCRIPTION:  0NE3  SMOOTHS  A  VECTOR  OF  DATA  USING  ONE  PASS  OF 
RUNNING  MEDIANS  OF  3  WITH  TU KEY'S  END-POINT  RULE.  WHEN  SMOOTHING 
(X,Y)  PAIRS,  THE  USER  SHOULD  USE  AS  THE  ARGUMENT  VECTOR  THE  Y 
VALUES    ORDERED    ACCORDING    TO    THE    MAGNITUDE    OF    THE    X    VALUES: 

VECTOR+YlbX] 
IT    IS    ASSUMED      THAT      THE      X       VALUES       ARE      EQUISPACED.  REPEATED 

SMOOTHING    BY    MEDIANS    OF    3    CAN    BE    ACCOMPLISHED   WITH    A3R.  OUTPUT 

CONSISTS    OF    THE    SMOOTHED      SEQUENCE.  A       PLOT      OF      THE      SMOOTHED 

SEQUENCE    IS    GIVEN    BY    TYPING: 

SCAT    0NE3    VECTOR 
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REGRESSHCW 
SYNTAX i       Z+Y    REGRESS   X 
PARAMETER t 

^INTERCEPT*    DETERMINES    WHETHER    OR    NOT    AN    INTERCEPT    TERM 
IS   TO    BE    INCLUDED.         MNTERCEPT=1      GIVES      AN      INTERCEPT 
TERM,    AND    MNTERCEPT=0    GIVES    NO    INTERCEPT.     (DEFAULT      IS 
1.) 
GROUP t       RELATIONS 
SUBPROGRAMS:       FMT    AND   SCAT 

DESCRIPTION:  REGRESS      DOES      A      MULTIPLE      REGRESSION         ANALYSIS 

RELATING  THE  DEPENDENT  VARIABLE  Y  TO  A  SET  OF  CARRIERS  X.  THE 
LEFT  ARGUMENT  Y  IS  A  VECTOR  OF  SIZE  N.  THE  RIGHT  ARGUMENT  X  IS 
AN  N  BY  K  MATRIX  CONSISTING  OF  N  OBSERVATIONS  ON  EACH  OF  K 
VARIABLES  OR  A  VECTOR  OF  SIZE  N  IF  K=l .  OUTPUT  CONSISTS  OF  AN 
ANCVA  TABLE,  RmSQUARE,  STD.  ERROR,  REGRESSION  COEFFICIENTS  (I, it 
FIRST  COEFFICIENT  IS  THE  CONSTANT  TERM  IF  AI NTERCEPT= 1 . ) , 
TnSTATISTICS,  V ARI ANCEnCCV  ARI ANCE  MATRIX,  DURBIN ~W ATEON 

STATISTIC,  AND  A  VECTOR  OF  PREDICTED  Y  VALUES  AND  RESIDUALS. 
THERE  IS  AN  OPTION  THAT  ALLOWS  THE  USER  TO  INPUT  A  VECTOR  OF  X 
VALUES  AND  USE  THE  REGRESSION  EQUATION  TO  FORECAST  Y  VALUES.  THE 
USER    CAN    ALSO    OBTAIN    A    SCATTER      PLOT      OF      THE      RESIDUALS.  Wh 

EXECUTION  TERMINATES ,  THE  PREDICTED  Y  VALUES  AND  THE  RESIDUALS 
RESIDE    IN    THE    N    BY    2  MATRIX    Z. 


STATISTICSHOW 
SYNTAX i       STATISTICS    VECTOR 
GROUP i       STATS 
SUBPROGRAM:       FMT 

DESCRIPTION:  STATISTICS  DETERMINES  THE  MEAN,  VARIANCE,  STANDARD 
DEVIATION,  COEFFICIENT  OF  VARIATION,  LOWER  AND  UPPER  QUARTILES, 
MEDIAN,       T RIME  AN,       MIDMEAN,       MID  RANGE,  RANGE,  MEAN         ABSOLUTE 

DEVIATION,  INTERQUARTILE  RANGE,  SKEWNESS,  AND  KURTOSIS  FOR  THE 
DATA    IN    VECTOR. 
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SCATHOU 

SYNTAX i       SCAT    ARRAY 
PARAMETERS x 

(1)  HID*    CONTROLS    THE    HORIZONTAL    SIZE       OF      THE       DISPLAY 
(DEFAULT=30    CHARACTERS) , 

(2)  DEP*    CONTROLS    THE      VERTICAL      SIZE      OF      THE      DISPLAY 
(DEFAULT=lb    LINES), 

(3)  NDIVX,    NDIVYr    NUMBER    OF    UNITS      ON       X~       AND       Y-AXES, 
RESPECTIVELY    ( DEFAULT  =  4,4)  . 

GROUP x       RELATIONS,    SMOOTH ,    COMPARISONS 

DESCRIPTIONx  SCAT  PRODUCES  A  SCATTER  PLOT  OF  THE  DATA  CONTAINED 
IN  ARRAY,  THE  ARGUMENT  ARRAY  CAN  BE  A  VECTOR  OR  A  MATRIX  WITH  AS 
MANY  AS  9  COLUMNS,  IF  A  VECTOR  OF  SIZE  N,  THOSE  VALUES  ARE 
PLOTTED  VS,  THE  INTEGERS  1  TO  N;  IF  A  MATRIX,  THE  SECOND,  THIRD, 
ETC,  COLUMNS  ARE  PLOTTED  VS,  COLUMN  1  ON  THE  SAME  AXES, 
DIVISIONS  ON  THE  AXES  OF  THE  PLOT  ARE  NOT  EXPLICITLY  PRINTED, 
EXCEPT  AT  THE  EXTREMES  OF  THE  PLOT,  THE  USER  CAN  CONTROL  THE 
RESOLUTION  OF  THE  PLOT  BY  MODIFYING  THE  PARAMETERS  WID,  DEP, 
NDIVX    AND    NDIVY ,  PRINTING      TIME      INCREASES      DRAMATICALLY       WITH 

RESOLUTION;  THEREFORE,  OTHER  PLOT  PROGRAMS  SHOULD  BE  USED  IF  HIGH 
RESOLUTION  IS  DESIRED,  FOR  A  SINGLE  GROUP  OF  DATA,  THE  NUMBER  N 
(2*0*9)  WILL  BE  PRINTED  IF  N  POINTS  LIE  CLOSE  TOGETHER  ON  THE 
DISPLAY,  FOR  MULTIPLE  PLOTS  ON  THE  SAME  DISPLAY,  THE  LETTER  A 
REPRESENTS    GROUP    1,    B    REPRESENTS    GROUP    2,       ETC,  FOR      2       POINTS 

CLOSE    TOGETHER   THE    LETTER    WILL    BE    PRINTED    WITH      AN      UNDERSCORE, 
THE    DISPLAY    CANNOT    HANDLE    3  OR      MORE      POINTS      CLOSE      TOGETHER      IF 
THERE    ARE   MULTIPLE    PLOTS, 
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SHOWRESHOW 
SYNTAX:       SHOW RES    MATRIX  OR  SHOW  RES    MEDPOLISH    MATRIX 

SUBPROGRAM:       AND 
GROUP:       TWOWAY 

DESCRIPTION:  SHOWRES  AIDS  IN  THE  ANALYSIS  OF  RESIDUALS  FROM 
MEDIAN  POLISHING  A  TWO 'WAY  TABLE.  THE  USER  SHOULD  SET  RESIDS=1 
IN  MEDPOLISH  IF  SHOWRES  IS  TO  BE  USED.  THE  OUTPUT  OF  SHOWRES  IS 
A  MATRIX  CONSISTING  OF  THE  SYMBOLS  OOoxX  DEPENDING  ON  THE  SIZES 
OF  THE  RESIDUALS.  LET  R  REPRESENT  THE  RESIDUAL,  M  THE  MIDSPREAD 
OF  ALL  RESIDUALS,  L  THE  LOWER  QUARTILE,  AND  U  THE  UPPER  QUARTILE. 
THE    CODES    ARE    PRINTED    AS    FOLLOWS: 


0: 

R<L~1.5M 

o  : 

L-l.SMZRZL-M 

O  J 

L~M<R<U+M 

x  : 

U+M<R<U+1.SM 

X: 

R>U+1.5M 

SPLITHOW 
SYNTAX:       SPLIT    A3R    VECTOR 
GROUP:       SMOOTH 

DESCRIPTION:  SPLIT  DOES  ONE  PASS  AT  DIVIDING  MESAS  (PAIRS  OF 
ADJACENT  POINTS  WITH  A  COMMON  VALUE  WHICH  IS  A  LOCAL  MAX  OR  MIN) 
USING  TU KEY'S  END-POINT  RULE.  IT  IS  USED  IN  CONJUNCTION  WITH 
A3R. 
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STEMLEAFHOW 
SYNTAX:       STEMLEAF    VECTOR 
GROUP:       DISPLAY 
PARAMETERS: 

(1)  WIDTH'    CONTROLS    THE    WIDTH    (CHARACTERS    PER    LIVE)    Of 
THE    DISPLAY    ( DEFAU LT=7 0) . 

(2)  SCALE-    VARIES    THE    DEPTH       (STEM      INTERVAL)       OF      THE 
DISPLAY    IN    UNITS    OF    1,  2,  OR    0.5  TIMES    A    POWER    OF    10. 
SELECT    AN    INTEGER    FROM    1    TO    3    (DEFAULTED  . 

SUBPROGRAM:       LIT 

DESCRIPTION:  STEMLEAF  GENERATES  A  STEM  AND  LEAF  DISPLAY  OF  A 
VECTOR  OF  OBSERVATIONS.  THE  FUNCTION  AUTOMATICALLY  SCALES  THE 
DATA  USING  A  SCALING  ROUTINE  BASED  ON  THE  RANGE  AND  SIZE  OF  THE 
DATA  VECTOR.  BY  VARYING  THE  SCALE  PARAMETER  THE  USER  CAN  CHANGE 
THE  SIZE  OF  THE  STEM,  BUT  IT  WILL  NOT  NECESSARILY  BE  SCALED  BY 
THE  AMOUNT  SPECIFIED.  TRY  VARIOUS  VALUES  LIKE  1,2,  AND  3  TO  SEE 
WHICH  CHOICE  YIELDS  THE  BEST  DISPLAY.  IF  A  LEAF  CONTAINS  TOO 
MANY  CHARACTERS  TO  BE  PRINTED  ON  A  LINE,  THE  LINE  WILL  BE 
TRUNCATED  AND  THE  NUMBER  OF  TRUNCATED  CHARACTERS  WILL  BE  WRITTEN 
AT  THE  END  OF  THE  LINE.  IF  THE  STEM  INTERVAL  IS  TWO  TIMES  A 
POWER    OF    TEN,    THE    LEAVES    MAY    CONTAIN    THE      CHARACTERS       A,B,C,  .  .  .  ,J 

IN    ADDITION      TO       THE      DIGITS       0,1 9.     IN      THIS      CASE,  THE 

CHARACTERS    A    TO   J    REPRESENT    THE    NUMBERS    10  TO    19. 


TWICEHOW 
SYNTAX:       TWICE    VECTOR 
GROUP:       SMOOTH 

SUBPROGRAMS:       ONE3,    SPLIT,    A3RSR 

DESCRIPTION:  TWICE  SMOOTHS  A  SEQUENCE  OF  DATA,  THEN  SMOOTHS  THE 
RESIDUALS  AND  ADDS  THE  SMOOTHED  RESIDUALS  BACK  TO  THE  SMOOTHED 
DATA  TO  OBTAIN  THE  FINAL  SMOOTH.  THE  SMOOTHING  IS  DONE  USING 
A3RSR. 
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nut: FAULT    VALUES    OF    THE    PARAMETERS 


DEP 
1  5 

DEPTH 
20 

EPSTLON 
0.01 

LENGTH 
50 

MISS 
"99999 

NDIVX 


4 

4 

3 

0 

2 

1 

1 

1 

1 

30 

70 

1 


NDIVY 

?J(JAP 

iWRMEFFECTS 

MUM 

OPTION 

REST  OS 

SCALE 

RSELECT 

WID 

WIDTH 

^INTERCEPT 

VI 
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APPENDIX    B 


92 


V  L+A    AND    B;C;D 

[I]  +(  (  (2<pp4)  v3<pp£),0*pp£)/     17     3 
[2]  B<-,B 

[3]  ->(  (  (3=pp£)Ai*ipp£)  ,  2  =  pp/1)  /    17    7 

[4]  4*-,  4 

[5]  -►(A/((piO*l,0),l*tf«-lp~2<|>p£i)/16 

[6]  A-(({D*pA)LU[ pA)*l)pA 

[7]  +  (l*ppM/9 

[8]  B<-UipB)[  il=pB)*lppA)  9l)pB 

[9]  -*•(  (A/£l*l,lpp/5)  ,  l*D<-lp~2<t>ptf  )/     16     11 

[10]  fl*(((3  =  ppB)pl),(lpp/l),lp<l>pfl)pfl 

[II]  +(3=ppb)/14 

[12]  L^(  (  (f^lp<t>p>3)p0),(lp<t>p£)pl)\B 

[13]  ->QxppL[  ;  \C]^/4 

[14]  ZXl,  (  (C^lp<t>p/1 )  p0  )  ,  (~l  +  lp<t>p£)pl  )\B 

[15]  -*-0xppL[;  ;l  +  if  ]«-/! 

[16]  -►o=pi>,42?Gtf0£/vr7ls  of     mo     /3A&'  .'vcr  conformable 

[17]  '/3/V  ARGUMENT    OF      AND      IS    Of   IMPROPER    RANK..1 


[1] 
[2] 


AN  01/ A  ;Y  iK;  BTSS  ;  MSR  \R\MAX\  BBSS  ;  MSBL\MSE  ;  VTi'S  ;  fi'tff  ;  TDE  ;  FSS ;  y  b  A  :  ;  .  ; 

T\BLOCK 

INPUT 

SUMSQ 


V/33/f[Q]V 
V    Z-A3R    X\Y1\Y2\XI\X2\N\C 
[1]         Z+OMEd    X 
[2]  +1*(  ppp*«-Z)xO*  +  /  |  Z*X 

V 


V,43/?6'S[U]V 
V    Z«-/43tfSi?    X 
[1]  X«-/33£    X 

[2]  Z-AZR    SPLIT    X 

[3]  ->2x(  pppX«-Z)xQ*  +  / 


Z*X 
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VBOXPLOTLUlV 
Z+BOXPLOT    XiL;U;B 


11] 

X-i  ilE~20+U*L)*(  {U-[ / 9X)* 

[2] 

0.1><L0.5  +  10xL;  (OT  LENGTH-*  6 

[3] 

+9*\0=U~L 

[4] 

Z^LENGTH    FILL    X 

[5] 

b+  •  "in  i" 

[6] 

B«-£,'  *|*|*  *x   ooooooooo 

[7] 

B+B9  '  _l L 1 l_      23456789 

L8] 

->0xppz^-(  3*LENGTH)pB 

[9] 

'RANGE    OF    DISTRIBUTION    IS 

LENGTti*L<-L/  9X  )+9XxLENGTHt1 

)p  '     '  ;0.  lx|_0.5  +  10x£/ 


»[i+z] 

®®®®®®®®® • [ 1+2 ] 
23456789  '  [1+Z] 

ZERO1 


[1] 
[2] 
[3] 

[4] 


ICBISQUARElWtf 

T   CHIS^UARE   X;DF;C 

-►3x1  ((Pn=Pn 

->0xpu^'  LENGTH    ERROR:    X    AND    T    MUST    BE    THE    SAME    SIZE 

£«-+/(  U-T)*2)*T«-Tx  +  /;j: 

}CtiI*SQUARE    VALUE    =     '  ;C;f       ^    =     ' ;(p*)*l 


VCODERESlUlV 

V  Z+CODERES    XiIiNlQiLtfiUW 

[1]  /!KpY«-YU  *«-,*] 

[2]  .3*-0.5xY[L0. 25x3 +/Vxi3]+Y[fO. 25x1 +/i/xx3] 

[3]  L^*«[l]»Z?^[3]-wCl] 

[4]  [/*^[3]+Z> 

[5]  Y«-U>Ltf)  +  U>4[l])  +  U>#[3] )+X>UU 

[6]  Z^f -»o++» [l+y] 

[7]  Z<-(  (pZ[;l]),2xpZ[l;])p(,Z)    ANU,(pZ)p' 
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VCOMPARElUlV 

V  Z+CCMPARE    XiUiLiYiMiCiJ iP 

[I]  X+(HE~20+U«Q+L)  *(  (U+U+-[  /Y)«DEPTH*L+1/Y«-{MISS*,X)/  tX)+X*DEP?H-l 
[  2  ]  JM  *  12?" 2 0+£' tL  )  x  ( V  m DEPTH  * L  )  +MISS*DEPTH«  1 

[3]  Z«-(0£,Px//,6xr*-pX[l;  Dp1     ' 

[U]  J+0 

[5]  <W  +  1 

[6]  P+DEPTH    FILL    PliP+(P*M)  I P+Xl  ;J  1  1 

[7]  Z[;~2  +  6x<7>'     |  "x_T  |  'll-'-Pl 

[8.1  Z[  ;  'l  +  btJl*-'    *T  =  i=     |*      ooooooooo    «®«®®««®®  '  [  l+p  ] 

[9]  Z[;6xj:]«-'     |"i_t|  23456789       2  34  5  6  7  89  »  [  1  +P  1 

[10]  +b*\J<C 

[II]  Z-^eZ 


V  CON  DENSE [QlV 
V  R+CON DENSE    X ;I iN ;NR W iQ iZ 
[1]    tf«-\0 
[2]    UTCOND    X 
[3]    +  0xi(/l/£'M*2) 
[4]    P«-(/VP,2)pZ 


vcontingencyluIv 

v  contingency  f ;r \c ;t ;e ;c hi2 ; df 

[1]  R*-  +  tF 

[2]  T++/C++/F 

[3]  E+C°.*R*T 

[4]  CHI2++/,   (  (Fn£')*2)Tf 

[5]  0FW(pF)"i 

[6]  'CHI'SQUARE    VALUE-     I;r«I2;t   £F  =  »  ;0F 
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VCONDENSElUlV 

V  R+C ON DENSE    X  ;I  ;/V  ;NRiU  \Q  I  Z 
[1]    fl«-iO 
[2]    UTCCND    X 
[3]    +  0x\WM*2) 
[4]    R-(NR,2)pZ 


VC0NTINGENCYLU1V 

V  CONTINGENCY    F  ;R ;C  ;T  ;E  -tCHI2  ;  DF 

[1]  fl«-  +  ^F 

[2]  T++/C++/F 

[3]  E+Co.xR±T 

[4]  CRI2+-*/  %{{F»E)*2)*E 

[5]  DF+*/{pF)*l 

[6]  'CHImSQVARE    VALUE-     1;CHI2;1       DF=     ' ;DF 


VCORRELATIONlQlV 

V  R+CORRELATION   W ;Z ;C ;S ;Cg;MEANS ;V AR 

[1.1    C^(^Z)  +  .xZ-HV«(pJ/)p(Afff>l55^  +  ^)*l+p("2  +  l,l,pf/)p{/  +  pCfi',-t  ' 
[2]    F^'BF8.2'  FWr(0,\p5)  ,[l](ip5),f  *5o.x5-h(7^/?5^  +  ^Z*2)* 

0.5 
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[1] 

[2] 

[3] 

[4] 

[5] 

[6] 

[7] 

[8] 

[9] 

[10] 

[11] 

[12] 

[13] 

[14] 

[15] 

[16] 

[17] 

[18] 

[19] 

[20] 

[21] 

[22] 

[23] 

[24] 

[25] 

[26] 

[27] 

[28] 

[29] 

[30] 

[31] 

[32] 

[33] 

VFMTLD1V 
V    OL+-1    FMT    R\S\W\LiGiX\T\K\J\MiQ\P\D\1S\0\LiBiViCUih 

0«-g«-h-m«-p/?«-(  ir~2tp/?)pfl 

02>((1=1+M)+     1     0     xA?+-:W  +  2t/^l<pC//<-££,  '  ,  '  )pA«-*0123456'789.  ' 

+ffM  (2I2«-0=2V)  v/«-l>PlS<-,£ 
L0:-»-~\r/v(  *P+HxQ  =  pK*-pX<-1     '  )  av/(  '^  '  ,(9^'H]'  )eS 

-+(LQ  +  (\/+-0  =  pS+J±S)+\B=Ml2]+-l)  ,L.~(  \*B+0  +  .  =K)  ,P*~' A' eK+K  AJ-S  \  \  '  ) 

iS 

+E  +  *pS<-'Tb'XT    DELIMITER1 

+L3-3xx(pG^A:  =  A>(£e~l  +  A)  /tf)  Lf/«-pX«-(  pK<-U  itf  )  t/O  +  (  -  (  4»: )  \0)+X 
&:+(D+~l+G+KeA)/L3>-2x(pK)*W+l<0+'XA1eK  +  (~r;(i  •     ,  '  ) /# 

+L3x  i  (2j*  +  /G)*xM[2]«-10i  I  l-2\\(2*«-|  1-Gi  0)  +  # 

+L3«4>6>,  -(L<-,2?PJ,£/0 /xtf«-10i  |  1-Ai(  |  1-G+B\  '  .  '  )+6«-(  l-(fcu)  i  0)+4 

4«-(  i  +  p/*(  (irp4)L  (AfCl ]-//),  v)+4)4>4 

-►L4-i~ltL,4«-l  +  p/?«-(0     1    xpp*i?[;  l^[2]««-gL^[2]rQx/AO])+/f 

P«-P*10*2>ll0»|P+0  =  P 

+L3xiO=J'«-+/5'«-(  'fl'eA:)T«0=P^-(L0.  5+.Vx,p)*/i/«-l0*2}«-10i|l-A\G  +  i; 
L4:+(pl+pL)/P-ppA^(l    0    *pG-J  pT\'~  '  )pJ+J  tO*-v /T4-Q>P+b /P 

+  (*L+{0[  LxJ  +  'Z1  eK)[  .  x~T<-(;F  +  <^l  +  [10®ir  \P)>0+L-W -D+0+~2iL)  /  L/ 1  ,  I1  , 

I 

+E  +  xpS<-'  FIELD    WIDTH1 

+L4+1+1  (  (J[2]^Lv.  <0)+O^l  +  10T  .  <|  Z>(£/,L)+7,^10=|P)>y-£)  +  O  + 

3 

71«-~£/«-P[77ilt.7]^L^plpA'«-,£'T  ,,+0"'[-/p2-xZ/],A[l  +  <S?(0pl0)TlL] 
F:+Uv2>£x~'Tletf  )/I,/V^pA^A[ll,l+$(Z;plO)Tr,Vxl|  |  p  ]  ,  X 

D+,  (-2V)  +  (Z?r  .x$X[;2+2?]*l  +  A)<>.<2?«-iZ?-l 

A  ^A/pAM^/xpA"  <-,*>'     ' 
J:^(J^e7vo  =  +  /0^orL-0)/7  +  p2?*-pP^G,  A[l  +  ^(LplO)TL  I  P] 

P^/JpC  ,0«-G<|>0)\(  ,0^0°  .<(  -G)4>\L+G«-l  +  pG)/,P 

^Z?--iJvL^~TL'eA:,P[J7\Z?^l'M^pP^P,^;  ]«-'*' 

P^*p(,4>0)\(  ,0<-~Xt~0)/,P 

+  (~/?)/£-2l2+l,Z^OpP+B^(0,*«-tfxl-2xL)+P 
220:0tf«-(pK«-(~l+2?«-O,  (M[2]Lp0)p0<-(  '  ,  f  =CH  )  I  \  pCH  )  pCH  )$CH 

2?«-,(M[2],;n+    0    "l    +  (M[2],2^p(,4>0°.Si4«-r/2?«-l+0-~l<l>D)U 

->(L0-/Ax5),pOL^OL,  ((l  =  ltAO+Afxi,tf)P2),,P 
EiK+'NO    VALID    E%     7,     02?    P    PHRASE1 

?FM7    PROBLEM  1  9K;  (  l,p5)pS 
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VFILLCDIV 

V  Z+W    FILL    KiNiD;Xl;X2;XZ',X^iYiIiJiS 

[1]  U-pX-rXliX+-,X] 

[2]  Z^0.5x^[L0.25x3+/l/x13]+^[r0.2  5xl+.Vxl3] 

[3]  S«-(Z[l]*0.5  +  1.5x£),  (Z[  1]  *  0.  5  +0«-Z[  3  ]r,Z[  1])  ,  Z  [  1  ] 

[4]  <?^L0.5+5,Z[2],Z[3],  (Z[3]  +  0.5+Z?),Z[3lf0.5  + 

1,5x2? 

[5]  Xl+-U2<£[l])/;r2+-U<.?[2])/XH0.5+X 

[6]  *2«-U2>5[l])/*2 

[71  X4«-U3>.?[7])/X3*-U>5[6]  )/X 

[8]  X3-f(X3<5[7])/^3 

[9]  y+-(  U>S[2]W<5[6]  )/* 

[101  Z«-Vp0 

[111  Z[I>20  +  9[+j4(~l,Xl)°.=I«-\(5*-0rSU)[ll 

[121  Z[I1^10^9[+/("l,^2)o.=I^5[ll+\l+5[2l^5[ll 

[131  Z[I>10  +  9L+*(~l,*3)o.=I+-~l+5[6]  +  ll+5[7]*£[6] 

[14]  Z[I]^20  +  9L^/IUo#=I^"i-r5[7]  +  \l-tVn1?[7] 

[15]  +18xi0=py 

[16]  z[«r+-~ity]+-z[i+-i+y]+-9 

[17]  Z[Ii-i0reT*r  +  l]+-8 

[18]  Z[5[3]  +  i0rc°[5]n5[3]  +  l]^6 

[19]  Z[5[3],5[5]>0 

[20]  Z[£[4]>1 

[211  Z[5[3],5[5]]^Z[1?[3],5[5]]+  2  4 


VIiVPiT[aiV 
V  INPUTiJ ;M;Z 
[1]    'ENTER   NUMBER    OF    TREATMENTS.1 
[2]    ^Ox/OQ 

[3]    'ENTER   IS.CZQ.E.    OF    NUMBER    OF   DBS.    FOR    TREATMENTS    1    TO    '  ;K 
[4]    Y+(  (MAX+[7n+L\)  ,K)pO 
[5]    +7xx(/C=p/V) 

[6]    +4xPpG^' RE-ENTER    COUNT    VECTOR,    ONE    ELEMENT    PER    TREATMEN'. 
[7]    +0xiU<.J+-J  +  i) 

[8]    x  ENTER    ';#[</];'  OBSERVATIONS    FOR    TREATMENT    '  ;«7 
[9]    +7xP  ppy  [;</>(  W/3Z,l)pZ,(.M^W/l^TIpZ^D)p0 
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V    F    KS    G ;  Z  ;  /V  ;  X ;  W  P  ;  Hi  G  ;  A/  4  X ;  5 1  u  :v ;  - 
[1]         ^♦■tf^iw^OxA-f-pZ^ZCiZ^P.Gl 
[2]         -MlM  (A  <  .-<-' i-i ) 
[  3  ]         NF+-NF ,  ( -i-/F<A>Z  [  ■•'  ]  )  *  PF 
C  4]  UF+liFt  (-<■/  t<X)  ipF 

[51         NG+(+/G<X)ipG 

Lb]         SIGiv  +  xivF+JjFo  .  -NG,  ( i-/£<*  )  *pG 
[7]         £>«-£,  #X  a*  xi+(  {NF=MAX+[/iiF+-\*JF)/FIGi1) 

L8]  ->2,/ifff;0 

L9]  /lit  ','*M   F(  a' )-■:;(  a  )    =  »  ;("/£ 
[JO]       '>;J/V    F(A)-r;(A)    =';L/0 


VLIA/F[Q]V 

V  Z+LINE   WiXiYiNiMiP;IiJiKiXliYliX2iY2;FiSLiYI 

[I]  +3xi2=ppV 

[2]  tf«-$(2,pf/)p(ipf/),V 

[3]  F^(A^pA'+-A'[P^A/^[;l]])ptf+-5L^0 

[4]  jri*0,5xlCl^L0,5xM+l]+^[J*r0,5xl+Af^L0,  5+tf*3] 

[5]  A2^0.5xA[/V^lrI]+A[A/  +  ln!/] 

[6]  yi^0,5xZ[Z]  +  (Z^Z[iz*-Wty-Ht/[P;2lrF])[J] 

[7]  J2^0.5xZ[I]^(Z^Z[iz^(-W)ty])[i/] 

[8]  F*-(SWLf(Y2«Yl)*A'2-A'l)xx' 

[9]  a>a>i 

[10]  -»6xXA<2 

[II]  YI+-0.5xY[L0.5xW  +  l]  +  (y+-y[^Y+-(V[P;2]nF])  [fQ.  SxN+1] 
[12]  'SLOP Ei     *;££;'    Y ^INTERCEPT i     »;YZ 

[13]  Z+fili2]»F*-YI+SL*Wlill 

[14]  -►16xlpSPLPf2'  =  i 

[15]  +  0,  ,Z^(2,AOpF,Z 

[16]  Z«-$(2,/l/)prf[;l],Z 


V    R+LIT    A 
[1]  £«-,  »  0123456  7  89'  [li-((H-Ll0»>Opl0)T4l 

V 
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VL5LJ/V£'[G]V 
V  Z+LSLINE   X;N;B 
[1]    +Al*\2=ppX 
[2]    X"-${29N)pX+{\fil-pX),X 

[  3  ]  All  b+X  [  ;  2  ]  S  ( $  (  2  ,  N  )  p  (  ,  X  C  ;  1  ]  )  ,  ( 0+1  +  p  X  )  p  1 ) 
[4]    'SLOPE:     '  ;BCl];'  X -INTERCEPT :     ' ; 6 [ 2  ] 
[5]    Z-$(2,,V)pU[;l]),A:[;2]'(fi[2]+X[;l]x£[l]) 


VMEDPOLISHlUlV 

V  Z-MEDPOLISH   X\R\C \N ;RI \CI \RA\C A;S \RG ;RE \CE \TV 

[I]  i?«-pX[;l] 

[2]  /V-ffxC  +  pXCl;  ] 

[3]  *«-X'in/^0.5xZ[L0.  5*/i/+l]  +  (Z«-ZUZ-,*])[ro.  5*0  +  1] 

[4]  RE+RpO 

[5]  f£«-CpO 

[6]  D-SWI,* 

[7]  Z-U,C)p(,*)U,*  +  (  i/?)o.xCpl*ffG-(r/,^)-L/,^] 

[8]  X+X*(RI4-+/Q.  5*Z[;L0.5*r  +  l]+ZC;ro.5*r  +  l])o.*rpl 

[9]  RE*-RE  +  RI 

[10]  Z*-$(C,ff)p(  ,$X)U,$X  +  (#pl  )°.x(  iC)*AG] 

[II]  W'(ftpl)o.x:i*+/o.  5*Z[L0.5*i?+l;  ]  +  Z[ro.5*i?  +  l;] 
[12]  CE+CE+CI 

[13]  +6*\EPSIL0N<1*  +  /  I  JvS 

[14]  Z<-X 

[15]  +11*\RESIDS*2 

[16]  Z<-$(2,/V)p(  (,flfi'°.*f  £)*?/),  ,Z 

[17]  +  2l*\l\IORMEFFECTS=0 

[18]  RE<-RE*RA<-(+/RE)iR 

[19]  CE<-CE*CA-(+/CE)*C 

[20]  TV+JT+flyl+fd 

[21]  'TYPICAL    VALUE:     '  \TV 

[22]  '/?0tf    EFFECTS:     '  ;i?£ 

[23]  'COLUMN    EFFECTS:     ';CE 
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V    MSTATS   WiZiCiSiC&ihEANS'%VAR8-%NiJ 

LI]  f-*-($Z)  +  *xZ«*-V*(  pW)pMEANS+-(-rtW)  *dNl  +  p  ("2  +  1 , 1  ,  pV  )  pW-fpCH,-*-  '  ' 

[2]  'BF8.2 '    FMT(Q,  \pS)  ,  [  1  ]  (  i  pf )  ,r*5<> .  x£«-(  VARS++J>Z*2)  * 

0.5 

[3]  'LW 4/1/5         ffl,6F8.2'    FAir    AiF/lA/5 

[41  '  DSI-fl.Z?  EV    3,£F8.2»    FAi?    S«-5*  ( /INI )  *  0  .  5 

[5]  Z^t</<-0 

[6]  ■*£'/V^xi(p5)<e/-*-^  +  l 

[7]  r+-X\_iX<-dl;J~n 

[8]  -►6,Z^Z,r-«-o#5xr[L0.2  5x3+/l/xl3l+r[r  0.2  5xl^.Vx!3] 
[9]       EftD\Z+{ (pS) ,3)pZ 

[10]  'fflL.    4#rLSB,flF8.2  '    FAi?    Z[;l] 

[11]  ^MEDIANS    Q,flf8.2'    FA??    Z [ ; 2 ] 

[12]  W.    QifrL^Q.Bf'S  .2  •    FA;?    Z[;3] 


VNUMSUMLUlV 

V  NUMSUM    X;N iQiTliCg 

[I]  Ci? 

[2]  •  NUMERICAL    SUMMARY1 

[3]  70p'n' 

[  4  ]  N+pX+X [  ^*- {X*MISS )/X+-,X] 

[5]  'SAMPLE    SIZE    =     '  ;/V 

[6]  70p»n' 

[7]  Q«-0.5x;aL0.12  5x7+/Vxl7]+*[r0.125xi+A/x17] 

[8]  n*QC4],0,CC«f],0,0,(  (Q[2]+e[6])v2),Q[2],0 

[9]  2,l*-ri,Q[6],(Q[6]nQ[2l),((Q[l]+QC7])*2) 

[10]  ?l^?lie[l],0J5[7],(Q[7]wQ[l]),  ((J[1]+Z[A/])t2) 

[II]  ri*ri,X[l],0,n/in,X[/VW[ll 
[12]  2»1  +    4    5    pTl 

[13]  CH+'MIDPTS,  tLCQ/Q/MlN,  J1EDIAN, ,VPQ/8/MAX3  .SPREADS , ' 

[14]  '£F10.2,ffi     I     ID1    FMT    Tl 

[15]  70p'^' 


VCNEHiQlV 
V    H+-ONEH    X;AiB;C 
[1]         4<-0,X,0 
[2]         B«-*,0,0 

[3]       r+o.o.x 

[4]  H+-    2*{C  +  B->-2xA)  t  4 

[5]  /^[11,2  +  fl.XCp*] 
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V0A/£3[LJ]V 

V  Z+0NE3    X;Yl;Y2iXl;X2iN;C 

[1]  Yl«-(3x*[2])»2xX[3] 

[2]  Y2«-(3xX[,V»l]  )^2x*[~2+AKpA'] 

[3]  X1*-Y19~1+X 

[4]  Z2«-1  +  X,Y2 

[5]  z«-(yi,x,Y2)[(  wv)  +  i  +  (c=xi<;o*U<x2)=r<-xi<x2] 


V  P+PARTIAL    W \C \CIS \M\Y \X\BETA 

[1]    Z^iHpV 

[2]    'ENTER    COLUMNS    Of    MATRIX    TO    BE    PARTIALLED    OUT.' 
[3]  Ali+A2*\Q*  +  /~(C+-U)€.M 
[U]    Y«-tf[;C,JiV-*-(~MeC)/M] 

[5]    +Q,pP-CORRELATION    R*-Y~X  + .*BETA*-Y\&X-1 ,  (2t(py[;f]),l)pi/[;C] 
[6]  A2:+Al,p'UNACCEPTABLE    VALUES,    ENTER    INTEGERS    FROM    1    TO     » ;jtf 
V 


VSHOWRESlQlV 

V  Z+SHOWRES    X\Y  ;N  \Q\D\LW \LF \UW \U F 

[l]  AKpY-YUY-,*] 

[2]  4^0.  5xy[L0.25x3+4Vxi3]+y[ro.25xl  +  Vxi3] 

[3]  iy«-Q[l]-Z?^[3]-(j[l] 

[4]  £/»/«-£[  3  ]  +  *> 

[5]  LF*-Ltf*D*2 

[6]  UF-Ud+D*2 

[7]       y^(i>LF)  +  (A:>Lv)  +  (^>ft[i])  +  (i>^[2])  +  (A'>c[3])  +  (^>(/v)+i>i;i;' 
[8]       z^'ooooooxr ci+y] 

[9]  Z-(  (pZ[;l]),2xP2[l;  ])p(,Z)     4/i/Z>,  (  pZ  )  p  »      ' 
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[1] 

[2] 

[3] 

[4] 

[5] 

[6] 

[7] 

[8] 

L9] 

[10] 

[11] 

[12] 

[13] 

[14] 

[15] 

[16] 

[17] 

[18] 

[19] 

[20] 

[21] 

[22] 

[23] 

[24] 

[25] 

[26] 

[27] 

[28] 

[  29] 

[30] 

[31] 

[32] 

[33] 

[  34] 

[35] 

[36] 

VREGRESSLQ1V 
V  Z+Y    REGRESS   X ;N ;K;C \XPXINV \XPY ;BETA\RSS ;TSS ;S2;ESS \JID\B 

A >  (  2  t  (  p  X.  )  ,  1 )  p  X 

X"-(0,  l*LINIERCEP1)\l,X 

XFXINV+&(§X)+.  *X 

BETA+XFXINV+  .*XPY*-(§X)  +  .  *Y 

RSS+( (§BETA)+.xXPY)-C+( ( ++Y )*2 ) *N-p , Y 

ESS+(TSS^( (§Y)+.*Y)^C)^RSS 

S2«-,555*(iV»l  )»£«-(  p,BETA)^MNTERCEPT 

CR 

'  ANOVA' 

£H+' SOURCE ,DF, SUM    SQUARES, ME  AN    SQU ARE , F~ RATIO1 
i  i 

'^REGRESSI0N\n,in,BE16.n  '    FMT(K),  (,RSS),  (  ,RSS*K)  ,  (  9RSS*K)*S2 

£&+■'  ' 

'ffl      RESIDUALS, IH  ,BE16.H'    FkT{  ( /INI )  »#  ) ,  ( ,  £SS) ,  52 ,  0 

'illTtfTML  fi, 14, 5516. 4  »    F/aT  ( /I/-1 ) ,  ( ,TSS) ,  0  ,  0 

t     » 

'i?  SQUARE:        •  ;  ,RSS*TSS 
'  STD    ERROR:     » ; ,52*0. 5 
£&*•' COEFFICIENTS ,T  STATISTICS1 

'515.4'    FMT$(2,p,S5!Z7i4)p(,557,/0,  (,B5IT/0*(1     1    S/-52  xjfpA;  v/  )  * 
0.  5 
'50    yOi/    f//ltfT    /I    PRINTOUT    OF    THE    V ART. 'AUCE-COV ARIANCE    MATRIX?' 

+ai*\  'y '*i+b 

%VARIAHCE-*COVARIANCE    MATRIX:     ',£5«-,! 

•512.4'    5/rfT    / 
41:  '  DUR  BIU -4  AT  SO  11 1     ,;(+/((14-,C)^('l  +  ,C))*2)i  +  /(  ,C«-Y  '#+.  x5£T>i )  * 

2 

Z^(2,£V)p(,^+.x£5T/J),  ,C 
Bl:'DO    YOU    si  Alii    TO    FORECAST    A    VALUE    FOR    Y?  • 

^Im'/^I+E 

'ENTER    X    VECTOR    (';£;»  VALUES)1 

'FORECAST    OF    Y    VALUE:     »  ;  (C<-(lvbIliTERCEPT)  4-l,G)  +  ,*BETA 

'VARIANCE    OF    FORECAST    ERROR:     '  ;S2*1+C+ .*XPXINV+ ,*§C 

+51 
Clz'UO    YOU    WANT    TO    SCAT    RESIDUALS    VS.    PREDICTED    Y?  ' 

+  0xi  »/v  »  =  ltiL 

DEP-Q.5xiiID*-L/7  0,  (  T  /  (  (  0  .  7  5*.V  )  ,  3  0  )  ) 

5C/U1  Z 
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VSCATlQlV 

V  tf+SCAT   ZiNiXiYiCiR',UiS;L;I;J;K;UT;CLiGiDiB;AiOiV 

I]  -►3*i(2  =  +  /2  =  AOv(x/AO>  +  //V-pZ 
2]  Z^(2,pZ)p(ipZ),Z+-,Z 

3]  y*Z[;l+ir^"l+(pZ)[2]] 

4]  /?«-pZ«-*«-,Zi;;l] 

5]  L-U<-S<-2pd 

6]  J +  1  ■<■  0*p  ( D+N DIV X,N DIV Y  )  ,B+WID,DEP 

71  £'2*«-10*LlO«rL«-lE~2  0f(  (£'[  J  1H  /Z  )  «S[  J  ]H  /Z  )  *£>[«/"  1 

8]  SlJl+-UT*LSU]*UT+UTll  +  MCLnUT+(l    2    5)xUT] 

9]  UlJl-UTxTUlJ]*UT 

10]  L[e/]^l+GxL(fi[J.]wl)TG^(i'[J]«5[.Jl)UT 

II]  Z^.Y 

12]  +7*\3>J+-J  +  1 

13]  4«-($L)pO 

14]  X+l  +  LO.S  +  iLlllnDxiXnSllDiUlllwSLl] 

15]  y«»-l  +  L0.5  +  (I[2]-l)x(y,5[2])*l/C2l«5[2] 

16]  I«-l 

17]  +20*\1<C 

18]  ALXLI  ill  iXilll  +  lOlAlYlI  ill  iXHII  +  l 

19]  +  18  +  6*R<I<*-I  +  1 

20]  J+1 

21]  Z?^-0=7^/i[y[I;t7];Z[Z]] 

22]  4[Y[I;en;*[I]>(10xy>K+l)  +  (  U  +  l )  x  AT  =  l^  )  +  U«-3  5-  2  «J  )  *D 

23]  -»-21x\/teZ«-I  +  l 

24]  ■*21xiC^JP^+I-«-l 

25]  (?^(<t>Pi4)Liri^L0.5^(Lnl)x5i5«i' 

26]  i4[;C[l]]^/3[;6)[l]]  +  36xO  =  >3  [;(?[!]] 

27]  AlCl!]  11+AIOL2]  i1-r3S*0  =  ALCL2]  il 

28]  J/*-'     o23i\56789miiKj(JlIlH!iGGFlEEDDCCBB.AAw\'[l  +  9A> 


29 
30 


'RANGE    OF    Xt     »  ;S[  1]  ,£/[  1] 
'RANGE   OF   Y%     »  ;S[  2]  ,£'[  21 


V 
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[11 

[2] 
[3] 
[4] 

[5] 

[61 

[71 

[8] 

[91 

[101 

[111 

[121 


VSPL2T[Q1V 
V  Z+-SPLIT    XiPiQiRiSiTiNiliJ  iC 


3+X 
'2+1+X 
'l  +  2+X 

■3+a: 


S+3+X 

I+(((P<Q)*R>S)^(P>Q)*R<S)x 

G+(X+X[.l~[tX,XlN'\)lH-I+-(I*0 

R+{Z*Q)m2*XLl] 

P+XH+21 

T+(3*S+XlI  +  H'])»2xXlI  +  5'\ 

I[I  +  2]^(i?xl=J-)i.(Px"l=J)  + 

X[It3]*(imxl=8f)t(Px"lsJ)  + 

Z<-l+~l+X 


-Zx\~3f,V«-p;n 


QxQ=J+(C=P<Q)n(Q<R)=C+P<k 
*xO=J+(C=P<S)  *(S*T)=C-*-P<T 


[11 

[21 

[31 

[41 

[51 

[61 

[71 

[81 

[91 

[101 

[111 

[121 

[131 

[141 

[151 

[161 

[171 

[181 

[191 


VSTATISTICSLUlV 

STATISTICS    R;XBAR;VAR;N ;Q;S ;Z 

XBAR+(+/R)iN+-p  ,ff 

VAR-{+/{ (R*XBAR)*2) )iN«l 

S+RlkR~\ 

Q-Q.S*SlL0.2S*3->-Nx\3~\+Sl[0 
A/£Mi  '  ;XBAR 
VARIANCE:     *;VAR 
STD.    DEV.i     ' iS+VAR*0.5 
COEFF,    CF    VARIATION \     ' 
L0V£/?  QVARTILEx     '  ;Q[ll 
£'PP£tf  QUARTILE:     »  ;G[3l 
MEDIAN  i     ' 
T RIME  AN \ 
MIDMEANi 
RANGE i     ' ; 
MIDRANGEx 


C  i 


2  5xl+.Vx  i3l 


*2J/1A 


Q[3]+2xQ[2] 

( (P>S[ll)A(A<g[3l) )/A 


;  G  [  2  1 

1  ;0.25*Q[ll 

'  ;  (fpZ)x  +  /Z 

cr/p)"L/p 

';(  (r/tf)fL//?)*2 

M£/]A/  ABSOLUTE    DE.IATICNt     »  ;  (  +/  |  fl*Q[  2  1  )  *tf 
INTERQUARTILE    RANGE:     f ;Q[ 3] "Q[ 11 

COEFF.    OF    SKEWNESSi     '  ;  (  (  +/  (  RnXBAR  )  *  3  )  *0«1 )  *S*  3 
COEFF.    CF    KURTOSISi     '  ;  (  (  ( ■►/  (  RwXBAR)  *4  )  frtf  -1  )  *5*u  ) 
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VSTEMLEAFlU~\V 
V  Z+STEMLEAF   X ;C ;R ;S ;SI ;I ;J ;F ;A ;V ;L ;AW ;AA  ;XW 

[1]  r4rlO*l«lld9R+-lE~  20+{XlpX]*(X+-XllX'-tXl)  111)  *$C  ALE 

[2]  5.Z>-(l  +  3x  +  /(i?xf  )>    25    50)W(p*)>    25    100 

[3]  XH0.5+#xf  xiO*  +  /SZ=     2    3    6 

[4]  F*+/(5I=i9)/    0.5    2    11    0.5    2    2    1    0.5 

[5]  A+-20p'012Z*S678$ABCDEFGHIJi  ,Z«-'  ' 

[6]  I*-FxUCl]*10xF 

[7]  Xrf+8+VW+WIDTH** 

[8]  ^'«     '[l+r/^[n>0] 

[9]  2>l  +  (  |V-H(/s(10xI)^x"i  +  i0xF)/^)7t10xL|I 

[10]  /4/K/Ul+(10    10)t|LI] 

[111  +((pAlLl)*WW)/Ll 

[12]  Z-Z%XWpS%AA,  f  |  f  .  OT  +  4CL1  )  ,  »  +  '  ,  (LIlr'(p^[L])r,f/y),2p' 

[13]  +L2 

[14]  LliZ«-Z  .K/p^./l/?,  •  1  '  ,^[L],^p»  » 

[15]  L2tZ^-I  +  Fxi^(J=0)xZ[l]<0 

[16]  -*8xiO<pX*-(pf/)+X 

[17]  Z«-(  (L(pZ)*XJ/),;W)pZ 


V5i/M5Q[D]7 

V  SUMSQiC iNUMBERiB ;BLDF 

[I]  7*"t-/Y 

[2]  TSS+(+/(+/Y*2)  ).^(  (+/T)*2)*NUMBER++/N 

[3]  BLSS+i+Zi  UB+  +  /Y)*2)iK)  )rf 

[4]  MSBL-BLSS*BLDF+MAX-1 

[5]  V7xi(0P2T00=2) 

[6]  BLSS+BLDF+Q 

[7]  WTSS+-TSS*BLSS  +  BTSS<-(+/  (  (7*2) *tf)  )-r 

[8]  MSR+BTSSiK^l 

[9]  TDF+NL'MBERnl 

[10]  F-MSR*MSE'-tfTSS*EDF'-TDFrBLDF'>-K-l 

[II]  BLCCK+-(B*K)»YBAR+{+/T)*NUMBER 
[12]  '  M<?7i4    TABLE1 

[13]  '       SOURCE  DF  SB  MS  F' 

[14]  Cflo-11 

[15]  'fl      T/ff^rWffA/rH.lS.FlS^.Fll^^S^1    FMT(Knl)  ,BTSS,MSR,F 

[16]  -*{0PTICN-1)  /Lb 

[17]  'CI   #L?rA:F(I],I8,F13.2,F11.2,F8.2  »  FMT  {MAX-D  ,BLSS  ,MSBL,{MSBL*MEE 

) 
[18]  Z^t'E   FFFflFffl.ig.FlS^.Fll^'  FM2*  EDF,WTSS,MSE 

[19]  'ffl   r^iTi4Ll!l.r9,F13. 2 »  FMI7  TDF,TSS 

[20]  *M«SQUARE    =    ffl,F5.3'  FMT (BTSS+BLSS) iTSS 

[21]  'D0I/FF/5LL  MFM  =  E.F10.21  FMT  YS/3F 

[22]  ^TREATMENT    EFFECTS   fl,F6.2*    FMT (T*N) *YBAR 

[23]  +(0PTI0N=1) /0 

[24]  'HBLCCX  EFFECTS    ffl,F6.2f  FM2"  flLCT/( 
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[1] 


t-TXICE    X;Y 
Z*-Y+A3RSH    X- 


Y*-A3kSh    X 


V!JTCOUDl[]~]SI 

V  UTCOND    X 

[I]  X«-(2p(pX),l)pX 
[2]  Z-» » 

[3]  l^0*/V/r*-(pX)[2] 

[  4  ]  N+pU<-U  LiU^iU  *MISS  )  /U-X  [  ;  7+  1 J  ] 

[5]  3^-0  .  5  xi/  [  L  0  . 1 2  5x  7  +/i/x  i  7  ]  +i/  [  T  0  . 1 2  5x  1 +/V  x  i  7  ] 

[6]  +7+2xlNUM*2 

[7]  Z*-Z,fi[4] 

[8]  ^4  +  llxfl/i?=J^J  +  i 

[10]  -»-ifx^it>J^J+l 

[II]  Z-Z  ,  £/  [  1  ]  ,  Q  [  2  ]  ,  2  [  4  ]  ,  2  [  6  ]  ,  £/  [  N  ]  ,  V 
[12]  >'+  +  13x/i/i?=J-»-J+l 

[13]  Z*Z,tf [i],«[l     2     4     6     7],£/[/n,fl 

[14]  ->4  +  16x/VJt  =  7«-J+l 

[15]  QJi-%  MEDIANS' 

[16]  +  0xpU^'F10.2  '     FWr    Z-(/i//?,  l)pZ 

[17]  Z«-(/V#,6)pZ 

[18]  Q£<-tMIN,L0Q9MEDIAN,UPii9i4AX9SIZE  ' 

[19]  ->0xpu^'Fl0.2  »    F^7    Z 

[20]  Z<-UA,8)pZ 

[  21  ]  £H4-'MIN,L08,L0Q,  MEDIAN,  UPQ ,UPB ,MAX 9SIZE  ' 

[22]  D*'F9.2'  FMT  Z 
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